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摘 要 


玻 色 - 爱 因 斯 坦 凝聚 (BEC) 是 一 类 涉及 物理 学 的 很 多 领域 
的 普遍 物理 现象 1924 年 爱 因 斯 坦 预言 了 当 温 度 低 于 临界 温度 
时 服从 玻 色 统计 的 理想 中 性 原子 气体 的 凝聚 现象 . 1995 ERE 
属 原子 ”Rb ,”Na 和 ?Li 通过 激光 冷却 、 静 磁 阱 与 蒸发 冷却 等 
技术 实现 了 近 理想 气体 的 BEC; 1998 年 实现 了 和 氢 原子 的 BEC. 这 
种 上 万 乃至 上 百 万 原子 构成 的 宏观 量子 系统 ， 以 相干 的 方式 在 
WE. 由 于 长 的 德 布 罗 意 波长 ， 原 子 完全 失去 了 孤立 粒子 的 特 
征 ，BEC 将 量子 现象 带 到 了 宏观 尺度 . 近 几 年 这 种 新 物 态 特 性 
的 研究 有 广泛 和 迅猛 的 进展 

在 阅读 了 许多 有 关 BEC 的 理论 和 实验 文献 的 基础 上 ,我 们 
主要 在 理论 上 做 了 一 些 工作 . @ 求 得 了 一 维 箱 势 阱 中 有 排斥 和 
吸引 相互 作用 情形 中 性 原子 BEC 波 函数 和 能 量 本 征 值 的 解析 解 ， 
相对 于 无 相互 作用 体系 ， 排 斥 体系 的 能 级 升 高 了 ， 原 子 分 布 趋 
于 分 散 ; 吸引 体系 的 能 级 降低 了 ， 原 子 分 布 趋 于 集中 ，@ 给 出 
了 求解 有 原子 间 相 互 作用 时 中 性 原子 BEC. 的 方法 ,估算 了 这 种 
非 理 想 气 体 BEC 的 临界 温度 随 原子 间 相互 作用 大 小 和 原子 密度 
的 变化 情况 ， 得 出 排斥 相互 作用 增 大 , 将 使 得 激发 态 能 容纳 的 
原子 数 减少 ， 临 界 温度 升 高 . @ 给 出 了 一 种 非 线 性 定 态 薛 定 谓 
方程 (NLSE) 数 值 求解 方法 ， 即 对 非 线性 微分 方程 进行 数值 积分 
上 时， 调整 波 函 数 的 边界 值 B (0) 和 本 征 值 8， 使 B(x) E x ->w 时 
收敛 于 零 和 归 一 化 条 件 都 得 到 满足 . 并 用 这 种 方法 求解 了 简 谐 


POOR AA IO AE ect ee net 
势 阱 中 中 性 原子 NLSE 的 基态 和 激发 态 的 能 量 本 征 值 B, MR 
数 . 讨论 了 参数 四 (0) ， 的 确定 与 积分 收 化 及 波 函 数 归 一 化 的 
KAR. 四 三 维 球 方 势 阱 和 简 谐 势 阱 中 具有 吸引 相互 作用 原子 体 
系 NLSE 的 数值 本 征 解 其 特性 是 类 似 的 ， 即 原子 的 空间 分 布 随 
相互 作用 的 增强 而 被 连续 压缩， 以 及 凝聚 原子 数 随 能 量 本 征 什 
变化 的 双 稳 态 曲 线 . 临界 点 将 本 征 态 分 为 亚 稳 的 凝聚 态 和 不 稳 
的 稠密 态 ， 并 给 出 与 实验 报道 相符 的 吸引 BEC 所 能 包含 的 最 大 
原子 数 Ne. 在 此 基础 上 我 们 又 给 出 了 - -个 严格 可 解 的 宏观 量子 
隧道 模型 ,计算 了 从 亚 稳 凝 聚 态 到 塌 缩 态 的 宏观 量子 隧道 速率 
分 析 了 吸引 BEC 的 稳定 性 . @ 根 据 最 近 的 BEC 的 动力 学 模型 和 
主 方程 ， 用 母 函 数 方法 求 得 玻 色 凝聚 体 主 方程 的 解 ， 并 计算 了 
凝聚 体形 成 的 增 率 与 起 伏 ， 增 率 曲 线 存在 一 平台 ， 而 平台 的 粒 
子 数 起 伏 则 是 超 泊 松 分 布 的 . 初始 增 率 7 随 最 终 平衡 原子 数 eq 
的 增加 而 逐步 减 小 ， 与 实验 结果 大 体 上 -- 致 . @ 给 出 N 原子 辐 
射 几率 的 准确 求解 模型 ， 求 得 BEC. 自发 与 受 激 跃迁 几率 的 增强 
因子 正比 于 用 27) (A 为 原子 的 跃迁 波长 ) 体 积 内 的 原子 数 ，， 
并 依赖 于 德 布 罗 意 波长 如. 当 德 布 罗 意 波长 2 与 原子 的 跃迁 波 
长 4 之 比 b8 = V2x hs 14=2.45 时, 辐射 速率 的 增强 因子 y' 达到 
它 的 极 大 值 /'=1.453. 增强 因子 y' 超 过 极 大 之 后 随 温度 的 进 一 
步 降低 而 减 小 ， 当 kT7 0. Ey oth, HERR Y EYE. 另 
一 方面 ， 当 8 很 大 时 、 能 级 的 Lamb 移 位 的 修正 5 趋 于 一 个 有 
限 值 -1.0. 

Xu GE(C-XDPEDHURSE. JERE E AS AE. 原子 激光 ， 
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Abstract 


Bose-Einstein condensation (BEC) is a ubiquitous phenomenon relating to 
many fields in physics. In 1924, Einstein predicted the phenomenon of 
condensation for an ideal gas of neutral atoms obeyed Bose-statistics under the 
temperature below critical value. Up to 1995, this idea were demonstrated in 
experiment of near-perfect gas of alkali atom % Rb ,? Na and "Li by 
loading laser-cooled atoms into a magnetic trap and successive cooling by 
evaporating, the BEC of atomic hydrogen was observed in 1998. A million or 
so atoms in the macroscopic quantum system are all evolving in a coherent way. 
They have completely lost their individual identity. The quantum phenomenon 
is brought in macroscopic scale by BEC. The extensive studies on the properties 
of a new state of matter have rapid developments in recent years. Based on a lot 
ofthe references relating to the experiment and the theory of BEC we have done 
some theoretical work in the following respects. (1) An analytic solution of the 
energy eigenvalue and wave function of neutral atoms with repellent and 
attractive interactions in a box potential are obtained. Compared with the System 
without interactions, with increase of interactions or number of atom, the energy 
levels of repellent system increases and the distributions of atoms approach 
decentralization, whereas the energy levels of attractive system decreases and 
the distributions of atoms approach centeralization. @ A general approach to 
the BEC of neutral atoms with interaction is presented. The critical temperature 
of imperfect gas BEC is estimated. With the increasing of interaction the 
number atom hold in exited states decreases and the critical temperature 
increases. (3) We present a numerical method to solve the stationary nonlinear 
Schrodinger equation (NLSE) with an extemal potential. by adjusting 
eigenvalue # and boundary value of wave function ?' (0) . such that the 
conditions of convergence and normalization of the wave function being 
satisfied simultaneously. We discuss the questions of convergence and 
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normalization of NLSE wave function. The accuracy of calculation is analyzed. 
@ The numerical solutions of neutral atoms with attractive interaction in a 
harmonic trap and a spherically symmetric well potential of depth infinite are 
calculated. And the characteristics of two solutions are analogous on the 
compressed atom spatial distributions and the bistability curves of condensate 
atoms versus the energy eigenvalues. The critical point divides the eigenstates 
into metastable condensate and dense state, and gives the maximum number of 
atoms with which an attractive BEC can hold, and which is in agreement with 
the experiment on the whole. Based on this we present an exactly solvable 
model for macroscopic quantum tunneling of BEC with attractive interaction, 
and calculate the rate of macroscopic quantum tunneling from a metastable 
condensate state to the collapse state and analyze the stability of the attractive 
BEC. (According to recent quantum kinetic theory and using the generating 
function method to solve the master equation of BEC, we evaluate the growth 
rate, statistical fluctuation of condensate atoms. It is obtained that there is a 


plateau in the growth rate curve and a super-Poissonian distribution of statistical 
fluctuation. It is observed that the growth rate y decreases progressively with 


the increase of the final equilibrium value of condensate 7, which agrees 
with the experimental result on the whole. (6) An exactly solved model for the 
emission by N atoms is presented, the spontaneous and induced transition 
rates are enhanced by a factor which is proportional to the number of atoms n 
in the volume 2? (2x?) (A is the transition wavelength of atom) and 
dependent on the de-Broglie wavelength Ag in a more complicated way. The 
enhancement factor y' attains its maximum »'- 1.453 at £y = J2x An /A4 
= 2.45 , then decays to zero for large £y ; and the energy level Lamb shift 
factor Ó' approachesto —1.0 forlarge £g. 


Key words Bose-Einstein condensation. nonlinear Schrödinger equation. 


atomlaser. optical property 


玻 色 - 爱 因 斯 坦 凝 聚 体形 成 的 动力 学 和 光学 性 质 研究 


H x 
绪论 1 
BEC 研究 的 简单 回顾 nnn ; 
实现 BEC HUF DOR n n 3 
碱 金属 原子 的 BEC SR ee 6 
BEC 研究 的 理论 模型 e MM 14 


BEC 的 应 用 前 景 及 研究 意义 …… 16 


本 文 的 安排 与 我 们 的 工作 m 17 
理想 气体 的 玻 色 - 爱 因 斯 坦 凝 聚 … sr 19 
均匀 玻 色 气体 的 BEC nm M 19 
简 谐 势 阱 中 的 BEC Mee 22 
简 谐 势 阱 中 坡 色 气体 在 有 限 温度 的 BEC ………… 25 


有 限 原 子 数 对 BEC 的 影响 eM 28 
as : [8] HE Xr BEC 的 影 A 30 
TP RT VEE PS GE AK] BEC ee MM 33 
排斥 相互 作用 对 BEC BR sic etaed 35 
引言 言 35 
近 理 想 Bose 气体 的 理论 描写 ee 35 
Bogoliubov 近似 微 扰 法 与 BEC -一 般 方法 ………… 39 
有 排斥 相互 作用 的 中 性 原子 的 态 函 数 …………… 4l 
中 性 Bose 原子 的 BEC 计算 … EREN 


与 局 域 密度 近似 模型 的 比较 
本 cR 


简 谐 势 阱 中 有 排斥 相互 作用 中 性 原子 的 BEC. ………: 55 


引 A LI 55 
JEER PE OE PS JT FE BE BOR ROI 


2000 年 上 海 大 学 博士 学 位 论文 


45 与 大 NM 极限 下 基态 能 量 的 比较 coerce ere terre eee 68 
4.6 小 a2 Me heh enhn hen nter e 71 
SHS ”吸引 相互 作用 对 BEC 的 影响 …… 72 
5S1 引 言 es 72 
52 负 散 射 长 度 的 计算 与 测量 n ees eee eee eee tenons 73 
5.3 “一 维 箱 势 中 有 吸引 作用 体系 的 能 级 与 态 函 数 …… 74 
54 ”三维 箱 势 中 有 吸引 作用 体系 的 基态 解 ………… 79 
55 有 吸引 相互 作用 的 中 性 原子 的 BEC 计算 een 82 
5.6 A PERS ME E AES 84 
5.7 吸引 性 BEC 的 稳定 性 分 析 TTE 88 
58 宏观 量子 隧道 速率 的 计算 … 90 
5.9 跃迁 速率 的 计算 IE e 9] 
5.10 小 结 MH e. 92 
BAR ”BEC 形成 的 动力 学 ee e 94 
6.1 引 A rrr rere ale de sere 4 
62 BEC 的 动力 学 模型 …………… .95 
63 ”占有 数 p 满 足 的 主 方程 … - 100 
64 主 方程 的 解 ee RI 101 
65 ”凝聚 体 原子 数 的 增益 与 统计 起 伏 的 计算 ven 106 
66 凝聚 增 率 与 平衡 原子 数 的 关系 …- … 109 
67 AN. GE ee nenne eorr nnn 111 
第 七 章 BEC 自发 与 受 激 跃迁 几率 的 增强 . 112 
Jo = esten iei cgaei iem 112 
7.2 单 原子 重 射 的 线 宽 与 能 级 移 位 …… 114 
73 NN 原子 的 自发 跃迁 速率 enm 117 
7.4 NW 原子 的 受 激 跃 迁 速率 nnn 124 
75 数值 结果 及 与 文献 结果 的 比较 …… .. 126 
7.7 小 缚 oH MIR - 128 
FIE ”结论 和 展望 - 129 
参考 文献 TT IR 132 
Jm ist HELD 143 


Ji €, XE VÀ JE HL ae A V JE n 0 20 7] A RC HE EE 


第 一 章 绪论 


玻 色 - 爱 因 斯 坦 凝 聚 实际 是 -类 涉及 物理 学 的 很 多 领域 的 普 
遍 的 物理 现象 ， 包 括 原子 分 子 物理 、 量 子 光学 、 统 计 物 理 和 凝 
聚 态 物理 等 . 1995 年 美国 的 JILA 和 MIT 小 组 分 别 用 碱 金属 原子 
”Rb A Na 通过 激光 冷却 、 静 磁 阱 与 蒸发 冷却 等 技术 实现 了 
BEC'?, Rice 大 学 的 研究 小 组 也 报道 了 超 冷 7Li 原子 进入 量子 
简 并 区 域 的 实验 现象 M，1998 年 实现 了 和 氢 原子 的 BECU. 这 种 上 
万 旋 至 上 百 万 原子 构成 的 宏观 量子 系统 ， 以 相干 的 方式 在 演变 . 
由 于 长 的 德 布 罗 意 波长 ， 原 子 完全 失去 了 孤立 粒子 的 特征 ，BEC 
将 量子 现象 带 到 了 宏观 尺度 . 近 几 年 这 种 新 物 态 特性 的 研究 有 
广泛 和 迅猛 的 进展 . 这 一 物理 现象 的 研究 会 对 将 来 的 物理 学 产 
生 深 远 的 影响 . 


1.1 BEC 研究 的 简单 回顾 加 


1924 年 6 月 24 H, 30 岁 的 印度 物理 教师 玻 色 (S.Bose) 寄 送 
- 份 手稿 给 爱 因 斯 坦 ， 并 附 信 说 : “尊敬 的 先生 ， 我 冒昧 地 送 上 
所 附 手稿 请 您 审阅 . 我 迫切 想 听 到 您 的 意见 . 我 试图 不 依赖 经 
典 电动 力学 来 推导 普 朗 克 ( 黑 体 辐射 ) 定 律 的 系数 8 x? 17 ， 办 法 
是 假定 相 空 间 最 基本 区 域 的 体积 为 志 ， 我 的 德 文 不 足以 翻译 此 
稿 . 如 果 您 认为 此 文 值得 发 表 并 安排 在 〈( 德 国 ) 物 理 杂 志 》 上 ， 
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我 将 不 胜 感激 "91 爱 因 斯 坦 亲自 把 手稿 译 成 德 文 ， 送 去 发 表 [， 
并 在 文 末 加 注 说 : “我 以 为 对 普 朗 克 公 式 的 推导 乃 是 一 项 重大 进 
步 . 所 用 方法 也 将 导致 理想 气体 的 量子 理论 ， 我 将 另行 前 述 .” 
爱 因 斯 坦 接着 发 表 了 文章 Bl， 其 中 第 二 篇 里 提出 了 原子 将 “ 凝 
K” 的 大 胆 预言 . 

原来 爱 因 斯 坦 在 1905 年 关于 光电 效应 的 著名 论文 中 ， 第 一 
次 把 普 朗 克 的 “能 量 量子 ”概念 推广 到 光 ， 虽 然 光子 一 词 直 到 
1926 年 才 出 现 . 玻 色 1924 年 推导 的 关键 , 是 用 光子 状态 的 计数 ， 
而 不 是 光子 计数 . 爱 因 斯 坦 立 即 把 状态 计数 的 思想 用 到 组 成 理 
想 气体 的 原子 ， 指 出 遵从 这 种 统计 的 气体 将 在 一 定 的 转变 温度 
下 发 生 凝 聚 : 部 分 原子 将 落 人 动能 为 零 的 最 低 量 子 态 ， 其 他 原 
子 则 组 成 “饱和 理想 气体 当时 两 种 量子 统计 的 存在 、 统 计 同 
多 粒子 波 函 数 在 粒子 交换 下 对 称 性 的 关系 ， 以 及 对 称 和 自 旋 的 
关系 等 等 ， 都 尚未 认识 . 难怪 爱 因 斯 坦 曾 把 电子 气体 也 列 人 可 
能 观察 此 种 凝聚 的 对 象 之 一 . 当时 相 变 也 没有 被 作为 一 类 深刻 
的 物理 问题 ， 人 们 自然 未 意识 到 BEC 乃 是 首次 由 统计 物理 学 推 
出 的 一 种 相 变 现象 . 

事实 上 人们 也 不 知道 到 哪里 去 寻求 这 类 凝聚 现象 . 直到 30 
年 代 ， 伦 敦 (F. London) 才 指 出 ， 超 流 和 超 导 现 象 可 能 是 BEC 
的 表现 . 不 过 这 两 类 现象 都 发 生 在 强 相 互 作用 的 体系 中 . 后 来 才 
知道 ， 超 流 液 氨 中 只 有 约 10% 的 原子 凝聚 . 超 导 与 BEC 的 关系 
要 经 过 电子 配对 ， 涉 及 更 复杂 的 相互 作用 . 由 于 粒子 间 强 相互 
作用 而 导致 BEC 相 变 的 纯 量 子 统计 特性 复杂 化 . 很 难看 成 是 纯 
的 BEC. 这 种 情况 促使 人 们 去 寻求 在 稀薄 的 、 弱 相互 作用 下 的 
玻 色 气体 中 实现 BEC. 1959 年 有 人 提出 自 旋 极 化 氢 原 子 气体 可 
能 是 BEC 的 候选 者 9 1980 年 ， 另 -… 种 候选 者 一 一 氧化 亚 铜 中 
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的 激 子 被 提出 WY 经 过 10 多 年 的 努力 ，1993 年 在 实验 上 观察 到 
了 激 子 的 BEC LR. 可 是 由 于 复杂 的 相互 作用 过 程 ，BEC f$ 
性 仍然 得 不 到 很 好 的 研究 ，80 年 代 中 期 ， 激 光 冷 却 和 捕 陷 原 子 
的 研究 已 得 到 长 足 的 进步 ， 几 个 研究 小 组 提出 了 冷却 的 碱 金属 
原子 可 以 形成 只 有 很 弱 相 互 作用 的 BEC. 在 不 断 解 决 实现 BEC 
的 一 系列 技术 难题 后 ，1995 年 碱 金 属 原子 的 BEC 终于 被 实验 
实现 了 . 1998 年 自 旋 极 化 氢 原 子 气体 的 BEC 也 被 观察 到 中 .这 使 
人 们 对 量子 简 并 气体 的 BEC 研究 进入 了 一 个 全 新 的 阶段 . 


1.2 实现 BEC 的 相关 技术 


1.2.1 激光 冷却 和 央 禁 技术 ! 


根据 共振 光 偏 转 原子 束 的 原理 ， 当 光子 被 原子 吸收 后 ， 便 
获得 动量 站 ， 然 后 以 自发 辐射 或 受 激 辐射 形式 释放 光子 ， 回 到 
基态 . 如 果 以 自发 辐射 形式 辐射 出 去 ， 其 辐射 方向 在 4 1 立体 角 
内 均匀 分 布 ， 给 予 原子 的 平均 动量 为 零 ， 故 原子 获得 所 吸收 光 
TAE RE. 只 要 激光 频率 o, 相对 于 原子 跃迁 的 频率 o, 为 红 
Feo, -wo <0， 在 激光 作用 下 的 原子 ， 会 由 于 Doppler 效应 被 
减速 ， 被 致 冷 ， 即 v,v? 随 作用 时 间 1 的 增加 而 下 降 . 但 这 只 是 受 
到 迎面 而 来 的 激光 辐射 压力 作用 的 结果 . 如 果 是 受到 正 反方 向 
传播 的 激光 作用 ， 这 原子 将 在 传播 方向 例如 x 轴 方 向 被 减速 如 
果 是 采用 6 ERG x yv. = 正 反 方向 作用 于 原子 ， 则 原子 将 被 禁 
MITE 6 束 光 作用 的 小 区 域内 ， 又 考虑 到 原子 的 速度 扩散 ， 原 子 
很 像 是 在 一 带 有 粘性 的 液体 即 光学 粘 胶 (Optical molasses) 中 运动 . 
另外 ， 当 激光 的 偏振 随 空间 坐标 变化 ， 亦 即 存在 偏振 梯度 时 ， 
对 与 之 相互 作用 的 原子 呈现 出 阻力 ， 而 且 这 时 的 阻尼 系数 当 原 
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子 的 速度 v 一 0 时 ， JLF 通过 这 种 激光 冷却 和 内 


禁 技 术 可 获得 大 量 的 高 密度 的 超 冷 ， 数 日 为 ~10" ， 密 度 为 
10! -102 em, 温度 为 几 十 EK. 


1.2.2 静 和 磁 阱 技术 0 

按照 麦克 斯 韦 的 电磁 理论 ， 具 有 固有 磁 矩 的 粒子 可 以 被 办 
禁 在 磁场 的 极 小 处 ， 不 均匀 磁场 下 对 磁 偶 极 矩 w 的 作用 力 为 : 
Fy --V(u- B). 可 用 静 磁 阱 来 囚禁 具有 磁 偶 极 矩 的 原子 . 结构 
最 简单 的 静 磁 阱 是 反 向 雍 姆 霍 兹 线圈 形成 的 四 极 阱 ， 其 中 心 处 
Sik WA. 当 线圈 间距 等 于 其 半径 的 125 倍 时 ， 沿 径 向 的 阱 深 
与 沿 轴 向 的 阱 深 相同 ， 其 特点 为 因 禁 力 强 和 因 禁 时 间 长 ,由 
于 这 种 静 磁 阱 的 磁场 强度 在 磁场 中 心 处 为 替 ， 在 其 附近 磁场 方 
向 变化 剧烈 ， 原 子 经 过 磁场 零点 时 ， 它 的 自 旋 取向 可 能 发 生变 
化 ， 对 于 自 旋 取向 反 转 的 原子 ， 势 场 零点 为 势能 最 高 点 ， 原 子 
会 逸 出 阱 外 . 由 于 存在 上 述 的 “漏洞 "， 严 重地 限制 了 阱 中 原子 
密度 的 增加 .Cornell 等 人 利用 时 间 旋 转 势 阱 (time orbiting 
potential，TOP) 来 堵 些 “漏洞” 中， 即 在 一 个 大 的 球 四 极 场 的 基 
础 上 状 加 一 个 以 凡 MHz 旋转 的 小 的 横向 磁场 而 形成 的 磁 阱 . NE 
转 频率 的 选择 原则 是 : 旋转 频率 足够 低 ， 使 原子 自 旋 取向 能 够 
缓慢 地 跟随 它 所 在 点 的 磁场 方向 变化 ; 但 对 原子 的 空间 运动 来 
说 ， 旋 转 频 率 又 足够 高 ， 以 致 于 原子 在 空间 的 运动 主要 由 势能 
的 时 间 平 均 来 决定 . 加 上 小 的 磁场 后 ， 使 四 极 场 势 阱 随时 间 在 
空间 旋转 ， 其 时 间 平 均 为 -- 个 椭 球 面 简 谐 势 阱 ， 势 阱 最 小 点 附 
近 非 零 ， 且 变化 平滑 . 因而 自 旋 取向 反 转 的 问题 大 大 减少 . MIT 
的 Ketterle 等 人 则 利用 聚焦 于 磁场 零点 的 高 斯 光束 所 产生 的 排 
斥 势 来 堵 此 漏洞 中 ， 这 个 光学 寒 由 一- 束 Ar’ 激光 束 (514nm) 形 成 ， 
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其 功 浇 为 3.5mW， 束 腰 为 304 m， 它 在 原点 使 ”Na 原子 的 基态 
产生 7MHz 的 光 频 移 势 ， 原 子 会 受到 排斥 力 的 作用 而 无 法 进入 
磁场 零点 附近 (激光 强度 最 大 ) 的 区 域 ， 于 是 极 大 地 抑制 了 原 
了 通过 漏洞 逃逸 出 去 的 几率 . 同时 ， 由 于 所 用 的 激光 频率 远离 
钠 原 子 的 共振 频率 ， 克 服 了 由 共振 光子 散射 而 引起 的 加 热效应 . 

另外 还 有 ， 用 永久 磁铁 以 “loffe” 组 态 形成 的 静 磁 阱 49 和 
利用 直流 电磁 铁 组 成 的 静 磁 阱 即 “Cloverleaf " 阱 (71， 可 有 效 地 解 
决 普通 四 极 磁 阱 中 心 存在 的 “原子 泄漏 ”问题 . 
1.2.3 蒸发 冷却 技术 9 

蒸发 冷却 是 有 选择 地 把 磁 阱 中 能 量 较 高 的 原子 释放 出 来 ， 
然后 剩 下 的 原子 通过 弹性 碰撞 重新 达到 温度 更 低 的 热平衡 ， 如 
此 反复 不 断 降低 原子 气体 的 温度 . 在 实现 BEC WE, RH 
冷却 是 由 一 个 射频 磁场 来 完成 的 9 在 磁 阱 中 ， 能 量 较 大 的 原 
子 可 达到 磁场 较 强 的 地 方 ， 产 生 的 塞 曼 分 裂 也 较 大 . 可 选择 射 
频 场 频率 ， 使 这 些 原子 从 因 禁 自 旋 态 ( 例 m=+1 ) 跃迁 到 非 办 
禁 的 自 旋 态 ( m+1 ) 而 逸 出 磁 阱 ， 通 过 把 射频 场 频率 慢 慢 恋 
低 ， 迫 使 更 多 的 能 量 较 高 的 原子 逸 出 磁 阱 . 于 是 ， 阱 中 的 原子 
密度 和 弹性 碰撞 儿 率 增加 ， 温 度 变 低 ， 最 终 的 温度 和 相 空间 密 
度 取决 于 最 后 的 射频 场 频率 
1.24 BEC 检测 技术 [3 

日 前 观测 BEC 的 形成 多 采用 共振 吸收 成 象 技术 (YY， 用 这 种 
技术 可 以 确定 原子 的 数 日 、 密 度 、 温 度 以 及 原子 的 空间 分 布 ， 其 
具体 过 程 为 : 突然 关闭 势 竺 ， 让 发 生 凝 聚 的 原子 云 白 由 扩散 ， 
然后 在 不 同 的 延迟 时 刻 用 共振 脉冲 光 来 探测 . 由 于 原子 对 共振 
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光 的 吸收 ， 企 探测 光 中 会 产生 阴影 区 ， 由 CCD 装置 对 透射 光 
成 象 ， 对 图 象 进行 数值 化 处 理 ， 可 以 得 到 原子 云 在 每 一 点 的 光 
学 厚度 . 对 由 此 获得 的 一 系列 飞行 时 间 (time-of-flight, TOF) 
象 进行 遂 点 校正 ， 以 修正 由 探测 光 的 偏振 度 和 饱和 效应 引起 的 
偏差 ， 可 以 得 到 扩散 原子 云 的 二 维 速度 分 布 . 它们 包含 着 原子 
的 许多 热力 学 信息 ， 如 : 分 布 曲线 下 所 包围 的 面积 ， 正 比 于 总 
的 原子 数目 ; 在 零 速度 附近 出 现 的 罕 特 征 峰 ， 其 峰值 曲线 下 所 
包围 的 面积 的 积分 正比 于 处 于 体系 基态 的 原子 数目 ; 从 扩散 原 
子 云 的 平均 半径 和 扩散 时 间 可 以 得 到 原子 的 平均 扩散 速度 和 特 
征 能 量 等 特征 参量 . 冷 原 子 温度 常常 通过 飞行 时 间 的 测量 来 完 
成 ， 用 共振 荧光 方法 测量 距 冷 原子 团 中 心 一 定 距离 处 的 荧光 强 
度 的 时 间 演 化 ， 可 以 得 到 冷 原子 的 速度 分 布 ， 由 此 推出 温度 Po 
由 于 利用 共振 光 探 测 ， 原 子 会 强烈 地 散射 共振 光子 ， 从 而 引起 
对 原子 的 加 热效应 ， 因 此 这 种 共振 吸收 成 象 技术 对 形成 的 BEC 
有 一 定 破坏 性 . 


1.3 碱 金属 原子 的 BEC 实验 


1.3.1 HO fff] BEC 实验 


JILA 小 组 首次 实现 碱 金属 原子 BEC 实验 的 核心 装置 是 
磁 光 阱 (MOT: magneto-optical trap) 方 柱 形 玻璃 容器 室 ， 两 个 
平行 反 向 电流 线圈 产生 非 均 匀 静 磁场 ， 零 场 强 点 正好 在 气 室 中 
心 ; 三 对 两 两 相对 的 激光 束 沿 三 个 坐标 轴 交 会 在 气 室 中 心 ， 激 
光 的 频率 略 低 于 铭 原 子 的 所 用 跃迁 频率 ( “ 红 ， 失 谐 ) . 激光 
线 宽 小 于 1MHz， 气 室 真空 度 达 ~133x10-? Pa. 它们 利用 磁 光 
BEA BEIC He P B zi A n pil SE RA HUE T, 然后 通过 调节 磁场 
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梯度 和 激光 频率 ， 把 原子 云 进 - 步 庄 缩 并 冷却 ， 接 着 关 掉 所 有 
激光 , 在 Ims 时 间 内 加 上 形成 TOP ( time orbiting potential ) Bf 
所 需 的 四 极 磁 场 和 旋转 磁场 ， 继 续 用 TOP 阱 囚禁 原子 ， 再 用 
射频 场 对 原子 样品 进行 燕 发 冷却 . 在 70s 时 间 内 ， 射 频 场 的 频 
率 和 旋转 场 强度 都 减 小 ， 射 频 场 的 频率 变化 时 ， 水 平 截面 内 原 
子 速度 分 布 出 现 明显 变化 ，BEC 开始 发 生 在 射频 场 频率 处 于 
4.25-4.23MHz 之 间 的 时 候 . 当 射 频 场 频率 大 于 4.23MHz Hj, 
原子 速度 分 布 为 平滑 的 单一 高 斯 分 布 ; 当 射 频 场 频 率 等 于 
4.23MHz 时 ， 原 子 速度 分 布 中 心 开 始 出 现 一 个 尖锐 的 峰 : 当 
射频 场 频率 小 于 4.23MHz 时 ， 原 子 速度 分 布 曲线 出 现 两 部 分 : 
一 部 分 为 宽 的 平滑 的 高 斯 分 布 ， 对 应 于 正常 相 原子 ， 而 另 一 部 
分 为 位 于 中 心 的 罕 峰 ， 对 应 于 凝聚 相 原子 . 当 射 频 场 频率 为 
425MHz 时 ,原子 云 中 有 2x10* 个 原子 ， 密 度 为 2.6x102 cm, 
温度 为 170 nK. 


T 


T 


1.3.2 钠 原 了 的 BEC 实验 


MIT 小 组 在 “光学 塞 孔 中 实现 Na 原子 的 BEC. 过程 如 下 3: 
经 激光 冷却 后 的 Na 原子 装 人 到 磁场 梯度 为 130G/em 的 “ 光 
学 塞 孔 ” 磁 阱 中 ， 峰 值 密度 约 有 10'"cm-:， 温 度 ~200 uK. 被 
内 禁 的 原子 通过 射频 场 进 一 步 蒸发 冷却 . 在 7s 内 ， 射 频 场 频 
率 从 30MHz 变化 到 1MHz ， 同 时 磁场 梯度 增加 至 550G/cm ， 
然后 降低 到 180G/cm. 凝聚 态 的 观测 是 先 用 共振 寺 
5=1 一 上 =2 的 激光 将 原子 抽 运 到 严 =2 态 ，10hs 之 后 用 调 
谐 于 F=2 一 请 =3 的 激光 进行 吸收 成 象 ， 探 测 光 成 象 CCD 
装置 ， 当 温度 高 于 15AK 时 ， 囚 禁 的 原子 云 呈 椭圆 形 . 当 温 
变 低 于 ISUK 时 ， 原 子 去 在 势能 最 低 点 处 分 离 成 两 个 目 模 ; 
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当 射 频频 率 wy 小 于 0.7MHz 时 ,可 观察 到 速度 分 布 对 称 性 发 生 
明显 变化 ， 大 于 此 频率 时 ， 分 布 是 球形 的 ， 如 同 热 运动 的 原子 
， 低 于 此 值 时 ， 速 度 分布 包 括 -个 椭圆 形 核 . 当 射 频频 率 进 
arenes 其 亮度 增加 ， 球 形 部 分 变 暗 . 这 个 椭圆 形 原子 云 
即 为 发 生 凝 聚 的 部 分 ， 球 形 原子 云 对 应 非 凝 聚 部 分 ， 椭 圆 图 形 
是 由 于 磁 阱 中 囚禁 势 的 不 对 称 引起 的 . 低 于 此 跃迁 频率 时 ， 出 
现 双重 速度 分 布 ， 罕 分 布 对 应 凝聚 部 分 ， 宽 分 布 对 应 非 凝聚 部 
分 ， 凝 聚 部 分 的 速度 分 布 不 再 是 高 斯 型 ， 两 重 分 布 的 宽度 不 一 
FÉ. 由 TOF 图 象 ， 得 到 原子 数目 N=5x10， UE fü m E 
ne =4x10 cm3， 凝 聚 态 寿命 大 约 1s， 温 度 T=2 pK. 
fg "Cloverleaf " BEC 实验 过 程 如 下 20: 利 
用 塞 曼 减速 和 上 暗 SPOT Db. E 2s 内 ， 有 (5-9)x10? 个 原子 
被 因 禁 ， 密 度 是 Sx10ucm-3:， 在 Is 内 快速 关闭 冷却 光 ， 打 
JF “Cloverleaf” — ( 阱 轴 向 曲率 最 大 值 是 1253Gcm-， 径 向 
梯度 170Gcm-!， 轴 向 偏 压 场 100G ). 原子 进入 阱 后 ， 得 到 
处 于 所 =1,m, = 一 1 态 的 原子 数目 大 约 2x10” 个 ， 峰 值 密度 
6x10" cm， 然 后 减 小 偏 压 场 到 1G， 接 着 用 射频 诱导 进行 
蒸发 冷却 . 射频 频率 从 30MHz 变 到 IMHz 的 时 间 是 15-26s . 
在 射频 场 关闭 后 ， 原 子 在 阱 中 热平衡 时 间 大 约 是 100ms， 
最 后 关闭 磁 阱 ， 对 原子 云 进 行 成 象 ， 从 原子 云 的 TOF 图 象 
可 以 看 到 ， 原 子 云 从 最 初 “铅笔 状 " 开始 扩散 ， 非 凝聚 衣 
分 均匀 扩散 ， 凝 聚 部 分 主要 沿 X 轴 扩 散 . 飞行 时 间 和 CCD 
成 象 反 映 了 原 于 云 的 速度 分 布 . fü \ 了 于 或 等 于 相 变 温度 
时 ， 正 常 部 分 速度 分 布 是 均匀 的 ;在 温度 T=2 uK 时 ， 
到 凝聚 态 峰 值 密度 nm =1.5x10 cm", BRAK FR 
N.=5x10 ， 是 光学 寨 孔 阱 中 数 日 的 10 4%. 从 相互 作用 能 与 
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的 扩散 独立 二 正常 成 分 ， 表 明 止 常 成 分 的 扩散 不 影响 凝聚 
态 的 扩散 . 在 射频 场 工 作 状态 下 ， 凝 聚 态 的 寿命 为 20s; 去 
掉 射 频 场 时 ， 压 命 减 小 到 ds. 这 是 由 于 背景 气体 掠 人 射 碰 
SET AAS. MER SUR THOA OH AE E E A ERE , 
LES BE Sb. 
最 近 的 钠 原 子 在 不 对 称 简 谐 势 阱 ( o, = c, = 2050 rad/s, 
€. =170 rad/s ) 中 的 BEC 实验 局 ， 其 柱 密度 分 布 如 图 1-1 所 示 . 给 


Optical Density 
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Zum) — 41.7 56.1 66.9 


图 1-1 钠 原 节 在 不 对 称 简 谐 势 阱 中 的 BEC FERRE o 
(e, =w, =2050 rad/s. «o, =170 rad/s) 


实 线 : 凝聚 原子 数 分 别 为 No =8x107. 35x105, 85x10 
均 场 的 计算 结果 . 相应 于 峰值 处 的 柱 密度 分 别 为 422、1072， 和 1826 
atoms/ p mò. EER N, AE x 10% 后 的 计算 结果 .箭头 给 出 三 方向 下 
Thomas-Fermi 表面 位 置 . 实验 数据 实心 圆 图 于 的 竖 线 到 示 统计 误 准 
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出 FEE TORT BEC 的 空间 分 布 的 定量 测量 ,图 中 曲线 表明 : 
灾 验 结果 与 平均 场 理论 很 好 地 -… 致 ;: 可 在 1% 以 内 确定 凝聚 体 的 
Thomas-Fermi 表面 . 
1.3.3 Hee fff] BEC 实验 
Rice 大 学 的 R.G.Hulet 等 人 报道 了 锂 原子 进入 同一 量子 
简 并 区 的 实验 现象 Pi， 其 实验 过 程 如 下 : 锂 原子 首先 利用 塞 
曼 减速 器 减速 ， 然 后 进入 loffe 阱 中 心 区 域 . 三 对 调谐 于 2S 
F=2, m, = 2 到 2Pj4, F-3, m, = 3 跃迁 频率 的 激光 相向 沿 
轴 向 人 射 到 阱 中 心 处 . 通过 多 普 勒 冷却 ， 锂 原子 温度 降 到 
200 LK， 峰 值 密度 ~7x10" cm-， 接 着 利用 微波 场 进一步 藻 
发 冷却 . 在 蒸发 冷却 5 分 钟 后 ， 气 体 发 生 相 变 ， 产 生 BEC, 
温度 T=300nK， 原 子 数目 N=1x10;. 一 个 低 于 原子 共振 频 
率 ， 失 谐 为 12~18MHz 的 高 斯 光束 被 作为 探测 光 来 测量 原子 
的 数目 和 温度 . 由 于 阱 的 弹性 常数 是 非 对 称 的 ， 图 象 呈 椭圆 
形状 ， 他 们 在 实验 上 发 现 ， 当 >15 时 ， 围 绕 吸收 峰值 处 


出 现 一 个 “ 光 晕 ”， 这 是 由 于 探测 光 与 形成 凝聚 的 原子 发 生 了 
衍射 效应 ;而 而 当 六- < 1 时 ， 衍射 效应 消失 . 后 来 ， 他 们 通过 


相 衬 成 像 技 术 测量 了 形成 凝聚 的 锂 原子 数目 3， 如 图 1-2 所 
AN. 观察 到 量子 简 并 区 的 条 件 是 在 原子 数 N 为 6800~135000 
和 温度 了 为 120-350 nK, 实验 上 确定 的 简 并 区 的 半径 及 为 
2.5-5 pm. 在 所 有 情况 下 凝聚 数 No 都 是 很 小 的 ，R = 5 pm 
时 求 得 的 N。 约 为 1300， 这 个 值 当 只 =4 hm 时 减 小 到 1000 和 
R=2.5 um 时 的 650. 
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Signal 


Radius (pm) 


图 1-2 WEERA EE TARIA 

竖 直 轴 表 示 信 号 强 度 相对 于 探 针 强度 的 大 小 ， 数 据点 由 所 测量 的 图 象 信 
号 得 出 ， 实 线 是 按 玻 色 - 爱 因 斯 坦 统计 分 布 所 计算 的 结果 , 短 虚 线 是 将 凝 
聚 原子 去 掉 问 样 分 布 的 计算 结果 . 长 虚线 是 经 典 ( 高 斯 ) 分 布 的 计算 结果 . 
考虑 到 成 像 系统 的 分 辩 率 (假定 是 4 po m 的 有 效 分 辨 率 ) 的 限制 . 在 上 述 曲 
线 的 计算 中 包括 与 贡 斯 函数 卷 积 . 在 图 (a) Duft i +191MHz 温度 
在 309 nK. 9.0x 10 个 原子 的 分 布 . 凝聚 原子 数 在 ~1 的 量 级 .表明 气体 
刚 接 近 简 并 . 图 (b) : 探 针 失 谐 -130 MHz. 按 信号 由 强 到 弱 依 次 对 应 的 
总 原子 数 和 温度 分 别 为 1.01x105 和 304 nK、2.6x104 和 193 nK. 
6.6x10 和 122 nkK， 它 们 相 忆 的 凝聚 康 了 数 分 别 为 500、810 和 270. 
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1.3.4. FL BHL USC Pf BEC 实验 


随 着 激光 冷却 和 约束 原 了 的 技术 向 前 推进 ， 再 加 上 蒸发 冷却 
与 做 阱 约束 ， 终 于 实现 了 中 性 原子 氧 的 玻 色 - 爱 因 斯 坦 凝 聚 . 
MIT 的 Kleppner 小 组 最 早 研究 所 原子 的 BEC， 在 1991 年 就 已 获得 
100 u KAY AEGAN 8 x 10° cm 原子 密度 PH， 该 密度 所 对 应 的 临界 
温度 为 30hK， 应 该 说 离 实现 BEC 条 件 不 还 了 ， 可 是 由 于 一 对 自 
旋 反 平行 的 氧 原子 在 第 三 个 氨 原 子 的 碰撞 下 会 复合 成 毛 分 子 
HT «Hi >H,+4.6 eV. 复合 速率 正比 于 密度 的 三 次 方 ， 从 而 限 
制 了 密度 的 增加 ; 发 热 则 妨碍 进一步 降温 ， 这 些 均 在 很 大 程度 上 
推迟 了 和 氧 原子 BEC 的 实现 . 在 经 过 儿 年 的 努力 后 ，MIT 的 Kleppner 
小 组 在 1998 年 终于 实现 了 氧 原子 的 BEC i 技术 仍 是 激光 冷却 与 
约束 ， 并 结合 蒸发 冷却 、 磁 阱 等 技术 .达到 的 握 原 子 密度 为 
10? em?, iE 7950 AK， 处 于 凝聚 态 的 原子 数 为 108 ， 约 10 售 
于 钠 原子 凝聚 体 的 原子 数 . 氧 原子 以 其 很 弱 的 相互 作用 ， 是 最 早 
选 作 进行 BEC 实 验 的 原子 ， 但 真正 实现 BEC 要 比 碱 金属 原子 晚 了 
JLE. 在 最 终 实现 氢 原 子 的 BEC 实 验 中 ， 所 用 的 探测 方法 也 与 以 
前 的 办 法 “将 原子 从 凝聚 体 中 引出 来 "， 并 测 其 速度 分 布 中 不 一 
PE, MERIO FIKRI. 即 氨 原子 同时 吸收 两 个 光子 ， 由 1S 
人 态 跃 迁 到 2S 态 ， 实 验 中 加 电场 时 2S 态 与 靠 得 很 近 的 寿命 很 短 的 
2D 态 混合 ， 原 子 由 2P 态 衰变 并 辐射 出 Lyman a 线 ， 通 过 对 
姑 光 了 于 的 计数 检测 出 经 由 双 光 子 吸 收 跃 迁 到 2S 态 氢 原 子 数 .图 
1-3 为 计数 率 对 激光 失 谐 的 双 光 子 谱 . 按 实验 的 安排 所 原子 吸收 
双 光 子 可 分 为 两 种 情形 . (1) 吸收 两 个 反 向 运动 的 光子 , 原子 获 
得 的 净 动 量 为 零 ， 没 有 光子 的 反弹 移 位 或 Doppler 如 宽 ， 如 图 中 
那个 塞 近 原点 “0” 的 尖峰 . (2) 吸收 同 向 运动 的 两 个 光子 ， 双 光 
了 的 动量 及 反弹 能 均 被 原子 所 接收 ， 使 原子 产生 约 6.7MHz 
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B13 BR SURO DOR gH! 
图 中 空心 圆 图 代表 无 凝聚 样品 的 谱 ; 实心 圆圈 为 有 凝聚 特征 的 谱 . 
原点 附近 的 尖峰 来 源 于 吸收 两 个 反 向 运动 的 光子 ， 没 有 光子 反弹 
移 位 和 Doppler 加 宽 . 右边 的 宽 低 峰 来 源 于 鹃 收回 向 运动 的 两 个 光 
Ta 实 线 是 正常 气体 在 T = 40 u K 时 有 反弹 移 位 、Doppler 加 宽 
的 高 斯 线 型 . 右 侧 较 窑 的 Doppler 加 宽 详 (实心 圆圈 ) 代表 凝聚 体 
的 窒 的 动量 分 布 ， 而 左 侧 无 Doppler 移 位 的 峰 开 示 高 密度 原子 引 
起 的 能 级 移 位 . 图 中 专 失 谐 坐 标 是 相对 无 微 扰 原子 的 无 Doppler 激 
发 、 所 有 频率 者 是 相对 243nm 的 激发 辐射 
的 移 位 并 有 Doppler 加 宽 轮 廊 . 在 图 中 原点 “0” 的 右边 ， 这 加 宽 
体现 了 未 凝聚 原子 气体 的 温度 . 图 中 的 原点 “0” 是 双 光 子 激光 
MU 2c, 等 于 1S-2S 跃 迁 频率 处 (对 应 的 跃迁 波长 为 243nm)， 对 于 
第 (D) 种 情形 ， 原 子 吸 收 了 两 个 反 向 运动 的 光 了， 没有 反弹 能 引 
起 的 移 位 及 Doppler 加 宽 , 计数 峰值 应 严格 与 原点 “0” 重合 ， 但 
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实际 上 还 是 向 左 移 位 了 18kHz， 这 是 由 邻近 的 高 密度 原子 ( 达 
108cm253) 引 起 的 所 原子 能 级 移 位 ， 移 位 正比 于 原子 密度 . 我 们 
还 注意 到 在 零点 左边 约 400 ~ 500 kHz 处 的 计数 峰值 ， 也 是 高 密 
度 原子 引起 的 能 级 移 位 ， 对 应 的 原子 密度 约 为 气体 原子 密度 的 20 
多 信 ， 而 这 就 是 高 密度 原子 凝聚 体 的 明证 . 我 们 还 注意 到 在 右 
侧 的 Poppler 加 宽 谱 峰值 附近 也 发 现 高 密度 原子 凝聚 体 . 当然 重 
要 的 是 原点 左边 的 约 400~500 kHz 处 的 峰 ， 真 正 代 表 了 无 
Doppler 移 位 的 双 光 子 吸收 谱 ， 这 是 一 个 线 宽 接近 于 自然 线 宽 的 
超 冷 原子 源 . 


1.4 BEC 研究 的 理论 模型 加 


研究 玻 色 凝聚 原子 体系 行为 的 自然 出 发 点 为 Gross- 
Pitaevskii 理论 形式 的 弱 相 互 作用 玻 色 子 理论 ; 这 种 平均 场 方法 
提供 了 描述 BEC 及 有 关 现 象 更 接近 和 较 简单 的 方程 ， 这 个 理论 
适 于 描写 稀薄 气体 中 二 体 相 互 作用 的 大 多 数 效应 . 文献 [26] 用 数 
值 计算 方法 讨论 了 1 维和 3 维 谐振 子 势 下 基态 波 函数 的 时 间 演 
化 ， 特 别 是 当 势 阱 随时 间 变 化 或 受 扰动 的 情况 . 该 研究 还 模拟 
了 形成 BEC HURT TE A A Sil BARAT BC A FS Bl ADL DIT IK 
函数 的 动力 学 演化 

有 许多 文章 研究 了 吸引 型 原子 体系 的 BEC 问题 ， 一 般 认 为 
具有 负 散 射 长 度 的 锂 原子 均匀 体系 不 能 形成 稳定 的 BEC 28271, 
这 是 因为 它 的 s 波 散 射 长 度 亏 是 负 的， 原子 间 的 吸引 相互 作用 
将 引起 凝聚 体 的 崩塌 . 然而 ， 文 献 [26,28,29] 预期 在 充分 少 的 原 
子 数 或 非常 弱 的 相互 作用 条 件 下 凝聚 体 应 是 稳定 的 ， 当 这 些 条 
件 不 满足 时 的 衰变 和 崩塌 也 有 几 位 作者 讨论 P?30. 

弱 相 互 作用 的 稀 沙 气体 的 BEC 的 理论 描述 已 有 很 长 的 历 
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史 ， 并 能 解释 大 多 数 实验 结果 !31. 相 比 之 下 ，BEC 形成 的 动力 
学 的 完全 描述 还 没有 建立 起 来 . 对 于 凝聚 形成 所 需要 时 间 的 长 
短 ， 早 期 的 预计 是 无 限 长 的 54， 这 是 根据 玻 尔 兹 曼 方 程 ， 凝 聚 
不 可 能 从 零 生 长 . 这 样 ， 一 个 分 离 的 成 核 过 程 一 定 在 成 核 之 前 52 
和 成 核 之 后 中 引入 到 描述 动力 学 的 玻 尔 兹 曼 方程 . Stoof 05591 提 
出 ， 在 短 的 相干 阶段 凝聚 成 核 并 根据 动力 学 方程 生长 . Kagan 等 
人 5 讨论 了 准 凝 聚 体 的 形成 和 相位 涨 落 . 最 近 ， 玻 色 气 体 的 全 
量子 力学 的 动力 学 理论 已 提出 并 用 于 模型 化 凝聚 的 形成 过 程 Bal， 
BEC 实验 上 的 实现 已 证 实 凝聚 体 是 在 有 限 的 时 间 内 形成 的 ， 同 
样 ， 两 凝聚 体 之 间 的 高 对 比 度 干涉 的 观察 证 实 凝 聚 体 在 有 限 的 
时 间 内 演化 为 长 程 相干 1. 凝聚 生长 速率 与 平衡 时 凝聚 原子 数 的 
关系 的 实验 结果 5 提示 出 : 发 生 凝 聚 的 两 个 阶段 ， 快 速 弛 瑰 产 
生 一 个 过 饱和 热 云 ， 随 后 在 热 云 中 进行 凝聚 的 慢 速生 长 . 目前 
的 理论 结果 仅 与 凝聚 第 二 阶段 的 实验 结果 一 致 . 

元 激发 的 研究 是 量子 多 体 理论 的 重要 任务 之 --， 对 于 玻 色 
流体 ， 元 激发 在 理解 其 超 流 特性 时 起 关键 作用 . BEC 实现 以 后 ， 
所 激发 的 最 低 模式 的 频率 的 测量 和 波 包 传 播 的 观察 都 已 得 到 . 
理论 工作 的 出 发 点 是 由 含 时 的 Gross-Pitaevskii 方程 导出 的 耦合 
方程 a 可 用 于 计算 激发 的 本 征 频 率 ， 它 的 引入 是 用 于 研究 均匀 
玻 色 气体 中 的 涡 旋 线 激发 . 这 种 方法 与 在 Bogoliubov 近似 下 哈 
密 顿 的 对 角 化 是 等 价 的 ， 文 献 [41] 在 均匀 玻 色 气体 中 准 粒 子 振幅 
为 平面 波 的 条 件 下 求 得 了 声 子 色散 关系 . 在 简 谐 热 阱 条 件 下 ， 
文献 [42] 计 算 了 捕 陷 气体 的 激发 . 在 捕 陷 气体 所 展示 的 各 种 激发 
中 ， 偶 极 模式 特别 引 人 注 意 ， 这 种 振动 模式 相应 于 系统 质心 的 
运动 . 文献 [43] 在 Bogoliubov 理论 框架 下 详细 讨论 了 这 种 偶 极 激 
发 的 特性 . 
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要 进行 BEC 光学 性 质 的 研究 ， 必 须 建 立 BEC 与 光 相 互 作 
用 的 理论 . 原子 BEC 是 特殊 状态 下 的 原子 系统 ， 其 与 光 相 互 作 
用 必然 有 不 同 于 普通 原子 系统 的 特点 . 文献 [44] 研 究 了 无 限 空 间 
中 均匀 的 BEC 对 光 的 散射 , 结论 是 , BEC 强烈 地 反射 近 共振 光 ， 
因此 在 激发 谱 中 共振 频率 处 会 出 现 一 个 禁 带 . 对 于 外 势 场 中 的 
BEC， 并 考虑 原子 间 的 相互 作用 ， 文 献 [45,46] 讨 论 了 强 的 脉冲 
激光 的 散射 ， 并 提出 可 利用 这 种 散射 光 作 为 探测 BEC 的 一 种 方 
法 . 而 文献 [47,48] 讨 论 了 弱 的 连续 光 的 散射 ， 指 出 其 散射 谱 是 一 
个 宽 的 非 罗 仑 兹 线 型 的 背景 ， 加 上 共振 频率 处 的 一 个 罕 峰 这 
一 系列 工作 都 显示 出 高 于 临界 温度 和 低 于 临界 温度 使 光 的 散射 
有 很 大 的 不 同 . 这 也 逐步 建立 了 较 完 整 的 BEC 与 光 相互 作 用 的 
理论 (448391. 

紧 束 缚 的 BEC 对 光 的 散射 的 研究 So， 表 明 这 时 BEC 比 通 
常 相 互 独立 的 原子 散射 要 少 ; 而 当 BEC 足够 大 时 ， 散 射 增加 ， 
并 主要 沿 前 向 散射 . 文献 [51.52] 研 究 了 非 共振 光 的 散射 ， 他 们 的 
计算 显示 了 BEC 的 散射 谱 具 有 双 峰 的 特性 ， 对 BEC 的 其 他 光 
学 性 质 ， 如 自发 辐射 性 质 呈 5 约 ， 荧 光 性 质 5S9， 非 线性 光学 响应 6569 
近 几 年 也 都 得 到 了 一 些 理论 结果 . 关于 EBC 的 原子 光学 和 非 线 
性 原子 光学 方面 的 研究 也 取得 了 令 人 触目 的 进展 . 文献 [57] 报 道 
了 根据 双 光 子 Raman 散射 技术 (一 级 Bragg 衍射 ) 产 生 的 高 度 准 
直 的 准 连续 原子 激光 . Lu Deng 5 报道 了 物质 波 的 多 波 混 频 实验 ， 
为 非 线 性 原子 光学 的 研究 打下 了 基础 . 
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BEC 的 实现 ， 提 供 了 在 体系 基态 上 宏观 布 居 的 样品 原子 ， 
其 物质 波 可 由 单 -的 宏观 波 函 数 来 描述 . 利用 特定 的 手段 从 了 
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中 相干 输出 的 原子 的 运动 行为 将 服从 单 模 波动 方程 的 描述 . 近 
期 最 引 人 注 目的 成 就 是 “原子 激光 "的 问世 5， 即 类 似 于 激光 
的 相 十 原子 束 输 出 . 但 是 ， 由 十 原子 具有 静 质 量 和 原子 间 的 磁 
接 ， 会 使 原子 激光 的 相 十 性 复杂 化 I， 原子 激光 的 获得 可 望 有 
效 地 提高 现 有 原子 钟 的 准确 度 : 相干 原子 束 的 平均 速度 很 低 ， 
有 利于 减 小 与 多 普 勒 效应 有 关 的 不 准确 度 . 原子 激光 的 男 一 个 
可 能 的 应 用 领域 是 纳米 技术 . 目前 ， 纳 米 技术 是 利用 扫描 隧道 
显微镜 和 原子 力 显 微 镜 来 搬移 原子 和 光 刻 蚀 的 方法 来 产生 微小 
的 特征 量 . 由 于 原子 激光 可 以 聚焦 并 像 普 通 激光 一 样 准 直 ， 因 
此 可 望 利 用 原子 激光 通过 适当 的 变换 直接 形成 所 需要 的 花样 
此 外 ， 原 子 激光 具有 明显 的 干涉 效应 ， 原 子 本 身 又 具有 非 零 的 
静 质 量 、 有 自 旋 和 内 能 等 特点 ， 因 此 原子 激光 还 可 能 极 大 地 提 
高 重力 加 速度 和 其 他 一 些 物理 常数 的 测量 精度 . 还 会 有 难以 预 
料 的 用 途 将 给 科学 技术 带 来 丰富 多 彩 的 应 用 前 景 . 


1.6 本 文 的 安排 与 我 们 的 工作 


第 一 章 简要 回顾 了 BEC 研究 的 历史 , 介绍 了 实现 BEC 的 相 
关 技术 和 碱 金 属 原子 和 氢 原 子 BEC. 的 主要 实验 结果 和 近 几 年 
BEC 的 理论 研究 工作 . 第 一 章 给 出 无 相互 作用 玻 色 气体 的 基本 
性 质 ， 引 入 振子 长 度 ， 能 量 单位 和 临界 温度 等 . 讨论 外 加 势 场 
对 BEC 的 影响 ， 以 及 体系 有 限 原子 数 的 影响 和 维 数 的 作用 . 第 
三 章 研 究 排斥 相互 作用 对 BEC 的 影响 ， 首 先 根 据 该 Gross- 
Pitaevskii 方程 得 到 定 态 非 线性 S 方程 ， 给 出 解 有 排斥 相互 作用 
情形 中 性 原子 玻 色 - 爱 因 斯 坦 凝聚 的 一 般 方法 . 并 应 用 于 有 排 
斥 相 互 作用 的 中 性 原子 体系 的 解析 求解 及 玻 色 - 爱 内 斯 坦 凝 聚 
计算 .第 四 毫 给 出 非 线性 定 态 薛 定格 方程 (NLSE) 的 : -种 数值 求 
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解 方 法 ， 并 求解 简 谐 势 阱 中 有 排斥 作用 的 中 性 原子 体系 NLSE 
的 基态 和 激发 态 解 . VETER TE E 2S BERE 1573 FE UD Ui ER BO SC 
与 归 一 化 问题 ， 并 对 计算 精度 进行 分 析 . 第 五 章 研究 吸引 相互 
作用 对 BEC 的 影响 ， 首 先 计算 有 吸引 相互 作用 的 中 性 原子 体系 
的 能 级 和 态 函 数 在 一 维 箱 势 中 的 解析 解 和 三 维 球 方 势 中 的 数值 
解 ， 分 析 原 子 间 的 吸引 作用 对 能 级 和 波 函 数 的 影响 ， 并 与 排斥 
相互 作用 情形 的 解 进行 了 比较 , 讨论 吸引 相互 作用 对 玻 色 - 爱 
因 斯 坦 凝聚 的 影响 . 根据 简 谐 势 阱 中 具有 吸引 相互 作用 原子 体 
系 的 非 线 性 定 态 薛 定 刘 方 程 NLSE) 的 基态 解 , 即 凝聚 原 子 数 随 
能 量 本 征 值 变化 的 双 稳 态 曲 线 ， 研 究 吸引 型 BEC 的 亚 稳 性 质 和 
宏观 量子 隧道 . 第 六 章 给 出 BEC 的 动力 学 模型 ， 用 母 函 数 法 求 
得 玻 色 凝聚 体 主 方程 的 解 ， 并 计算 凝聚 体形 成 的 增 率 与 粒子 数 
起 伏 的 统计 规律 . 第 七 章 给 出 一 准确 求解 入 原子 辐射 几率 的 模 
型 ， 计 算 玻 色 - 爱 因 斯 坦 凝 聚 体 的 自发 与 受 激 跃迁 几率 的 增强 因 
T. 第 八 章 给 出 结论 和 该 领域 进 -- 步 研究 的 前 景 . 
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第 二 章 ”理想 气体 的 玻 色 - 爱 因 斯 坦 凝 聚 


2.1 均匀 玻 色 气 体 的 BEC 
WAN 个 质量 为 m 的 无 相互 作用 的 玻 色 子 处 于 体积 为 
VAZI ( 三 维 箱 中 的 理想 玻 色 气体 系统 )， 其 哈密 顿 量 为 
H= $a (2-1) 


EP p?=p,-p,, p 为 第 i 个 粒子 的 动量 算 符 . 哈密 顿 量 
与 粒子 的 位 置 或 坐标 无 关 . 系统 的 能 量 本 征 值 是 如 下 单 粒 
子 能 量 (能 级 ) 之 和 

Ep =p? /2m (2-2) 
P= Pl, pe PAF A BATE: p=(2nh/Lyn, nf 
KAA 0 ALE EAM, L-y!^, 
该 系统 满足 玻 色 - 爱 因 斯 坦 分 布 2 


<1, >=>—_ 
p ue (2-3) 


其 中 <n, > 表示 处 于 动量 为 p、 能 量 为 en 的 能 态 的 平均 粒 
子 数 ，z 是 所 谓 “ 易 锡 度 "(fugacity)，==e%pay | <n, > 是 
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化 学 势 . 总 粒子 数 
N-«ny»*Y«n,» (2-4) 


p=0 

4V > okt, pA BRE MEE. XP pR E h oP 

RH, HY — v) fa? p, AERE, o 1 时 ， 
^ 


(2-4) 式 中 的 和 式 是 发 散 的 ， 因 为 对 应 于 s=0 的 那 一 项 趋 于 
无 穷 . 因此 ，e=0 时 的 这 一 项 可 以 像 整 个 和 一 样 重要 . 考虑 


到 这 一 点 ， 把 这 一 项 单 分 出 来 ， 把 其 余 的 求 和 以 积分 代替 ， 

并 利用 自由 粒子 态 密度 Dede - Qx) "VOm/ i? y" Jede, 4i 
N 1 z2 2r ype £? de 
y'wpiu o hoea € 


积分 下 限 仍 取 为 零 ， 这 是 因为 当 a =0 时 D(z) = 0， 这 时 对 
积分 无 贡献 . 第 一 项 中 的 2/1- :=) 则 代表 单 粒子 态 e=0 上 
的 平均 占有 数 Ng. 将 (2-5) 的 积分 作 变量 代 换 ， 可 进一步 
t 


N 1d z 1 
y^yi-; perO (2-6) 
AE MEAE xs d? agaa di 

其 中 g3,.(2) A = ff gm. aF f Slee , H 


81:200) < 83/20) = 2.612--, (OS 251), T(n) AMBRE; 
= 20h? /mksT 称 为 热 波 长 ， 与 质量 为 加 、 能 量 为 ks7 粒 
了 的 德 布 罗 意 波长 同 量 级 . 
当 W/F> Sa/ 下 时 系统 中 的 粒子 以 宏观 上 有 限 的 
HAT OBER Yc = 0 单 粒子 态 上 ， 此 即 玻 色 - 爱 因 斯 坦 凝 聚 现 
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Ri 对 十 给 定 的 粒子 密度 n= N/V ， 由 等 式 


NIV = g3,.(0/43 (2-7) 
定义 一 个 临界 温度 : 
2 2/3 2/3 
2nh- N N 
E = — 2- 
Lip Gard A) (2-8) 
By = 下 2.612 3 ， 是 一 个 只 跟 原 子 质量 有 关 的 常数 ， 当 


km 

T <T: 时 ， 对 应 热 波 长 4 将 足够 长 ， 使 得 必 Tm DEM. 
g3/2(z)， 由 此 看 出 ， 形 成 BEC 的 条 件 是 

À » (2.612V / N)'^ (2-9) 


此 式 的 物理 意义 是 粒子 德 布 罗 意 波长 超过 粒子 间 的 距离 . 
用 临界 温度 7 的 语言 表达 : 凝聚 区 即 为 7<7 的 区 域 . 在 这 


Ari EXE, HV ot, z-1 ， 由 (2-6) 和 (2-8) 可 得 到 
凝聚 数 与 温度 的 关系 
No/N =1-(T/T, ? (2-10) 


同样 ， 对 于 给 定 的 温度 7 ， 由 等 式 (2-7) 也 可 定义 一 个 临界 
KE n = (N/V).: 


2 


N mT Y? (T ye 
ps -260| "57 | . 2- 
E ir] E (5) oe 


确定 形成 BEC 的 条 件 用 临界 参量 可 表示 为 


2 


MEME oem 
V Vj. \Bo 
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这 表明 : 玻 色 粒子 体系 存在 一 个 临界 密度 ， 当 粒子 数 密度 
>me 时 ， 和 再 增加 的 粒子 将 进入 基态 ， 这 些 增加 的 粒子 将 不 
对 系统 的 热力 学 性 质 有 贡献 ， 此 即 BEC REN, 

由 以 上 的 分 析 可 以 看 出 ， 可 以 从 两 个 途径 实现 BEC，… 种 
是 降低 原子 气 团 的 温度 ， 使 低 于 给 定 密度 下 的 临界 温度 ; 
种 是 提高 原子 气 团 的 密度 ， 使 超过 给 定 温度 下 的 临界 密度 . 


2.2 简 谐 势 阱 中 的 BEC 


目前 实验 中 BEC 的 形成 是 借助 某 种 形式 的 原子 势 ， 对 
于 碱 金 属 原子 BEC 实验 中 表征 其 磁 阱 的 重要 特点 就 是 束缚 
势 ， 它 可 近似 为 如 下 二 次 形式 


Vex =F lois? *oly! +22") (2-13) 


在 不 考虑 原子 间 相 互 作用 的 情况 下 ， 这 种 体系 的 几乎 所 有 
预期 结果 都 是 解析 的 . 多 体 哈 密 顿 就 是 如 下 形式 单 粒子 险 
密 顿 之 和 
Ennyn: = x + ho, +(ny tIio, + 
: A (2-14) 
(n. + Sho. 


式 中 onon) 是非 零 整数 . 1 (2-13) de NA 
作用 玻 色 子 的 基态 rory) 可 通过 把 所 有 粒子 都 放 在 最 
低 的 单 粒 子 态 而 得 到 (nm, = Om, = Om, = 0), PeZ- 


AN ; - 
PTEE) nes) ILectr) > xx B go) t Fx 


n 
is) 
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给 出 
-Porr toys? mz!) 
evo - (297 Y entm Q-15) 


这 里 引入 了 谐振 频率 的 几何 平均 

oy =(0,0,0- )3 (2-16) 
这 样 ， 密 度 分 布 函数 可 写 为 w(r) = Njo, 它 的 值 随 N 
值 而 增加 . 但 原子 云 的 大 小 与 N 无 关 . 定义 如 下 的 谐振 子 长 度 


ar (GU (2-17) 


mor 

它 相应 于 (2-15) 式 的 平均 高 斯 宽度 ， 通 常用 它 表 征 系统 的 空 
间 尺 度 . 在 有 限 温 度 ， 仅 原子 的 一 部 分 占据 最 低能 态 ， 其 
余 为 在 更 高 能 量 的 激发 态 上 的 热 分 布 . 热 云 的 半径 比 er Kk. 
假定 kgT »» hoy 和 具有 经 典 的 麦克 斯 韦 - 玻 耳 兹 曼 分 布 
na (r) exp[-Voa(m)/kagT] 的 热 原子 云 密 度 近似 条 件 下 可 以 大 
致 估计 其 大 小 . 如 果 V(r)=(1/2)mw?2r* ， 高 斯 宽度 是 
Ry =ay(kgT/ho,)'? , Wb, Ry tha, KS. 使 用 玻 色 分 
布 函数 对 这 个 估计 没有 明显 的 改变 . V.Bagnato 4% ANER 
地 讨论 了 普遍 寡 函 数 势 的 情况 ， 例 如 U(r) = str /1a)? 势 场 中 
的 临界 温度 、BEC 基态 布 居 比 例 为 


Ba (2-18) 


j BN 
T, =| 一 一 一 一 
12022? (2m)? 
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sn. (T«T) (2-19) 
上 上 面 的 讨论 揭示 出 简 谐 势 阱 中 的 BEC 在 密度 分 布 的 中 心 区 会 
有 尖锐 的 峰 出 现 . 图 2-1 给 出 球 对 称 简 谐 势 阱 中 5000 个 无 相 
互 作用 原子 在 温度 7= 0.97, 凝聚 和 热 原 子 密度 分 布 曲线 ， 它 
相应 十 柱 密度 ， 即 沿 一 个 方向 积分 之 后 的 密度 ， 
nE) = dxn(x,0, =) ， 这 是 典型 的 测量 量 ，x 方 向 为 用 于 原子 
云 成 象 光束 的 方向 . 若 直接 画 凝 聚 密度 和 非 凝聚 密度 之 比 
则 在 中 心 的 值 会 更 大 . 


400 
2 
E 300 Ho 
5 
uv 
z 
5 200 | 
5 
5 
100 4 
ny 
0 
-8 4 0 4 8 


Pel 2-1 球 对 称 简 谐 热 阱 中 5000 个 无 相生 作用 诛 子 在 温 度 
T = 0.9T。 BER AMIS TR MESS AB MER. PL OER UR DE MIESE 
WNT ib. SERS RIE HL a, Ma? ye. HEIN 
-化 为 原子 数 〈 取 自 文献 [631) . 


通过 基态 波 函 数 的 傅 里 叶 变 换 ， 可 以 计算 凝聚 体 中 的 
原子 的 动 晤 分布， 对 于 理想 气体 为 高 斯 塌 ， 中 心 在 零 动量 ， 
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EIER] Far. 热 云 的 分 布 在 动量 空间 中 更 宽 . 应 用 经 典 
分 布 函数 可 求 得 这 个 宽度 比例 于 Gs T) 7. 实际 上 简 谐 势 中 
的 理想 气体 的 凝聚 和 非 凝 聚 粒子 的 动量 分 布 有 与 图 2-1 密 
度 分 布 和 nm) 完全 相同 的 形式 . 

凝聚 体 在 坐标 空间 和 动量 空间 都 有 一 个 罕 蜂 出 现 这 一 
特征 是 很 独特 的 ， 这 对 阱 中 玻 色 气体 凝聚 的 实验 和 理论 分 
析 都 有 重要 价值 . 这 不 同 于 粒子 凝聚 进入 零 动 量 状态 的 均 
匀 气 体 情况 ，BEC 不 可 能 在 坐标 空间 中 显示 出 来 ， 因 为 凝 
聚 的 和 非 凝 聚 的 粒子 都 充满 了 相同 的 体积 . 


2.3 简 谐 势 阱 中 玻 色 气体 在 有 限 温度 的 BEC 
在 巨 正则 系 综 有 限 温 度 了 粒子 总 数 可 由 下 面 的 求 和 得 到 


N= Y (expl(e, uu, 7 40/ kgT] - ^ (2-20) 


而 总 能 量 为 

E= Y tox uu 72 kaT]-W (2-21) 
其 中 必 是 化 学 势 . 在 某 给 定 温度 以 下 最 低能 态 的 密度 变 成 
宏观 的 ， 这 就 相应 于 BEC 的 开始 . 临界 温度 、 基 态 布 居 比 
例 和 其 他 热力 学 量 的 计算 都 可 以 在 适当 的 能 谱 s,, 根据 
(2-18) 和 (2-19) 式 而 得 到 . 


^ 


EE (2-22) 
IHD = (w, +o, +w-)13 是 势 阱 频率 的 算术 平均 .将 它们 
(KAR Ast. alti 
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N-No= 
? P exp[A(@,n, +w,n, *o.n.)/ kgT]-l 


(2-23) 


为 了 计算 这 个 和 式 的 解析 表达 式 ， 通 常 假定 当 N > o Ht fE 
级 间隔 变 得 越 来 越 小 以 至 于 求 和 可 用 积分 代替 
7 dn, dn, dn. 

Eee arer +@,n, +@.n,)/kgT)-1 
这 个 假定 相应 于 激发 态 的 半 经 典 描写 . 计算 的 有 效 性 决定 
于 对 于 对 (2-23) 式 求 和 有 贡献 的 激发 态 能 级 是 否 远 大 于 由 振 
子 频率 所 决定 的 能 级 间隔 . 如 果 阱 中 的 原子 数 很 大 和 
kp7 >> ho7 半 经 典 近 似 被 期 望 是 一 个 好 的 近似 . 通过 比较 积 
分 式 (2-24) 和 求 和 式 (2-23) 的 值 来 检验 . 

积分 式 (2-24) 通 过 变量 代 换 可 以 容易 求 得 


(2-24) 


3 
N-No= “cl 和 (2-25) 


xh co) E SE & BK, or 是 几何 平均 (2-16) 式 ,从 这 个 结果 
可 得 到 BEC 的 转变 温度 . 事实 上 ， 在 转变 点 使 No 0, 我 
们 得 到 


N 
£a 
将 上 面 的 表达 式 (2-26) 代 入 (2-25) 式 ， 可 以 得 到 当 T<T BERE 
聚 比 例 对 7 了 的 依赖 关系 为 


3 
No _,_| T_ 
ye 5 (2-27) 


V3 
Kg. -ol ) = 0.94hw, N! (2-26) 
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同样 的 结果 也 可 以 重 写 (2-24) 式 通过 对 能 量 的 积分 得 到 


iz ay ee p(e)de 
d 0 exp(fe) - 1 


其 中 pe) 是 态 密度 .可 由 能 谱 (2-13) 式 求 得 
P(E) =(I/2)(hovo) E. 将 这 个 值 代入 (2-28) 式 可 再 求 得 结果 
(2-25) 式 . 由 积分 E= f plee de expe) -1] 可 求 得 系统 
总 能 量 (2-21) 式 的 结果 为 


4 
E  _3¢(4)( T 
PEN CL ey 这 过 -29 


(2-28) 


MoT, ÇO 


MAER RAER, HAN EL th ER. 以 上 这 些 
结果 可 以 和 第 一 节 均 匀 玻 色 气体 的 理论 结果 进行 比较 ， 
可 通过 半 经 典 近 似 计算 的 另 一 个 量 是 热 粒 子 密度 nj Cr). 


nr) 与 凝聚 密度 nm(r) = Nolpo(r)| 之 和 给 出 总 的 密度 
n(r) = n(n) + p(t). YET «T, 和 热力 学 极限 下 ， 热 粒子 密 
度 可 通过 对 动量 空间 的 积分 求 得 ， nj(r)= 
[d PTA fexp(Be(p 7) - V. stp epr) = p! /2m 4 Vou (1) 
是 相 空间 的 半 经 典 能 量 . 结果 是 

ny (e) = AP gi (e P 0) (2-30) 
Hp A, - uh. mkgTY? 是 热 波长 .函数 gis) T js 
BA gal) = 3 az" 1n* .通过 my(r) 的 空间 积分 可 再 -… 
求 得 与 (2-27) 式 一 致 的 结果 . es 
空间 的 热 粒子 分 布 : 
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ny (p) = (A, moy, )^ gs (oHm EA (2-31) 

以 上 的 分 析 指 出 理想 气体 存在 能 量 的 两 个 尺度 : 转变 
温度 AyT, ALF H3 fie t ju] BA hop. 从 (2-26) 式 可 清和 你 地 看 刘 
ksT IAF ho,. 在 合适 的 势 阱 中 ，W 值 的 范围 从 几 干 
到 几 百 万 ， 转 变温 度 是 joy 的 20 到 200 fi. 这 也 意味 着 半 
经 典 近 似 在 感 兴趣 的 温度 范围 内 可 得 到 很 好 的 期 望 结果 . 
由 午 势 给 定 频率 or /2r) 典 型 的 范围 是 10 ELAR, A 
IE hho, AILS nK BY at BR. 


2.4 有 限 原 子 数 对 BEC 的 影响 


能 进入 势 阱 中 的 原子 数 并 不 真正 是 宏观 的 ， 到 目前 为 止 
所 进行 的 实验 中 最 大 的 凝聚 原子 数 约 为 107. 因此 ， 热 力学 极 
限 从 未 严格 达到 . 第 一 个 影响 是 热力 学 函数 缺乏 突变 性 . 因 


而 ， 实 际 上 当 温 度 降低 时 可 观察 到 最 低能 态 的 宏观 占有 数 有 
很 突然 地 变化 ， 如 图 2-2 所 示 . 这 种 相 变 实际 上 是 相对 于 
Nw 期 待 值 的 近似 ， 但 这 种 影响 很 小 ， 对 于 有 限 粒 子 数 体系 
的 相 变 、 临 界 温度 等 术语 仍 是 相当 重要 的 . 

为 了 求 得 无 相互 作用 玻 色 气体 的 热力 学 性 质 ， 需 要 知 
道 玻 色 分 布 函数 中 的 单 粒子 能 谱 : 如 在 巨 正 则 系 综 中 由 求 
AIX (2-23) ERE WMP HM, RAL EN -> om 极限 ， 对 
于 给 定 的 原子 数 和 温度 ， 化 党 势 是 原子 数 和 温度 的 函数 求 
和 式 可 以 数值 地 进行 . 用 这 种 方法 得 到 的 凝聚 比 例 
NoT) N ， 比 热力 学 极限 预期 值 (2-27) 略 小 . 图 2-3 给 出 的 即 是 
10° 个 无 相互 作用 粒 FER ECL PRESE. 计算 表明 !k49 有 限 原 
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子 数 的 影响 仪 当 w 小 于 10 才 明 显 . 有 限 原 子 数 对 凝聚 函数 
(2-27) 式 的 修 止 ， 在 大 入 极限 已 有 解析 计算 1， 结果 是 


No a-[ T) ... 306(2) T | ii (2-32) 
N (S) 2e/ AGP? T 


这 种 影响 随 NCC ERE. ERC S EE A ERE o 5 
JUGE EI o, ZF RA PRE, i SB He T SE 
" 


_ 0- © A+2 


oO o 24^? 


0 0.2 0.4 0.6 0.8 1 12 1.4 1.6 1.8 
TIT? 
图 2-2 凝聚 粒子 数 比 例 No / N BETIT, De eO 

出 点 为 文献 [201 的 实验 结 来 ， 虚 线 为 (2-27) 式 的 理论 曲线 . 

24 N = 1000 时 ，(2-32) 式 的 期 待 值 已 与 通过 对 简 谐 振子 
哈密 顿 的 激发 态 求 和 所 得 的 严格 结果 难以 区 分 . 在 T. 附近 
窄 的 区 域 必 须 考虑 高 阶 修正 ， 才 能 得 到 严格 结果 ， 如 图 2-3 
所 示 ， 但 他 们 部 与 热力 学 极限 Nu ov =1- (六 人 下 略 有 不 同 
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0 0.2 0.4 0.6 0.8 1 


图 2-3 X V BL TRC No / N Bf 17 /的 变化 关系 
Fe ct A B tS F ER XF Eg (09 (8) 3 BE ON = 1000 4 CT 
格 的 量子 力学 计算 ， 实 线 相应 于 (2-32) 式 的 预期 ， 虚 线 为 
热力 学 极限 (2-27) 式 的 理 沦 值 曲线 


有 限 原 子 数 影响 会 降低 凝聚 的 比例 ， 因 此 导致 转变 温 


度 的 降低 (与 N 一 四 的 极限 比较 )， 使 (2-32) 式 左边 为 零 ， 可 
以 计算 临界 温度 的 相对 改变 量 !(G 49 


oT, TL(2 E 2, i 2 
oe -_ O80) SN VS 50.73 N75 (2-33) 
T.  2ow|£Q)] Oho 


2.5 空间 维 数 对 BEC 的 影响 


到 目前 为 止 ， 所 讨论 的 都 是 三 维 空 间 中 理想 玻 色 气体 


的 特性 . 尽管 在 势 阱 三 个 方向 上 的 频率 可 以 不 同 ， 凝 聚 体 
依赖 于 温度 的 有 关 结果 是 假定 了 ksT 远大 于 所 有 的 振子 能 


dt he.. ha,. he. 为 了 观察 降低 维 数 带 来 的 影响 ， 我 
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们 应 该 去 掉 - 一 个 或 两 个 方向 . 
一 维和 二 维 玻 色 气体 的 统计 行为 表现 出 很 独特 的 性 质 . 
在 有 限 温度 一 维和 二 维 的 均匀 气体 的 BEC RAAE RES, 
因为 热 涨 落 使 凝聚 体 极 不 稳定 . 由 于 BEC 状态 的 理想 气 
W, TESERORTE. xdg. n(p)e[exp(p? /2mkgT) - , 
显示 出 1/p? 红 外 发 散 ， 因 此 ， 在 热力 学 极限 下 ， 在 一 维和 
二 维 情况 对 积分 fd pn(p) 产生 发 散 的 贡献 . 在 有 限 温度 一 
维和 二 维 的 有 相互 作用 的 均匀 气体 的 BEC 也 不 可 能 发 生 [G71. 
当 存 在 简 谐 势 场 ， 由 于 态 密度 所 展示 的 不 同行 为 ， 热 涨 落 
的 影响 被 强烈 地 抑制 . 事实 上 ， 对 于 均匀 气体 态 密度 ple) 
RAMA O°, CERAM, CAME, AE 
ple) < e”!， 因 此 积分 (2-28) 式 在 2D 情形 也 是 收敛 的 . dH 
应 的 临界 温度 由 下 式 给 出 


1/2 
kl = ron] (2-34) 


式 中 ozp -(0,0,)". 人们 首先 注意 2D 情况 相应 于 
N >œ Ail o) 一 0,， 但 积 No3p 保 持 为 常数 的 热力 学 极限 . 
为 了 在 实际 3D 势 阱 中 获得 2D BEC, 应 该 选择 第 三 方向 的 
频率 w- 足 够 大 使 条 件 hos << korp e ho, 得 到 满足 ， 这 意 
味 着 对 阱 的 变形 有 严格 的 要 求 . 简 谐 势 阱 中 2D BEC 的 主要 
特性 和 对 非 均匀 气体 的 拓展 已 被 详细 讨论 ‘63. 

在 VD 情况 下 ， 即 使 有 简 谐 热 存在 也 不 可 能 发 生 BEC, 
因为 在 式 (2-28) 中 的 积分 是 对 数 发 散 . 这 意味 着 如 果 积 
No 保持 为 定 值 ，1D 的 BEC 的 临界 温度 将 趋 于 零 . 事实 
上 1D 情况 下 理想 气体 的 临界 温度 可 被 计算 I 
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kaTip = hey, PUT (2-35) 
式 中 心 = o.. 尽管 在 热力 学 极限 下 不 存在 BEC， 然 而 对 
有 限 粒 子 数 系统 在 合适 的 温度 间隔 内 ， 最 低 单 粒子 态 上 显 
示 出 大 的 布 居 数 . 此 外 ， 如 果 的 值 和 势 阱 参数 选择 适当 ， 
可 以 观察 到 与 最 低能 态 宏观 占据 数 有 关 的 有 趣 现 象 ， 它 发 
生 在 两 个 不 同 的 阶段 . 势 阱 的 有 关 参 数 同时 满足 条 件 
Tin < Ty fll hoo, < ka73y,， 这 里 7 与 (3.26) 式 给 出 的 临界 温 
Ei, o 是 地 平面 内 势 阱 的 频率 . Ty « T < Ty 间隔 
内 ， 仅 径 向 自由 度 内 凝结 ， 而 在 轴 向 自由 度 内 没有 凝聚 发 
E. 在 7< Ty, 的 更 低温 度 轴 向 也 开始 凝结 ， 并 且 基 态 全 面 
进入 宏观 方式 布局 两 段 BEC 的 例子 如 图 2-4 所 示 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 
TIT, 


图 2-4 EAR AME BE Co, = 56 x10* 0.0 PN = 10° 
个 原子 的 理想 气体 的 行为 

实 线 : BS (n, = 0.n，= 0.n. = 00. 上 原子 的 比例 

ER: 最 低 径 向 态 Cn, = 0.n，= 00 上 康子 比例 


LIA E IYE HEL BER WK JE R A A JI SE AN YE EAE! 

对 十 低 维 势 阱 情况 的 研究 Wl 表明， 如 果 低 维 势 阱 不 是 

“ 硬 壁 "的 ， 而 是 空间 渐变 的 势 ， 则 在 2 维 情况 下 ， 任 何洁 

函数 势 阱 中 都 可 以 形成 BEC; 而 在 一 维 情况 下 ， 则 只 有 在 
TER BUNT 2 的 势 阱 中 才 有 可 能 形成 BEC. 


2.6 非 简 谐 势 阱 中 的 BEC 


在 实现 BEC 实验 中 要 达到 所 要 求 的 低温 的 关键 步骤 是 
燕 发 冷却 ， 这 个 技术 的 原理 是 根据 阱 中 热 原子 的 损失 而 降 
低 湿度. 如果 人 们 能 够 以 可 道 的 方式 绝热 地 冷却 系统 Po7n 
将 会 打开 新 的 有 趣 的 前 景 . 通过 以 一 个 比 内 部 平衡 速率 更 
慢 的 加 率 绝热 地 改变 势 阱 的 形状 可 以 获得 气体 的 可 逆 冷 却 
用 于 研究 绝热 变化 影响 的 一 类 囚禁 势 有 如 下 笃 律 形式 

Vex (t) = Ar? (2-36) 
为 简单 ， 这 里 假定 是 球 对 称 的 .在 上 式 势 阱 中 BEC 的 临界 
温度 已 求 得 


143/246 
er -| NE 6n A? Soe 
PE | Om rar5zG72+5) 


(2-37) 


这 里 已 引入 了 参数 5 = 3/a ， 而 rOREKHUMAAK. 通 
过 设 5 = 3/2 和 4= mo; /2 我 们 可 以 回 到 各 向 同性 简 谐 势 阱 情 
形 的 转变 温度 . 设 5 — 0 得 到 箱 势 的 结果 . 在 (2-36) 束 缚 势 中 无 相 
互 作用 气体 的 热力 学 可 以 直接 计算 外 忆 例如 可 求 得 凝聚 比例 
NoIN aeri, ”与 可 逆 过 程 更 密切 的 就 是 粹 ， 它 在 绝热 
过 程 中 保持 不 变 . 在 以上， 系统 可 近似 为 经 典 的 麦克 斯 书 - 玻 
尔 兹 曼 气 体 ， 单 粒子 粹 有 如 下 形式 
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Perm s) T (2-38) 
JA ix 47; Fe ap DUE E DS E T/T, KF ah HY 
数 4. X FEOXSET EERE (Co 给 定 )，4 的 绝热 变化 ( 例 
如 简 谐 势 阱 的 绝热 膨胀 ) 并 不 能 更 接近 相 变 ， 因 为 比值 T/T 
保持 常数 . 但 通过 绝热 地 增加 可 使 比值 T/T 降低 ， 例 如 
从 简 谐 势 6 =3/2 变 到 线性 势 0, =3 ， 可 得 到 初 末 约 化 温度 
t=T/T, ZIM KAR = 0.70 人 .在 这 种 情况 下 一 个 在 两 倍 
临界 温度 ( 约 化 温度 t1 =2 ) 的 系统 可 被 冷却 到 接近 临界 点 
( 约 化 温度 1,=1.1). 利用 这 种 技术 (适当 地 改变 5 ) 应 
能 使 一 个 系统 从 没有 凝聚 的 高 温 相 绝热 地 冷却 到 有 和 较 大 比 
例 的 原子 在 凝聚 态 的 下 以 下 的 温度 . 应 用 这 种 绝热 变形 实 
现 BEC 的 可 能 性 最 近 已 成 功 地 进行 了 实验 探索 021. 
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第 三 章 排斥 相互 作用 对 BEC 的 影响 


3.1 引言 


在 Bose-Einstein 凝聚 实验 中 ， 中 性 原子 有 磁场 势 阱 作 
用 ， 原 子 间 有 有 效 的 排斥 或 吸引 相互 作用 ， 这 与 原子 的 
Bose-Einstein 凝聚 理论 好 针对 自由 中 性 原子 ， 也 没有 原子 与 
原子 间 的 相互 作用 ， 仅 仅 是 自由 原子 在 动量 空间 的 凝聚 不 
一 样 . 将 Bose-Einstein 凝聚 理论 推广 到 外 势 场 木 为 零 的 情 
况 并 不 困难 ， 但 原子 间 存 在 相互 作用 时 ， 一 般 来 说 不 易 求 
AR. 在 总 粒子 数 N 很 接近 于 凝聚 体 原 子 数 NS. B] 
N - No << No 时 ,有 Bogoliubov 近似 处 理 与 微 扰 法 4731， 当 
这 条 件 不 满足 时 ， 这 些 近似 与 微 扰 法 0 也 就 不 能 应 用 . 本 
章 的 目的 在 于 探求 不 受 N - No << No 条 件 限 制 且 g>0 的 中 
性 原子 的 Bose-Einstein 凝聚 的 一 般 方法 04( 下 面 简称 BEC 
一 般 方法 ) 为 了 简化 并 能 解析 求解 ， 并 应 用 于 势 场 V =0， 
但 有 排斥 相互 作用 g >0 情 形 的 中 性 原子 体系 . 求 其 解析 
解 ， 并 进行 玻 色 - 爱 因 斯 坦 凝 聚 的 有 关 计 算 ， 分析 原子 间 
的 相互 作用 对 玻 色 - 爱 因 斯 坦 凝聚 的 影响 . 


3.2 近 理 想 Bose 气体 的 理论 描写 
通常 用 近 理 想 气 体 来 描写 对 理想 气体 改进 的 物理 气体 
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模型 ， 它 是 非常 稀薄 的 粒子 系统 ， 粒 子 问 的 相互 作用 是 在 
粒子 间 的 有 限 范围 的 作用 势 ， 不 存在 二 粒子 束缚 态 . 气体 
兢 注 的 意思 是 能 将 粒子 间 的 相互 作用 作为 理想 气体 的 微 扰 . 
在 非常 低 的 温度 下 ， 这 种 系统 有 两 个 重要 的 长 度量 纲 的 参 
数 ， 热 波长 4 和 平均 粒子 间距 V .这 两 个 参数 在 大 小 上 是 
可 以 比较 的 ， 但 它们 一 定 要 远大 于 粒子 间 势 场 的 作用 范围 . 
在 量子 力学 中 ,一 个 粒子 不 可 能 定 域 在 它 的 德 布 罗 意 波长 
范围 内 ， 这 里 它 可 用 热 波长 来 代替 . 这 样 一 个 粒子 “展开 ”在 
比 势 作 用 范围 大 得 多 的 空间 范围 上 ， 在 任何 给 定 粒子 的 相 
互 作用 范围 内 找到 另 一 个 粒子 的 几率 是 很 小 的 . 因此 ，@ 尽 
管 粒子 之 间 的 势 可 以 有 很 大 的 值 ， 粒 子 所 感受 到 的 有 效 相 
互 作用 是 很 小 的 ; @ 相 互 作用 势 的 细节 不 重要 ， 因 为 展开 
在 空间 中 的 一 个 粒子 仅 感受 到 势 的 平均 影响 . 在 量子 散射 
理论 中 ， 通 过 一 个 势 场 产 生 的 低能 散射 不 依赖 于 势 场 的 形 
状 ， 而 依赖 于 来 自 势 场 的 单一 参数 ， 即 散射 长 度 互 ， 低 能 
散射 时 总 的 散射 截面 为 4x52?. 因此 二 可 大 致 地 认为 是 势 的 
有 效 直径 ， 也 可 以 说 对 于 低能 情况 来 自 一 个 势 场 的 散射 看 
起 来 像 来 自 一 个 直径 为 a 的 硬 球 的 散射 因此 ， 在 极 低 的 温 
度 下 可 借助 三 个 参数 4，V' 和 来 描写 这 种 非 理想 气体 
这 样 一 种 体系 的 哈密 顿 量 可 由 一 个 有 效 的 哈密 顿 量 代 
替 ， 这 个 哈密 顿 仅 有 散射 参数 明显 地 出 现 . 在 硬 球 势 模型 下 
利用 性 势 方法 求 得 在 低能 (4 > 0 ) 球 对 称 (s BOBO KE a 
s 波 相 移 5 之 间 的 关系 [371 


= .. tan ó, 
a = - lim 一 一 
k20 大 


对 十 硬 球 势 ， 散 射 长 度 就 是 硬 球 直径 . 对 于 慢 粒 于 在 球 方 


(3-1) 
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势 阱 的 散射 求 得 中 


tan dy = zl 


tan Kd © ] (3-2) 


Kd 

式 中 K2 =k? +K, KÈ = 2ml /加 ，W 是 阱 深 d 是 阱 宽 . 如 
将 上 一 0 的 波 函 数 表示 为 硬 球 势 波 函 数 的 形式 ， 散 射 长 度 
与 8 波 相 移 的 关系 仍 是 (3-1) 式 . 散射 长 度 和 与 势 阱 参数 
和 内，d 之 间 的 关系 可 表示 为 


i- 4m Kod E ] (3-3) 


当 势 阱 深度 Kod > x/2Wf, m QG-2)X&/4b tanóo <0, MA>0, 
表现 为 排斥 势 . 当 势 阱 的 深度 不 够 ， Kod < z/2 时 ， 阱 内 束 
缚 不 住 s 态 的 粒子 ，tan6o。 >0， 故 散射 长 度 是 负 的 ， 表 现 为 
吸引 势 . 

近 理 想 Bose 气体 的 条 件 为 : a(NIV)?D <1. 仅 考虑 
二 体 作用 的 玻 色 粒子 系统 的 哈密 顿 量 为 


H =Y Eaa, +E E< pi PlU | pp napa, 

(3-4) 

求 和 是 对 下 标 出 现 的 所 有 动量 . 粒子 的 产生 和 潭 灭 算 符 服 

A XA NUM à, -ata, =1. 假定 粒子 动量 很 小 ， 即 

pa /1h << 1， 以 及 激发 态 上 的 粒子 数 相对 基态 上 的 粒子 数 很 
小 Ml ag dg = No < N， 可 求 得 哈密 屯 量 的 表达 式 为 (35) 


N?U p! NU 

tnt “0 ata $ 0 A^ atat A^ 

H= - IES asd, + È lapan + apa", +24}4,) 
pom 21 ps 


2] 


(3-3) 
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根据 微 扰 理论 二 级 近似 公式 和 散射 振幅 的 玻 恩 近似 公 
式 可 求 得 


_ 4nh? a q+ Anh^à my eo) (3-6) 
p 


Uo = 


通过 算 符 的 线性 变换 使 哈密 顿 量 对 角 化 ， 从 而 求 得 能 级 . 对 
于 基态 能 量 有 


(3-7) 
在 7=0 时 , 气体 的 化 学 势 为 


OE,  4nh)?aN ， 32 |a^N 

“ev mV ECT nV 
a cha PT PARES 为 使 ge 和 /为 实数 ， 则 要 求 散 
射 长 度 5 > 0. 这 表明 上 述 公 式 仅 适 用 于 散射 长 度 互 > 0 的 
空间 均匀 的 Bose 气体 原子 体系 . 

为 研究 相互 作用 对 BEC 的 影响 ， 许 多 文献 研究 了 这 种 
均匀 近 理 想 气体 的 转变 温度 相对 于 理想 气体 转变 温度 To 
((2.8) 式 ) 的 变化 AT IT, - (T, -70)/To。， 得 到 非常 不 同 的 结 
果 . i ee 因为 空间 的 排斥 导致 动量 空 


] (3-8) 


间 的 吸引 ， 这 这 种 非 理 想 气 体 更 容易 凝聚 . Toyoda UC! 
“hee. 
aA 8V2r =3\1/6 
(CAT /To )royoda = Scop ) (3-9) 


El AT < 0. Gruter 中 根据 Monte Carlo 模拟 的 数值 计算 得 到 
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CAT /To )Gnter = 0.34(na*)°™ (3-10) 
最 近 , Holtmann U* 在 某 平均 场 假 定 下 给 出 
(AT1To)Homam = 070127) (3-11) 
K. Huang !7°) 应 用 状态 方程 的 维 里 展开 求 得 
sv2n vs 
AT IT.) nung = (3-12) 
OE = ag FEE He 


这 一 结果 与 Toyoda 的 结果 (3-9) 式 数值 相同 ， 都 正比 例 于 
Ya ， 但 符号 相反 . 而 式 (9-10),(9-11) 的 结果 的 数值 虽然 不 
同 , 但 都 近似 为 AT ~a 

对 于 束缚 的 非 均 匀 近 理想 气体 的 转变 温度 移动 AT, X 
献 [80] 求 得 的 结果 是 

ATIT = -13 -NS (3-13) 
ay 

有 限 原子 数 对 转变 温度 的 影响 (2.32) 式 也 产生 负 的 移动 ， 
临界 温度 的 测量 2 似乎 表明 有 负 移 动 的 出 现 ， 但 目前 实验 
的 精度 还 不 能 得 出 确定 的 结论 . 


3.3 Bogol iubov 近似 微 扰 法 与 BEC 一 般 方 法 


为 了 讨论 方便 , 我 们 参照 文献 [41], 给 出 基于 Bogoliubov 
近似 的 微 扰 法 概要 . ETER Hamiltonian Jj 


H = fa? rE*(- PN yat, gjd rere ee (3-14) 


HY, ¥* 4 Boson p f. HRE Boson 场 算 子 对 易 规 
JU. 现 考虑 温度 已 低 于 转变 温度 的 凝聚 相 ， 绝 大 多 数 原子 
为 凝聚 体 ， 只 有 少数 为 非 凝 聚 体 . 很 自然 将 多 表示 为 凝聚 
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tk o) C ROS EXER b OE TOf HK In (Bogoliubov 近似 ) 


È = P(r) + Gr) (3-15) 
(3-15) 式 成 立 的 条 件 为 N -No << No. Yr), à 可 表示 为 


(p) = VT erst ft) 


> f (3-16) 
G(r) = eS [u (a, - v ()a] 


(3-16) 式 即 Bogoliubov 变换 ，c, at 为 准 粒子 Boson HF, 3f 
满足 Boson 算 子 对 易 关系 ， 求 和 号 并， 不 包括 凝聚 体 . 
P(r) i EJERE Schrodinger 方程 即 Gross-Pitaevskii 方程 

[-Qm) à? v? - u + gel - 0 (3-17) 
当 求 得 凝聚 体 的 解 w(r) 后 ， 就 可 解 非 凝 聚 体 (激发 态 ) 的 
ujr), v,(r) 的 耦合 方程 

Lu; ELE v, = Eju, Q-18) 

Lv, -g|ivp u; =-E,v, 

L= “iv vso -u*2g|v | 

2m 

一 般 说 来 方程 (3-18) 有 正 能 解 u(r), v e), E, *3 ff fib f 
ur). vir), -E,. TERKEPEN, RIERA. u0), 
v, C) ii e AF BJ IE RRA 

[du u ~v ve] = 8, (3-19) 


按 下 式 求 粒子 密度 的 统计 平均 
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Ji (63€ VI WE 38 gE R MI TE D D sh JI AG F EE 


«aja, >=< aja; > -6 = Ont, 


(3-20) 
f, [e^ -ap 


BE, 
HAE ER ALA EI. MRTA v=- sn, 


由 上 所 述 ， 微 扰 法 就 是 在 假定 已 有 No 个 原子 已 处 于 
BEC 状态 前 提 下 通过 Bogoliubov 近似 与 变换 得 出 的 计算 非 
凝聚 体 粒子 统计 分 布 n. = n(T) (<< No) 的 方法 . 而 我 们 这 
里 要 给 出 的 一 般 方法 是 包括 凝聚 体 与 非 凝聚 体 在 -起 求 统 
计 平 均 . 先 求 出 系统 的 完备 的 态 函 数 及 能 量 本 征 值 吧 ,E， 
然后 按 一 般 的 统计 力学 方法 J 计算 配 分 函数 必 ， 再 求 粒子 


Bn- 并 判明 处 于 基态 的 原子 数 No 是 否 已 发 生 BEC. 
这 种 理论 模式 与 无 相互 作用 自由 原子 的 动量 空间 的 BEC 81 
基本 相同 ， 差 别 只 在 g> 0 情形 多 ,,E, 的 确定 与 在 g = 0 情形 
FE. 8 = 0 情形 ， 我 们 只 需 解 线性 莅 定 证 方程 ,而 g >0 
情形 ， 则 需要 解 非 线 性 莅 定 庙 方 程 (NLSE). 对 于 V0 情形 
可 数值 求解 ， 对 于 天 =0 情 形 可 求解 析 解 ， 此 即 下 节 所 要 讨 
论 的 . 
3.4 有 排斥 相互 作用 的 中 性 原子 的 态 函 数 
首先 参照 (3-14) 式 ， 将 包括 势 阱 VY 与 原子 相互 作用 的 

ER: BW 

[rure r= vc A La ay ZA YW dr (3-21) 
式 中 风 ,更 为 原子 的 波 函 数 而 不 是 算 拓 ， 对 (3-21) 式 进行 变 
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jy 
Sf P'HY d r-o (3-22) 
并 满足 归 一 化 条 件 
je dr-1 (3-23) 


应 用 未 知 乘 子 法 ， 便 得 出 原子 所 满足 的 定 态 ( 包 括 基态 
和 激发 态 ) 非 线性 醉 定 谱 方 程 


* | 
civ +V «gp = Ev (3-24) 
m 


这 方程 形状 就 是 Gross-Pitaevskii 方程 3-17)， 但 含义 不 一 
样 . 方程 (3-17) 是 用 来 描写 基态 即 凝聚 态 的 原子 ，g = Nou 
wha 

(No 为 凝聚 态 的 原子 数 ，v= 22 7. zy kn) 而 方 
程 (3-24) 中 的 g = Nu, N 包括 基态 与 激发 态 在 内 总 的 原子 
数 ， 不 论 凝聚 与 否 . 同样 史 也 不 只 是 描写 基态 ， 也 包括 激 
发 态 在 内 . 式 (3-24) 帮 助 我 们 定 出 中 性 原子 系统 的 本 征 态 Y, 
及 相应 的 本 征 值 E, ， 而 这 些 正 是 我 们 计算 配 分 函数 4 ， 并 
RN = - 给 计算 原子 数 所 需要 的 ， EFRORER OBA 
斯 坦 凝聚 取决 于 总 的 粒子 数 N 及 温度 了. 这些 罗 是 归 一 化 
AY. 方程 (3-24) 是 非 线性 的 ， 当 g o 0 时 ， 这 些 w, 趋 近 于 完 
备 正 交 系 列 它 . 为 简单 起 见 ， 我 们 先 讨论 一 维 运动 的 自由 
原子 ， 这 时 方程 (3-24) 可 写 为 


pog 
(-— 


a grep = Ev (3-25) 


+ 42° 


Ji (5 XE VAN P RE K JE IR 0 2 JI "RC VER UE 


4g 0n. Jr; FEG-25)u[?j Jy 


2 


+KP, =0 (3-26) 


d 
C 
dx 


方程 (3-26) 满足 边界 条 件 “ 当 =0. a, Y -07 (SEM IER A 


列 为 
~ er 
Y, = Esn kxk, = me 
a a 


= (hk, )* 


"n 


(3-27) 


2m 
MEFA, SY =N aY, NA Hal, YH 
无 量 纲 的 ， 方 程 (3-25) 可 写 为 


2 


d^ 2 2 2 
Ga 2 + (ky + 2k) - 2k 2 yp, = 0 (3-28) 
" n ~ >= EN: =. 
RP k= us Na-kLNa Niv2dx=1 . Ng 
0 


化 常数 . 换 变数 并 引入 新 参量 ， 将 (3-28) 式 写 为 


z 


G tla- ary, =0 (3-29) 
k? 
Rha’ = pGAlaeAlx 
or 
(3-29) 式 的 积分 为 
day, y jure a 
We] ed AT, (3-30) 
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v, 


" Í dY, 
640 PI 2279) 
这 是 第 KME. uH Jacobi MiP BOK RU 
V, - sme A) 
sn(E + dnk +4m'i KA?) = sn(£, A?) 


1 m2 
s dv, ? d 
人 (4 ) = ki ? k) a £ 2 
oya- rI- è disais! p 


K (27)=K(1~ 2?) 


(3-31) 
式 中 mm HR. MA qoe II, quj sé A!) 可 写 为 


L 
n& oq tne 
Yt sin CL (3-32) 
3 


3 2 
sn(č A) = 一 一 
AK Sle 2K 


SLIP SE SEM AR, POR ATI, EX 


Kad 1 | mg [4Na 

go g Na 
n, = = — = 3-33 
mrs ro ET HM 


FOUN 为 总 原子 数 ， 由 (3-28)，(3-29) 式 中 参数 的 定义 得 


22 
kia Hg = (ky, ss pa 

e Nac Va = n; Na (3-34) 
x A hea ul fa 


为 求 得 Na TAM ROA, UB GE 义 (3-28) 式 出 发 . 
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并 应 用 (3-31) 式 得 
z Jh tiia 
Na = —= (3-35) 
f EHARA Ne Ade 


0 


根据 x =0, a, = 0 的 边界 条 件 ， 故 J 妇 + 尼 ao=2nK ， 则 上 
式 为 


Na = Ink an A 二 一 大 ah (3-36) 
LE sn? (£,43)d€ Lh sn? (£, 3! )d£ 
k 一 
由 边界 条 件 ks + 局 a= = ak R- D Na 得 
ng — 2KÀ 
— =- (3-37) 
^ è VNan 
th (3-29), (3-34) 式 
EE 一 
k, = k,,k, =k, V Na (3-38) 


(3-36) 式 给 出 Na 与 的 函数 关系 Na = Nala). 代 人 (3-32) 式 右 
端 ， 当 左 端 n,,n 给 定 后 ， 由 (3-37) 式 便 能 定 出 X 的 本 征 值 
4,， 并 可 计算 出 ,及 本 征 值 , 以 及 能 量 本 征 值 


2 
^k) 20k! af (Inka? 
SY Cn PI PN E 


2m 2m 2m a 


"E (nn) (1 42) 
p , 


(3-39) 
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图 3-1 给 出 当 n, 给 定 后 ， 本 征 值 随 nn 的 变化 曲线 . 由 这 
曲线 可 读 出 本 征 值 元,45,4;,…. 图 3-2、 图 3-3 分 别 给 出 
En, E, BE n RIRA. 


当 4 很 小 时 ，K 5 二， 由 (3-37), (3-38) 及 (3-33) 式 给 出 


A 
ng = ——n 
i = (3-40) 
k, e Va = EAL my V8 - ™ 
A a | Na a a 


X (3-40) HIE Y 4 2> 0E] g >On, k, HF 3-27) M 

k= 了 二 .这 时 5, 趋 于 (3-27) 式 的 无 相互 作用 自由 原子 的 能 量 
a 

AREAL E,. 


图 3-1 AREE H An Bi n h 4e fe hh Ac 
ny = 0.5.2.5.10 C. fr tl £X, sh 
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图 3-2 Anh üETÉ Du £, Bt n (0) AE AE dh R, 
ng =0.5,2,5,10 已 在 曲线 上 标 出 


图 3-3 ARKEA En Bi n 的 安 化 曲线 ， 
My = 0.5.2,.5.10 已 在 曲线 上 标 出 . 
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3.5 中 性 Bose 原子 的 BEC 计算 
参照 文献 [81] Bose 气体 的 基本 公式 即 配 分 函数 5、 粒 
TUN 及 能 量 E 分 别 为 
&--Y o, n - e ^7) 
no-— go. 2e. G-41) 
l a e p 
将 此 公式 应 用 于 方 盒 (x=0,aiy = 0,a;z = 0,a ) 中 的 无 相互 作 
用 自由 原子 气体 B0， 便 得 
€ = -ob Y, Inl- e 7 P ot, (3-42) 
Par 
Xs Eu, = H) pg ng oA i TEAR 
相互 作用 的 自由 原子 的 能 量 本 征 值 . 当 原 子 间 存在 排斥 相 
互 作用 时 ， 原 子 的 能 量 本 征 值 应 由 Ep Eg E R Epro 等 
由 (3-39) 式 给 出 ， 即 5, = Ene 而 相应 的 4,,4,,4, 等 由 
(3-37) 式 给 出 . 这 时 (3-42) 式 应 改写 为 


-a-o( prep tqrey+re;) hx? 
=-O0 Y Inle ^ 7 eee )es——— 
par=l 2ma° KT 


(3-43) 
在 作 凝 聚 计算 时 ， 十 分 重要 的 是 将 (3-43) 式 中 的 基态 
P=9=r=1 分 离开 ， 即 将 6 表示 为 
C= -wo In(1- e 4795 ) ~ wo p p(l - e 4 7 tn 
paral h 


(3-44) 
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和 式 中 不 包括 p=g=r=1 的 项 . &z-e7, g 


$o 47.75 
> -a(pe,*q'&, "£,—-3M| 
Ws eG (pz i5 X zo 0g TU tym) 

= - i - 
i-z 


2 CER ~ 
paredes e TU pH rnm) 


(3-45) 
式 中 N,N 分 别 为 基态 和 激发 态 的 原子 数 . 若 用 化 学 势 uR 


gU, Wa--fu, u«3E,, z-eP"» <1. 


参照 本 征 值 曲线 图 3-3 及 (3-45) 式 N = T N, = 


2 2 2 
“OIP Ey ty eg+r £,-36) 


e, Y  — e 


parel fos 


为 1 )， 可 算出 当 o=0.1 时 ，Wi,N;, 随 z 的 变化 曲线 为 图 3-4. 


(计算 时 ， 将 wo 取 


2 2 2 
e 97? £p*q Eq *r E-361) 


图 3-4 o = 0.05 f. MARAT FON, (KH, n, = 06,2506 
PEAR Vo bg) HUE AS LF ON, (虚线 ) 随 = 的 变化 曲线 


。49 。 


2000 年 上 海 太 学 博士 学 位 论文 


4x10" 
3x10 


= ox10° 


1x10° 


0 A 
0 200 400 600 800 1000 


图 3-5 z=1 时 , BORSE N BO YE Ae HER, 
ny = 0,2,5,10 c fri x HRB. 


946 0.2 0.4 0.6 0.8 1.0 
(CLAN 


图 3-6 = = 1 上 时、 凝聚 原子 数 与 总 原子 数 之 比 v_y w MiO Meh 
FRH n, = 0.2.5J0 的 曲线 ,虚线 为 1 - (7/7, ^ =1- (o, 70)" hik 


= 50° 


UEEATDEULUPTITDERAISEAI I 

Mcd. NERIS. 但 当 = 一 1，NM 和 急剧 上 升 ， 远 
远 超 过 N, ， 这 就 是 发 生 凝 聚 的 情形 . 为 求 得 系统 发 生 BEC 
的 临界 温度 藉以 及 g 对 的 影响 ,我 们 计算 了 当 =: 二 1 时 
N; 随 o-! 的 变化 曲线 ， 即 图 3-5. 从 图 中 看 出 ， 当 温度 
(Teo), HP HR MAIER g > 0, 使 得 激发 态 能 容 
纳 的 原子 数 Nt g =0 情 形 减 少 了 . 现 取 定 总 粒子 数 
N=16963 时 ， 即 图 3-5 中 虚线 所 示 ， 与 这 些 曲 线 的 交点 定 出 


后 按 凝 聚 到 基态 的 粒子 数 N.=N-N,， 对 (o/co.)'=(7/T.) 作 
图 ， 如 图 3-6 所 示 . 由 此 图 看 出 ， 与 理想 情形 (将 
sn =E =E, =I it, GAR PRS ARDRE) 
1-(7/7. 凡 相近 ， 但 有 差别 . 故 g 的 影响 主要 体现 在 临界 温 
BET. coz! g 的 增 大 ， 亦 即 , 的 增 大 ， 对 应 的 临界 温度 工 
也 随 之 增 大 ， 这 与 图 3-5 显示 的 g 的 增 大 ， 激 发 态 能 容纳 的 
原子 数 w; 的 减少 是 一 致 的 .这 也 与 文献 [80] 给 出 的 均匀 非 
理想 玻 色 气体 的 转变 温度 相对 于 理想 玻 色 气体 的 转变 温度 
的 变化 关系 (3-12) 式 随 相 互 作用 的 增强 或 密度 的 增加 而 提 
高 的 结果 是 一 致 的 . 


3.6 与 局 域 密度 近似 模型 的 比较 


近年 T.T.Chou 等 在 文献 [84] 中 采用 局 域 密度 近似 (local 
density approximation LDA) 来 处 理 势 阱 中 原子 的 BEC ju] 
题 ， 并 在 文献 [85] 中 引进 局 域 逸 度 CF) = zexp(-BV), 式 中 2z 
为 原子 逸 度 z = exp(w/K7),w 为 化 学 势 . ROG RB ABE UR 
子 气体 的 坐标 -动量 分 布 


ET 
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h^ p(F, p) =<< n, >>= 


(3-46) 


我 们 就 上 面 的 密度 动量 分 布 公式 与 我 们 求 得 的 结果 进行 比 
较 . 首先 在 我 们 的 计算 中 已 取 定 势 场 上 = 0 代入 局 域 逸 度 公 
式 , 得 6=s=e“. 参照 3-45) 式 ， 故 得 激发 态 的 动量 分 布 为 


SR 
-a(paep+gaey+r2er-3c0 


Ppar 7 式 中 不 包括 p=g=r=1 情 


2 2. 2 
rer” £p*q £q*r £j 36) 


形 . 为 简单 起 见 ， 我 们 取 9 =r=l pan n=2,3,4,…， 并 


4Zz-ze^, Mi 


ge^? Ej. ze" e 9 6D 
Pn = = - = 
T 1-2e79 ên 1-£e e 9" (En) 
d 7 (3-47) 
~ -on? 
ze 2 1 
(l-o? (e, -1) D 
1-Fe™ -z.e7 


将 (3-46) 与 (3-47) 式 比较 ， 得 on? (e, -1) = A2? ap ,并 应 用 (3-33) 
式 进一步 写 为 


2nh? za N 
5 A Pr. 
om (e, -1 KT 4N ai 
RE Ua D. MET AN ar PS EET (3-48) 
n, hen 
2ma^ KT 


当心 ,2 给 定 后 ， 便 可 根据 解析 解 求 出 e, ,并 按 (3-48) 式 的 定 
义 算出 R. 图 3-7 给 出 R 随 ,的 变化 关系 ,n 为 参量 ,由 图 看 
出 & 的 值 基本 上 在 3 附近 变化 .与 4 的 偏差 可 能 是 硬 球 模型 
中 系数 4425 的 导出 是 应 用 了 球 对 称 的 s 波 函 数 假定 9， 而 
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我 们 立方 体腔 的 本 征 函 数 并 不 兵 球 对 称 性 . 若 将 硬 球 半径 了 
修正 为 同体 积 的 正方 体 边 长 #5， 则 了 =(x/16)3a， 亦 即 用 了 
代替 (3-48) 式 中 的 a， 则 (3-48) 式 中 的 Rx3.2， 这 样 与 解析 
解 更 为 接近 . 但 严格 说 来 R 应 是 nn 的 函数 ， 即 R= Ra), 而 
式 (3-46) 和 (3-47) 原 本 均 是 一 级 近似 结果 . 关于 凝聚 体 动量 
分 布 公 式 的 比较 情形 较为 复杂 ， 主 要 因为 LDA 用 的 能 级 
hoy 是 声 子 谱 ， 与 凝聚 体 密度 p, 有 关 ， 亦 即 与 XT 有关; 而 
我 们 的 能 级 p?s, 仅 与 8 和 腔 尺 寸 a 的 数值 有 关 ， 与 温度 无 关 


图 3-7 n=2.510.20 时 ，RR 随 的 变化 关系 
3.7 小 结 


本 章 提 出 了 解 有 排斥 相互 作用 情形 中 性 原子 玻 色 - 爱 
因 斯 坦 凝 聚 的 一 般 方法 . 即 是 包括 凝聚 体 与 非 凝 聚 体 在 一 
起 求 统计 平均 ， 根 据 系 统 的 态 函 数 及 能 量 本 征 值 ， 按 一 般 
的 统计 力学 方法 计算 配 分 函数 ， 再 求 粒子 数 ， 由 处 于 基态 
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的 原子 数 No 判定 是 否 已 发 生 玻 色 - 爱 因 斯 坦 凝 聚 . 由 求 得 
的 波 函 数 和 能 量 本 征 值 的 解析 解 ， 估 算 了 这 种 非 理想 气体 
BEC 的 临界 温度 随 原子 间 相 互 作 用 大 小 和 原子 密度 的 变化 
情况 ， 得 出 排斥 相互 作用 增 大 ,将 使 得 激发 态 能 容纳 的 原 
子 数 wz 减少 ， 临 界 温 度 趣 升 高 ; 并 与 局 域 密度 近似 模型 的 
结果 作 了 比较 . 
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第 四 章 ” 简 谐 势 阱 中 有 排斥 相互 
作用 中 性 原子 的 BEC 


4.1 引言 


碱 金 属 原子 ( ”Rb, ”Na, Li) 气体 的 玻 色 - 爱 因 斯 坦 
凝聚 实验 的 实现 是 通过 激光 冷却 、 蒸 发 冷却 磁 阱 中 的 原子 
而 观察 到 的 ， 其 中 多 与 钠 的 散射 长 度 亏 为 正 ， 当 原子 靠 得 
很 近 时 ， 表 现 出 排斥 作用 ， 对 于 实现 稳定 的 凝聚 态 是 有 利 
的 ,而 锂 原子 的 散射 长 度 世 为 负 ， 不 利于 实现 稳定 的 凝聚 
态 . 有 关 的 实验 现象 为 对 在 势 阱 中 中 性 原子 实现 BEC 进行 
更 为 广泛 而 深入 的 理论 探索 提供 了 可 能 . 情形 也 的 确 是 这 
样 的 ， 几 乎 与 此 同时 就 有 关于 磁 阱 中 中 性 原子 所 满足 非 线 
性 醉 定格 方程 (NLSE) 定 态 解 的 研究 S71, 稍 后 不 久 又 有 NLSE 
非 定 态 解 的 研究 I， 这 些 本 应 成 为 BEC 研究 的 基础 ， 但 初 
步 看 来 ， 该 文 似乎 存在 本 章 中 将 提 到 的 一 些 问 题 ， 计 算 步 
又 也 的 确 复杂 些 . 实际 上 文献 [5] 作 者 也 已 指出 文献 [87] 的 
“THER RK, RARE”. 还 有 另外 一 个 重要 方面 ， 即 现在 发 
表 的 NLSE 的 解 主要 是 基态 ， 对 激发 态 则 很 少 有 报道 ,从 全 
面 理解 BEC 的 形成 和 发 展 过 程 ， 对 基态 研究 固然 重要 ， 对 
激发 态 研究 也 是 不 可 缺少 的 ， 这 对 凝聚 态 的 稳定 性 研究 有 
重要 意义 . 在 本 章 中 给 出 一 种 数值 求解 NLSE 的 方法 ls， 讨 
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iE AE AR HE aE AS AB ae HT FE AO ie eR BOW CSA {ele eda 
给 出 简 谐 势 阱 中 中 性 原子 的 NLSE 的 基态 和 激发 态 解 ， 对 计算 
精度 进行 了 分 析 ， 与 大 N 极限 下 基态 能 量 进行 了 比较 . 


4.2 非 线 性 茶 定 谱 方 程 及 其 数值 求解 方法 


REGSE, FRES V) = 了 mor? 中 的 中 性 
原子 体系 ， 其 有 效 哈密 顿 为 


1 2d Qnh?a 


H- [e + + 


[PPY (4-1) 


为 了 导出 凝聚 波 函 数 满足 的 方程 ， "—— 
算 子 多 的 时 间 演 化 方程 为 


ov 
in - v, H]- 
X UI 
R (4-2) 
-X iy ul Tm «uev 
2m 


其 中 U,Q-4nf/a | m ,在 绝对 温度 T-0, 3ECGUS TET V(r) 就 
是 凝聚 体 波 函 数 ， 仍 用 风 (r) 表示 归 一 的 凝聚 体 波 函数 ， 则 
得 到 

LAG n? 


2 15 2323 
^r oa Vr.) 2 nar P(r,t)+ (4-3) 


NUol¥ (rf n) 


其 中 m 为 单 原子 质量 ，wy 为 简 谐 势 角 频率 ，AN 为 凝聚 体 的 
WU PA. 这 个 方程 就 是 Gross-Pitaevskii 方程 ， 即 非 线性 薛 
AE FN R(NLSE), © E 1961 年 由 Gross 和 Pitaevskii 分 别 
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独立 导出 的 .实际 上 第 三 章 已 经 通过 能 量 泛 函 的 变 分 导出 
了 这 个 方程 . 

为 求 定 态 解 , REUSE UU P (是 体系 的 化 学 
势 ) 代入 (4-3) 式 可 得 吧 (r) 满足 的 方程 


h? 
= awe) + moles) + NUP P(r) = ut) 
m x 


(4-4) 
对 球 对 称 的 s 波 函 数 ， 方 程 (4-4) 可 以 简化 ， 采 用 没有 相互 


uva 
作用 情形 简 谐 振子 基态 长 度 单 位 Em 和 谐振 子 能 量 


单位 iwy 进行 归 一 化 ， 令 : 


ro = (4-5) 
a ee 
2mo, | 
1 2 
NU, a EPO qw 
umm) oo 
4n 
2mo, 
式 中 
2Nà 
Ci = (4-7) 


与 文献 [88] 的 非 线性 常数 Ci = 2NaQmo, IAO 是 -- 致 的 . 
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于 是 由 式 (4-4) ~ (4-6) 得 


dp ^ peo -0 —— (48) 
dx^ 4 x 
归 -化 条 件 为 : 
amfer dr= ASTE =] (4-9) 
0 0 


式 中 各 量 均 是 无 量 纲 的 ， 而 (4-3) 式 中 的 NU, 等 也 包含 到 参 
量 Co 中 去 了 . AC > 0 时 ，(4-8) 式 就 是 简 谐 振子 方程 ， 其 


RAE A = n+ 二， 考虑 到 PO fE x -0 附近 应 为 有 限 的 ，n 


只 能 取 奇 数 ， 1 =1.3,5,7,. …，hn=1 为 基态 ，n =3,5,7,… 为 激 
发 态 . 若 Co 不 为 零 且 给 定 后 ，p 就 是 定 态 NLSE 方程 (4-8) 
的 本 征 值 ， 它 的 取 值 应 能 保证 : (i) 当 x > oft, D(x) 30; 


(ii) 4 x -> On, Y Gc CO. arg. 在 x > OME, D(x) 
x 
可 表示 为 : 


P(x) = O(0)x (x > 0) (4-10) 
在 实际 数值 求解 方程 (4-8) 时 ， 我 们 将 它 写 为 两 个 一 阶 方程 

do 

dx ^? 

E. o? (4-11) 

aa prac, Toyo 


参照 (4-10) 式 写 出 边 值 条 件 为 
(i) 34 x > oft 
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dx) = 0, y(x)| =0 (4-12) 


xo 


(ii) 在 x -> 0 附近 
PD, =P Oe. yo =O) — (4-13) 


当 9'(0) 给 定 后 ， 就 可 按 Runge-Kutta 方法 由 x=0 到 
x 一 oo 进行 积分 ， 若 1 WRAY, 34 x o0 时 方程 (4-9) 的 
积分 是 发 散 的 ,只 有 当 p 选择 恰当 时 ， 才 会 使 得 (x) 在 
x — ool] SCESE. 还 要 注意 到 ， 既 使 COCz) 在 x 一 oo 时 收 
AFTE, IA B(x) 也 不 一 定 满足 归 一 化 条 件 (4-9)， 为 了 
满足 归 一 化 条 件 还 要 适当 选择 边界 值 w(0D); o(0'3p 并 
不 存在 文献 [87] 给 定 的 关系 式 (0) - BAN y +1) (这 
里 NAPY 即 我 们 方程 中 的 参量 C。) ， 因 为 解 (x) YE x — oo 
处 收敛 于 零 并 满足 归 -- 化 条 件 就 已 将 OWS B. 完全 确定 
了 ， 就 不 可 能 满足 如 文献 [87] 给 定 的 上 述 关系 . 下 面具 体 计 
算 后 还 要 分 析 这 个 问题 . 


4.3 基态 解 


我 们 利用 前 面 的 方法 计算 了 基态 波 函 数 和 基态 能 量 ( 本 
征 值 6 )， 在 实际 计算 中 将 非 线 性 参数 C, 的 取 值 范围 定 为 
0.1~150， 上 略 大 于 文献 [88] 的 取 值 范 围 0.1~100. 求解 精度 为 
eps =10“ ， 按 (4.12) 及 (4.13) 式 即 


49? o) + y? (x) | s « eps (4.14) 


je coax -1i « eps (4-15) 
0 
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所 求 得 的 基态 波 函 数 的 数值 结果 如 图 4-1 所 示 ， 横 坐 
标 以 谐振 子 基 态 长 度 为 单位 . 可 以 看 出 : 基态 波 函 数 的 分 
Ay BE. Cy 的 增 大 逐渐 变 宽 为 超 高 斯 型 ， 对 Co = 100 进行 数值 
拟 合 ， 其 基态 波 函数 近似 为 025xexpf-(x*14.5734}， 此 时 的 波 
函数 的 宽度 为 4~5 谐振 子 基态 长 度 单位 . 按 文献 [2] 给 出 的 
实验 参数 ,对 Na 原子 , m = 3.841 x10 75 kg, 0, 221x235 s, 
则 谐振 子 长 度 单位 为 0.96pym， 凝 聚 原子 的 实际 空间 分 布 约 
为 2x0.96x4.5x10 um ， 按 文献 [2] 给 出 的 势 阱 体积 10-8 cm, 
可 算出 阱 的 大 小 为 10 em? ) ~20pm， 二 者 可 以 比拟 . 基 
态 能 量 随 C, 的 变化 曲线 如 图 4-2 中 实 线 所 示 ， 当 Co 增 大 ， 
即 凝聚 体 的 原子 数 增多 时 ， 其 体系 的 单 粒子 本 征 能 量 随 之 
增 大 ， 这 是 凝聚 原子 之 间 的 排斥 作用 的 反映 . 


图 4-1 非 线 性 常数 ,分 别 为 0.11.5.10.25.50.100.150 
的 BEC 基 念 波 函 数 ( 按 宽度 增加 的 顺序 ) 


基态 能 量 的 数值 拟 合 表明 它 与 非 线 性 系数 有 如 下 关系 
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或 写成 


有 = (Ce +1.5 2y5 (4-16) 
í in EA VIE 

WN) = hol Nā Mh | 2mo) ? 1,552) (4-17) 
8 

6 

ad 

2 

P pera EEE EE 

0 40 80 120 160 200 

e 


0 


图 4-2” 基 念 能 量 随 非 线 性 常数 Co 的 变化 
实 线 为 数值 解 结果 ， 虚线 为 (4-17) 的 计算 结果 


现在 分 析 文 献 [87] 的 算法 ， 他 们 所 用 的 边界 值 是 
BD'.(0)= YBAl+Co) ， 而 我 们 的 @B'(0) 是 作为 数值 求解 的 调 
整 参数 使 得 波 函 数 归 一 化 得 到 的 ， 为 了 与 文献 [87] 的 结果 进 
行 比较 ， 我 们 在 图 4-3 中 给 出 了 对 每 一 个 非 线 性 参数 Co 求 
得 的 定 态 解 的 (0) 与 Co 的 关系 曲线 ， 图 中 用 虚线 来 表示 ， 
按 文 献 [87] 公 式 wme(O=VZMI+Co) 算 得 的 四 。(0) 与 Co 的 关 


系 , 图 ! 


明显 ， 


让 用 点 线 表示 ; 以 及 比值 2'(0)/w' .(0) 用 实 线 表 示 . 很 
当 CGo 很 大 时 ，@'(0)/@' (0)->1， 文 献 [87] 的 公式 


DBD'.(0) 可 用 ， 但 当 Co 较 小 时 ，@2' (00102' (0) 远 远 偏离 于 1， 
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D'.(0) 是 不 可 用 的 . AHT d"2(0). YU ER BAS RT Be TCR UA 
-化 . 


图 4-3 i (E (0) 随 非 线性 常数 Co 的 变化 
虚线 是 我 们 计算 中 使 用 的 边界 值 D (0) ， 
点 线 为 文献 [87] 中 的 边界 值 @ ' (0) = JB I0 C,)- 
实 线 为 二 者 之 比 


现在 就 文献 [87] 中 图 1 的 参数 p= 人 -= 43 进行 实例 
10, 


计算 ， 比 较 我 们 的 与 Edwards Burnett (下 面 简 称 为 E.B.) 的 
计算 方法 . 按 我 们 的 方法 ， 根 据 上 面 图 4-2 的 B 与 Co 关系 
曲线 及 进 … 步 的 计算 可 得 , 当 =4.3 时 , Cy = 34372365, UE 
函数 B(x) 随 x 的 变化 如 图 4-4 中 的 实 线 所 示 ， 满 足 归 一 化 
JO? (x) dx = 0999993 . 按 E.B. 方 法 , e. 0) - JB C.) 及 
TE R E x AAR ABT, RIC, = 20.33642 ， 波 函数 
D(x) Bi x WEL AP 4-4 EAR, URN EUH E ART, 
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图 4-4 本 征 值 8 = 4.3 BEC HAR RH 
RRA D(x), BRAD. (x), 
ARIRIH SCIC D(x) (与 实 线 重 合 ) 
而 j@2(x)dx=1.69018. 但 这 只 是 E.B. 方 法 的 第 一 步 ， 即 找 


E = 7 9. Aor) ， 为 了 归 一 化 ，E.B. 将 四 修改 为 


四 = 和 we(No Aso). A = [o a 即 文献 [87] 中 的 (4-11) 
1 
X. 于 是 


esjrertar- (4. Jo? 94x = 

$ (4-18) 
B fo2(x)dx=1 

No C 

m, Ge orem. 代入 人 (4-18) 式 是 满足 的 . 波 


1 


T£ (4-7) 3K, 
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函数 G6.) LO) ER 4-4 中 用 点 划 线 表示 ， 与 我 们 


来 得 的 oo S5. 由 此 看 出 EB 方法 引用 不 适当 的 边界 条 
fox yo — ， 没 有 达到 月 一 化 ， 而 在 第 一步 中 得 到 修正 
很 明显 E.B 方法 通过 计算 一 个 参数 为 No 的 系统 ， 反 求 得 参数 为 

N 的 系统 的 解 ， 即 求解 用 到 的 参数 N 并 非 实际 系统 的 参数 N, ， 
而 N 在 计算 完成 之 前 是 不 知道 的 . 车 Ni 给 定 ， 则 需 计算 许多 No 
去 试探 ， 看 哪 一 个 与 实际 要 求 的 N 相近 ， 而 选择 该 No 试探 中 的 
对 应 解 . 而 我 们 的 方法 就 是 从 实际 系统 出 发 求 本 系统 的 解 ， 即 根 
据 给 定 的 Ni 由 (4-7) 式 求 出 Cu ， 然 后 在 波 函数 满足 归 一 化 及 在 
为 很 大 时 收敛 两 个 条 件 下 求解 方程 (4-11) 直 接 得 出 本 征 值 p 
及 相应 的 波 函 数 ， 故 不 存在 上 述 问题 


4.4 激发 态 解 


线性 化 的 GP 方程 (即将 含 时 的 GP 方程 经 过 Bogoliubov 
变换 得 到 耦合 方程 ) 的 解 是 凝聚 体 的 简 正 模式 ， 也 称 为 宏观 
量子 系统 的 元 激发 和 准 粒 子 . 集体 激发 则 是 服从 流体 动力 学 
方程 的 大 范围 的 密度 涨 落 . 这 种 准 粒子 的 寿命 由 形状 振荡 的 
阻尼 来 确定 ， 实 验 上 的 阻尼 时 间 约 为 230ms， 目 前 还 没有 理 
论 期 望 值 . BEC 元 激发 的 研究 的 最 终 目的 是 集体 激发 谱 的 完 
整 描 述 ， 包 括 准 粒子 寿命 ， 和 不 同 温度 与 高 激发 能 量 条 件 下 
的 行为 . BEC 的 激发 起 因 是 密度 涨 落 或 捕 陷 势 的 含 时 微 扰 . 无 
相互 作用 的 理想 气体 所 有 模式 的 频率 都 是 简 谐 势 频 率 vo 的 整 
数 倍 . 在 均匀 弱 相 互 作 用 密度 为 m 的 玻 色 凝 聚 体 中 ， 最 低频 
率 的 激发 使 是 在 7 = 0 以 声速 * = SnU 7m 传播 的 声 子 . 在 
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非 均匀 凝聚 体 中 最 低 的 激发 频率 正比 于 mw ， 与 mm 无 关 ， 内 
为 声速 和 凝聚 体 大 小 对 m 的 依赖 关系 恰好 相 消 .对 各 向 同 
性 的 简 谐 势 最 低 的 简 正 模式 的 频率 为 V2v !G1， 对 于 大 于 平 
均 相 互 作 用 能 的 激发 频率 ， 预 计 会 发 生 从 集体 到 单 粒 子 跃 
迁 的 行为 . 我 们 下 面 将 通过 定 态 的 GP 方程 给 出 原子 在 简 谐 
ADIN BAA LINGER, TERRI Rit P BEC 激发 中 
部 分 地 处 于 这 种 态 的 几率 是 存在 的 . 
根据 前 面 的 数值 解法 ， 对 球 对 称 的 激发 态 作 类 似 的 计 
算 ， 对 于 不 同 的 非 线 性 常数 Co 分 别 求 得 了 n=3,5,7,9,11,13 的 
HESR DOM p, 图 4-5 给 出 了 非 线性 常数 Co 分 别 为 
0.1,1,5,10,25,50,100,150 J] n=3 的 BEC 激发 态 波 函数 ， 除 有 一 
个 节点 外 ， 同 基态 波 函 数 类 似 ， 其 分 布 随 Co 的 增 大 逐渐 变 
宽 . 图 4-6 给 出 了 Co =100，n =3,5,7,9,11,13 的 BEC 激发 态 波 
函数 ， 都 具有 超 高 斯 分 布 特征 ， 节 点 数 为 ”+ 图 4-7 给 
Hi fon -35/191113 的 激发 态 能 量 D, 随 Co 的 变化 ,各 能 级 的 
本 征 能 量 随 Co 增 大 而 增 大 ， 但 在 相同 Co 条 件 下 ， 不 同 激发 
态 能 级 与 前 一 能 级 之 差 是 不 同 的 ， 随 能 级 升 高 增 量变 小 ,如 
图 4-8 所 示 ， 它 给 出 了 n=3,5,7,9,11,13 的 激发 态 能 量 B, 与 
前 一 能 级 p,_, 之 差 随 Co 的 变化 关系 .可 近似 地 表示 为 
^B,-B,-B,o-2-alm*, a t CLE X. 4 0, =150 时 ， 
a = 1.19. 
为 了 达到 给 定 的 求解 精度 eps=104， 对 o (0) 精度 调整 
到 10”~ 107 就 能 满足 波 函 数 归 -化 的 精度 要 求 (4-15) 式 . 
而 要 满足 (4-14) 式 ， 对 B, 调整 精度 要 更 高 ， 且 随 所 选 x — x 的 
边界 值 的 增 大 而 变 高 ， 因 为 Cv 的 增 大 和 激发 态 的 升 高 都 会 使 波 函 
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x 


E 4-5 非 线 必 常数 Co 分 别 为 0.145.10.25.50.100,150 的 
n-3Íí'] BEC 激发 念 波 函 数 ( 按 宽度 增加 的 顺序 ) 


图 4-6” 非 线性 常数 Co = 100. n =3,5,7,9,11,13 4 BEC 
激发 念 波 函 数 KY (0) 值 由 小 到 大 的 顺序 ) 
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0 20 40 60 80 100 120 140 
C, 


0 


图 4-7 n=3,5,7,9,11,13 的 激发 态 能 量 B, 随 非 线性 
常数 Co 的 变化 ( 按 自 下 而 上 的 顺序 ) 


0 20 40 60 580 100 120 140 


0 


图 4-8 n=3,5,7,9,1113 的 激发 仿 能 量 B, 与 前 - -能 级 5，， 
之 差 随 Co 的 变化 ( 按 自 下 而 上 的 顺序 ) 
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数 分 布 逐渐 变 宽 . 例如 基态 的 实际 计算 中 随 C6 的 增 大 选取 边界 值 
x=8~10 ， 而 对 n=13 的 激发 态 随 C6 的 增 大 选取 边界 值 
X=11~12. EPRICE MAR IRL Bp, 和 边 值 B'(0) 的 依赖 性 
随 C 的 增 大 变 得 越 来 越 敏感 ， 亦 既 Cu 较 小 例如 C。< 10 ， 较 容易 
调整 B, 使 积分 收敛 满足 (4-14) 和 (4-15) 式 ， 但 当 C > 50 时， 对 
B, 的 调整 就 须 更 精细 ， 才 能 达到 收敛 . 


4.5 SAN 极限 下 基态 能 量 的 比较 

由 不 同 Cu 的 基态 波 函 数 可 以 看 出 ， 随 着 凝聚 原子 数 N 
的 增 大 ， 中 心 密度 变 得 越 来 越 平 ， 波 函数 分 布 宽度 增加 Dd 
此 , NLSE (4-4) 中 的 动能 项 的 贡献 仅 在 边界 附近 较 明 显 , B 
相对 于 相互 作用 项 随 N 的 增 大 变 得 越 来 越 不 重要 . 为 研究 
大 NN 极限 的 情况 ， 我 们 略 去 NLSE (4-4) 中 的 动能 项 ， 有 


gnoir'we) + NU, JPG YE) = uiv) (4-19) 

1/2 u- 1 mor? 
X Pren-|— o. CAR Po -——— — , t 
E | ( ) NU, f 


这 个 区 域 4> Val). MER ror, PHP -0. 这 就 
是 著名 的 托马斯 - 费 米 近似 . 波 函 数 的 归 一 化 条 件 可 用 来 确 
定 化 学 势 与 粒子 数 之 间 的 关系 


2/5 2/8 
15NU, Y. eap I5Na | 
—— (= ma;) = ho, et 
81 2 8(h/2mco,) $ 


«n 


(4-20) 
当 W 很 大 时 ， 我 们 将 前 面 的 数值 解 拟 合 结果 (4-17) 式 中 的 1.5572 略 
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差别 . 这 是 由 于 在 所 计算 的 N 值 的 范围 内 , 动能 项 虽 随 N 的 增 大 
而 减 小 但 还 不 能 忽略 ， 如 图 4-9 中 实 线 所 示 . 
为 讨论 化 学 势 与 能 量 的 关系 ， 我 们 将 (4-1) 式 写成 如 下 
形式 


E = En + Epo + Eim (4-21) 


三 部 分 能 量 分 别 随 凝 聚 原子 数 N 的 变化 关系 如 图 4-9 所 示 . 
直接 积分 (4-4) 式 ， 可 求 得 


H = (Ekin + Epo + 2E in )/N (4-22) 


8.0x10 L6x10* 24x10* 3.2x10° 
N 


图 4-9 凝聚 体 各 部 分 能 量 随 NN 的 安 化 关系 
实 线 为 动能 Ek、 虚线 为 势能 Epo. SER HH HEM fie E 


实际 BEC 实验 中 的 检测 是 通过 突然 关闭 势 阱 ， 让 发 生 
凝聚 的 原子 云 自由 扩散 ， 用 CCD 对 原子 气体 的 探测 光 成 
象 ， 图 象 数值 化 处 理 可 以 得 到 原子 云 速 度 分 布 对 原子 速 
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度 分 布 的 积分 可 得 到 原子 的 动能 ， 这 个 能 量 被 称 为 释放 能 ， 

它 是 在 自由 扩散 刚 开始 时 原子 的 动能 和 相互 作用 能 之 和 1 
Era = Ey + Eg (4-23) 


因此 ， 当 N 很 大 时 ， 释 放 能 主要 来 自 相互 作 用 能 . 此 外 ， 
根据 wp= BE/16N ， 由 (4-17) 可 求 得 


E 15 十 
N 7 (ul ho? 21,55? o (424) 
1.57/2 
(ul ho, 5? 一 Lsa no) 
以 及 巾 式 (4-20) 可 求 得 
eas (4-25) 


fea E 5i E53 u zr RITE E BER 4-10 中 ， 对 于 所 计算 


ee Eee 
00 gox10' 16x10" 2.4x10° 32x10* 


o 


图 4-10 能量 Elu)/ Nu FA u kA 
实 线 为 数值 解 结果 . 虚线 为 (4-24)} 式 的 结果 . 点 线 为 为 (4-25) 式 的 结果 
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的 大 N, (4-24) 式 的 结果 已 接近 托马斯 - 费 米 近似 (4-25) 
的 结果 ， 即 略 去 动能 项 的 贡献 . 这 时 (4-23) 式 的 释放 能 7 
E. 2 


时 = 二 4-26 
Nc (4-26) 


4.6 小 £& 


Æ E iN T AFR NEE AS BR IS BB BBR A PE 
即 对 非 线性 微分 方程 进行 数值 积分 时 ， 调 整 波 函 数 的 边界 
E (0) MATE 2, E DG) Ex 一 o 时 收敛 于 零 和 归 一 化 
条 件 都 得 到 满足 ， 并 求解 了 简 谐 势 阱 中 中 性 原子 的 非 线 性 
定 态 薛 定 请 方程 基态 和 激发 态 解 . 计算 了 基态 和 激发 态 的 
能 量 本 征 值 2, ， 边 值 必 (0) 及 空间 分 布 . 讨论 了 参数 @'(0) ， 
的 确定 与 积分 收敛 及 波 函 数 归 一 化 的 关系 . 
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第 五 章 吸引 相互 作用 对 BEC 的 影响 


5.13l & 


势 阱 中 ”Li 原子 气体 的 玻 色 - 爱 因 斯 坦 凝 聚 实验 523 为 这 种 
宏观 量子 现象 的 研究 提供 了 新 的 素材 . 锂 原子 的 凝聚 机 制 倍 受 
重视 ， 因 为 锂 原子 的 散射 长 度 忆 是 负 的 ， 这 就 意味 着 原子 间 有 
吸引 相互 作用 . 一 般 认 为 锂 原子 均匀 体系 不 能 形成 稳定 的 BEC 
态 ""”1， 这 是 因为 原子 间 的 吸引 相互 作用 将 引起 凝聚 体 的 崩塌 . 
然而 ， 具 有 吸引 相互 作用 的 中 性 原子 体系 在 空间 束缚 条 件 下 可 
以 形成 亚 稳 的 BEC 且 凝 聚 体 的 原子 数 No 小 于 某 个 临界 值 N,, 文 
献 [26,28,29] 预 期 在 充分 少 的 原子 数 或 非常 弱 的 原子 间 相互 作用 
条 件 下 凝聚 体 是 稳定 的 ， 当 这 些 条 件 不 满足 时 的 凝聚 体 的 衰变 
和 前 塌 也 有 几 位 作者 讨论 Ps3"”， Kagan PI 从 稍 不 同 的 论点 出 
发 指出 : 有 只 有 很 大 负 散 射 长 度 的 俘获 气体 团 答 的 形成 是 不 稳定 
的 ， 原 子 间 的 吸引 使 其 存在 更 稠密 的 N 原子 态 . Stoof "3 讨论 了 
这 种 吸引 型 BEC 的 热 诱导 崩塌 ， 文 献 [28,30] 估 计 了 宏观 量子 隧 
道 (MQT ) 的 速率 . 由 MOT 速率 可 进一步 近似 确定 吸引 型 BEC 
的 寿命 . 然而 ， 随 着 实验 精度 的 提高 ， 可 靠 的 理论 计算 是 必要 
的 . 本 文中 我 们 给 出 有 吸引 相互 作用 的 中 性 原子 体系 的 能 级 和 
态 函 数 在 一 维 箱 势 中 的 解析 解 中 ， 分 析 原 子 间 的 吸引 作用 对 能 
级 和 波 函 数 的 影响 ， 并 与 排斥 相互 作用 情形 中 性 原子 体系 的 解 
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进行 比较 . 通过 直接 数值 求解 球 对 称 简 谐 势 阱 中 的 非 线性 定 态 
OE PST BR, 求 得 系统 的 能 量 本 征 值 和 波 函 数 随 凝 聚 原子 数 变 
化 关系 . 类 似 于 排斥 势 情况 ， 求 解 二 维 简 谐 势 阱 中 有 吸引 相互 
作用 体系 的 NLSE 的 严格 解 踢 ， 计 算 密 度 涨 落 对 BEC 化 学 势 的 
影响 . 给 出 了 一 个 严格 可 解 的 宏观 量子 隧道 模型 "4 计算 了 从 亚 
稳 凝 聚 态 到 塌 缩 态 的 宏观 量子 隧道 速率 ， 研 究 吸引 型 BEC 的 稳 


5.2 负 散射 长 度 的 计算 与 测量 
按 分 波 法 ， 弹 性 散射 的 相 移 的 计算 公式 为 131 


tan d, lp july Rr dr (5-1) 
式 中 6 为 / 波 相 移 ，j,(r) 为 球 贝 塞 耳 函 数 ，V(n) 为 中 心 势 ，R(r) 
为 径 向 波 函 数 . 将 渐 近 的 径 向 波 函数 Ry(r)~sin(hr - n3, fR 
入 (5-1) 式 即 是 相 移 的 积分 方程 . 一 旦 知道 了 原子 靠近 时 的 相互 
作用 势 V(r) 的 形式 ， 就 可 以 求 得 相 移 5 ， 从 而 由 (3.1) 求 得 散射 
Kia. 
文献 [96] 使 用 半 经 典 近 似 对 原子 间 相互 作用 势能 有 渐 近 形式 
V(r)~— = 的 原子 -原子 弹性 散射 进行 了 计算 ， 求 得 的 散射 长 度 为 


se [B te 
n-2} h(n-2) -l 
T =) (5-2) 
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posf JMV, METETA. 

对 于 处 于 基态 的 Li 原 了 了， 相应 于 Lih x Ez Rr a^; 态 的 
近似 势能 函数 V(r) RV, (7) 分 别 为 5 


6 
MGE p — DE 10" 10012885559 e7129 
r r 


(5-3) 

在 零 能 极限 下 求 得 8] 相 应 于 势能 函数 大 (rm) Cx! £z AS) 和 V(r) 

Ca? E; AS) 的 散射 长 度 分 别 为 36 ao 和 -17co ， ao 是 玻 尔 半径 . 

为 从 实验 上 确定 "Li 原子 散射 长 度 a， 文 献 [98] 通 过 二 光子 

拉 曼 过 程 探测 基态 分 子 能 级 的 方法 对 三 重 态 ?Li airy 势 的 

最 低 束缚 态 的 束缚 能 进行 了 测量 ， 测 得 离 解 能 De = 

333.73+0.02cm-! ， 从 而 得 出 处 于 所 = 2, m, = 2 三重 态 s 波 散 射 
长 度 的 值 为 (-27.3 土 0.8)ao . 


5.3 一 维 箱 势 中 有 了 吸引 作用 体系 的 能 级 与 波 函 数 


参照 第 三 章 (3.24) 式 ， 有 弱 相 互 作用 中 性 原子 体系 的 波 函 
数 所 满足 的 定 态 非 线性 薛 定 刘 方 程 为 : 


ctv +V 4 gl?) Y= EY (5-4) 
x 
式 中 g = Nu ，N 为 凝聚 原子 数 ， x = 4 2 zm. 


为 负 值 . 由 式 (5-4) 可 求 出 中 性 原子 系统 的 本 征 态 多 | 及 相应 的 本 
TE E,. 为 简单 起 匈 ， 我 们 先 讨论 一 维 箱 势 中 的 原子 体系 ， 这 
时 方程 (5-4) 可 写 为 
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nh? od? 2 
-一 +g¥ Y =E¥ , 0<x< 5-5 
2m dx? sl | ; "end (5) 


4v-NU,, HES-SyEX 


2 


CL qd -adyeadwiy, =0 (5-6) 
dx 
stp 2 -242 -2mE IW? , kom Nas kp Ña, N[v2ax=1, 
a i 
on k? 
N 为 归 一 化 常数 . 换 变数 并 引入 新 参量 a, 
sch 
E= kn -ka x，(5-6) 式 变 为 
d? 
(> +(1-27) +2272, =0 (5-7) 
dé? 
进一步 得 到 
ELA i =1- (1-7)? -A 
ar)" wy -av, (5-8) 
(5-5) 式 的 积分 为 


Yy 


dv, 


n 


= \-— (5-9) 
OY- + AY) 


因此 , EBAY, 可 用 第 一 类 Jacobi 椭圆 函数 来 表示 区 ,= sm&,i4). 


142? 


通过 虚 模 数 变换 sn(E, i) = TIU m 波 函 数 


Y, 可 表示 为 
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Y, =sn(EiA)=A 1 sd(--,A) = 


sn(£ A) (5-10) 


+A -sE T 
EX 


EL 
函数 sn(&',4?7) 有 下 面 的 性 质 


AIL (5-11) 


o yU- 1-27 ¥,") 


mm JRR, sE, À) 的 展开 式 为 


式 中 kd n- jT K(R)= K-27) ， 


a 
~ Ír & d " . Qj*Dne 
sn(ë' A>) = = ih (5-12) 
d AK Lp -qU" 2K 


Hpg, A nk, al m= 4—mga! mh Tk, 的 定义 ， 得 到 


UR Na=(k a/n) Na = n; Na (5-13) 
Th 


由 归 一 化 条 件 ， No 可 表示 为 


ke -kla 
z (5-14) 


根据 x=0,a, 里 =0 的 边界 条 件 , 即 sm( -k aA)= 0, 
ec 
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fit PARTE 


ATE -kè a= A 2nK,n-123,— (5-15) 
考虑 到 波 函 数 的 周期 性 ， 则 (5-14) 式 变 为 
AAxXKO 0 1 AK 


- (5-16) 

rz ru K y P (EdE 1-22 ur (d£ 
利用 量子 化 条 件 (5-153) 和 4 与 元 的 定义 式 ， 有 k=2nKX1a 和 
kè IK) -kD)- AI - 22). 代入 (5-13) 式 得 
_ 2K4 


Nan 


3 |a 


(5-17) 


当 nn 给 定 后 ， 由 (5-16)，(5-17) 式 可 求 出 本 征 值 亏 Je KÕ 2), 
便 可 计算 出 心 及 心 以 及 能 量 本 征 值 


Ros 3、 AUR (2K 
E, =— (k? -2k2 
1 ag 24 Pre (24) 


图 5-1 给 出 当 , 给 定 后 ， 本 征 值 4, 随 的 变化 曲线 . 由 这 


HER AT EH SAEI 2.25.2... 图 5-2 Hite, BE n AEE. 由 图 


看 出 ， —. - 维 体系 的 能 级 变化 情况 相反 , 吸 
引 作用 体系 的 能 级 相对 无 相互 作用 体系 降低 了 ， 基 态 与 低 激发 
态 能 级 对 原子 间 吸 引 相互 作用 强度 反映 灵敏 ， 且 随 强 度 增 大 而 
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变 小 以 致 变 为 负 值 ， 这 是 原子 间 吸 引 相互 作用 的 反映 . 


o 
5 
n 
s 
- 
Ss 
a 
è 
w 
3 


图 5-1 REA A, Bln ELIA 
Hl, 2, 3, 4 分 别 对 应 Ny =0.1, 0.25, 0.56, 1.0 


0 
ey 
2 
4 
| SS ae 
0 2 4 6 8 


图 5-2 KER e, Min 4E XE 
HERI, 2, 3, 4 HM n, =0.1, 0.25, 0.56, 1.0 


4 g>0, MA+Ont, Kon/2, V,, E, 趋 于 熟悉 的 无 相 
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互 作用 情况 下 的 波 函 数 寡 = V2/asin(ma/a 和 能 量 本 征 值 


E, =m nh? /2ma?. 


图 5-3 基态 波 函 数 随 mg 的 变化 关系 
RRI 2 3 4 分 别 对 应 于 ,二 V-4NE/za -= 0.75, 1.0, 125, 1,5 (吸引 势 
SUI, 2,3, 4 PRR n — V ANa 12a = 0.75, 1.5, 2.5, 5.0( 排 斥 势 ) 
虚线 对 应 于 me=0 (无 相 所 作用 ) 


在 图 5-3 中 给 出 了 体系 在 不 同心 (8 <0 ) 情 况 下 的 基态 波 函 
数 ， 为 比较 ， 也 将 第 三 章 g > 0 情况 下 的 基态 波 函 数 (3.31) 式 画 在 
同一 图 中 . 由 图 看 出 ， 随 相互 作用 强度 的 增 大 ， 吸 引 作用 体系 的 原 
子 空间 分 布 趋 于 集中 ， 而 排斥 作用 体系 的 原子 分 布 趋 于 分 散 . 


5.4 三 维 箱 势 中 有 吸引 作用 体系 的 基态 解 


根据 方程 (5-4)， 半 径 为 的 球 方 势 阱 中 s WARO O) 满足 
的 方程 可 写 为 : 
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-B 1 0 yehe ve) EvO) (5-19) 
2m r? dr? 


采用 无 量 纲 的 长 度 、 能 量 单位 ， 仿 : 
r=ryx, B= 2ma! E/h? 
Co =2NG/ ro (5-20) 
P(r) = (4125) ^ d(x) 
得 到 : 
2 
í 4 2-£+9-€,0%9 jen -e (5-21) 


dx) xdx 


波 函 数 的 规 一 化 条 件 为 : 

an [t ny? dr= {Por dx=1 (5-22) 
HAR D(x) 满足 的 边界 条 件 是 : 

DO =0, DBO), -(0,0' (x)... =£, (6-107) (5-23) 


应 用 4.2 节 中 的 数值 计算 方法 ， 方程 (5-19) 满 足 条 件 方程 (5-22)， 
(5-23) 的 数值 解 画 在 图 5-4 ~ 图 5-6 中 . 

图 5-4 给 出 非 线性 系数 c。 随 基态 能 量 8 的 变化 关系 ， 当 
B=-0.225 时 ， 非 线性 系数 有 一 极 小 值 Cu,，= -0.65276， 由 此 和 
方程 (5-20) 可 以 推出 临界 原子 数 Ne, 它 是 吸引 型 BEC 能 容纳 的 
最 大 原子 数 . 为 与 实验 比较 ， 取 实验 数据 中 17F145nm ， 
n=54m， 可 得 出 Ne =roComin /#5=1080. 这 与 简 谐 势 阱 中 临 
界 的 凝聚 原子 数 1300 059 是 同一 量 级 . 

对 每 一 个 Cu> Cuw。 都 对 应 具有 不 同 本 征 值 的 两 个 本 征 态 ， 如 
图 5-6 中 的 波 函 数 曲 线 1 与 1 分别 对 应 B=5.14 和 -10.27 ， 但 它们 
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对 应 同一 个 韭 线性 系数 C。= -0.45 ， 或 说 它们 有 相同 的 原子 数 . 这 
表明 给 定 原子 数 的 吸引 相互 作用 体系 存在 双 稳 态 . 对 于 能 量 本 征 值 
p> -0.225 的 态 , 凝聚 体 原子 数 随 本 征 能 量 8 的 增加 而 减少 . 而 对 
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图 5-5. AES FON 与 能 量 木 征 值 8 的 关系 
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5-6 不同 非 线 性 系数 Co PAE BATE 8 的 基态 波 函 数 


T B«-0225 的 态 ,凝聚 体 原子 数 随 8 的 减少 而 减少 ， 但 原子 
的 空间 分 布 越 来 越 集中 ,如 图 5-6 所 示 , 波 函数 曲线 3', 2', 1 分 
别 对 应 B = -2.5788, - 5.2523, -10.2274. 我 们 的 这 个 结果 类 似 于 
简 谐 势 阱 中 吸引 相互 作用 体系 非 线性 薛 定 刘 解 的 稳 态 和 不 稳 
态 分 支 U. 因此 球形 箱 势 中 吸引 型 BEC 在 压缩 的 空间 分 
布 、 临 界 原子 数 和 凝聚 体 的 亚 稳 特性 等 方面 类 似 于 简 谐 势 中 
的 同类 体系 . 


5.5 有 了 吸引 相互 作用 的 中 性 原子 的 BEC 计算 


由 5.3 节 的 结果 可 以 看 出 ， 在 一 维 情形 吸引 势 的 BEC 并 不 
受 限 制 ， 而 5.4 节 三 维 箱 势 和 后 面 简 谐 势 阱 中 的 数值 结果 都 表 
明 , 吸 引 作 用 使 得 三 维 非 线性 薛 定 刘 方 程 在 当 非 线性 系数 或 BEC 
原子 数 要 充分 小 时 ， 方 程 才 有 解 . 为 了 定性 地 研究 原子 间 吸引 
相互 作用 对 BEC 的 影响 ， 仍 近似 地 用 :三 个 - 维 解构 成 三 维 情形 


PE 
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的 结果 ,但 应 当 考 虑 三 维 情 形 中 原子 数 或 非 线 性 系数 的 限制 ， 为 
此 我 们 取 基 态 能 量 为 零 ， 由 求 得 的 能 级 公式 (5-18) 可 以 看 出 ， 这 
24, =1， 可 求 得 这 时 相应 于 mw =0.56. 因此 ，n 在 0~0.56 
的 变化 范围 内 研究 其 对 BEC 影响 . 在 作 凝 聚 计算 时 ， 参 照 第 三 
章 排斥 势 情况 的 方法 ,基态 和 激发 态 的 原子 数 分 别 为 


NET «r3, -38) 


N,-0Q > 一 (5-24) 


e -a( pis, qu +r2er -348) 


式 中 o= z=ef HSE) 化 学 势 y 238, Pe; pq, 


9,r=1,2,3,…,) 是 体系 的 第 j EASE CIS EMA, dust Ys p 
不 包括 p=g=+=1 的 项 . 实际 计算 时 ， 我 们 将 oo 取 为 1. 为 求 
得 系统 发 生 BEC 的 临界 温度 7 以 及 g 对 7 的 影响 ， 我 们 计算 
了 当 = 一 1 时 wz 随 c ”的 变化 关系 .计算 结果 表明 : 在 同一 温 
度 下 ， 随 吸引 势 的 增强 ， 激 发 态 能 容纳 的 原子 数 N; 略 有 下 降 ， 
这 与 排斥 相互 作用 情况 类 似 ， 同 样 因原 子 间 的 相互 作用 ， 使 
激发 态 (n >1) 能 量 与 基态 (n=1) 能 量 之 差 增 大 了 ， 这 从 图 5-2 
"Poe, 随 n 的 变化 关系 可 以 看 出 来 . 我 们 取 定 总 原子 数 为 
N=1250 ， 改 变 o 使 按 (5-24) 式 求 得 的 激发 态 原子 数 N, 为 
1250， 这 时 对 应 的 c” 就 是 凝聚 的 临界 温度 ， 用 cz! 表 示 . R 
$F n,770, 0.45, 0.56 对 应 的 cz 分别 为 cz! =96.9, 97, 97.15, Bi 
吸引 势 的 增强 略 有 增 大 . 即 在 所 考虑 的 范围 内 g 对 临界 温度 
的 影响 很 小 . 
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5.6 简 谐 势 阱 中 吸引 作用 体系 的 基态 

在 充分 低 的 温 ei ee 的 条 件 下 ， 弱 相 
互 作用 臻 色 子 的 诗 聚 体 可 通过 单一 波 函 数 Y(r) 来 描写 ， 原 子 间 
的 相互 作用 由 s- 波 散射 长 度 5 来 表示 ， RICO 可 通过 求 下 
面 Gross-Pitaevskii 能 量 泛 函 的 极 小 而 得 到 P9]: 


EW) = far’ (r)x Jy + Emo? rt 2 wy Wir) 


(5-25) 
AYE ae AS fe XC FRAY ES AA AE EE PE BE SS 
方程 : 


My v lmol rrr) 
2m (5-26) 
Sa e yr) = po) 
AH a 是 化 学 势 、 ¥(r) 的 归 一 化 为 : 
Jaro = ws (5-27) 
其 中 No 是 凝聚 体 中 的 原子 数 . 考虑 到 BEC 基态 波 函数 是 球 对 
称 的 ,， 作 变 换 ， 使 长 度 和 波 函 数 无 量 纲 化 : 
z(h/2mo, ) ^ x2 ax, B-ulho, 
V2 (5-28) 
Co =2Noala, V(r)- (v, /4na*) Ox) 


这 样 ， 方程 (5-26) 有 下 面 的 形式 
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[eit 2-3-c he joey =0 (5-29) 
= x dx 4 
外 (5-27) 和 (5-28) 式 ， 波 函数 @(x) 满足 的 归 一 化 条 件 为 : 

[e^ cox? dx=l (5-30) 

0 
I eS IC GO) 满足 的 边界 条 件 是 DG] | = 0,0). -0 和 
D'(x)| =0 P(x), = 90), (£=10%). 

再 次 应 用 4.2 节 的 数值 计算 方法 ， 我 们 求 得 了 基态 波 函 数 

D(x) 和 能 量 本 征 值 = ul ho, 的 数值 解 ， 画 在 图 5-7~ 图 5-9 rp. 


图 5-7 Ala] Cy A B MEAIB 


图 5-9 给 出 的 基态 能 量 6 随 非 线性 系数 Ce 的 变化 关系 有 两 
个 特点 , 第 -， 对 十 所 有 的 非 线 性 系数 C。( 或 说 凝聚 原子 数 )， 
基态 能 量 P<1.5. 这 表明 原子 间 的 吸引 相互 作用 使 体系 的 能 量 
FRET. 第 二 ， 在 6=0365 处 非 线性 系数 CW 有 一 极 小 值 
Comin =-1.62625. 由 此 和 方程 5-28) 可 以 推出 临界 原子 数 w.， 它 
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是 吸引 型 BEC 能 包含 的 最 大 原子 数 ，Ne = -1.626 x a /12a. 


1.00 


5-8” 波 函数 的 宽度 9 随 木 征 值 8 MAE 


图 5-9 非 线 必 系数 Co 与 能 量 本 征 值 之 间 的 关系 
内 插图 为 能 量 本 征 值 P 随 凝聚 原子 数 与 临界 原子 数 之 比 的 变化 关系 
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当 阱 的 宽度 增加 时 ， Ne 将 随 K 线 性 增加 ， 这 -结果 似乎 与 

均匀 气体 中 不 会 发 生 吸 引 的 BEC HFE, 其 实 , 当 a o oo 时 , BE 

REEN /a FE. N 和 凝聚 密度 之 间 的 权衡 正 是 设计 实 
验 时 所 要 考虑 的 . 

根据 Co 的 极 小 值 和 实验 数据 局 |aF145 nm，w= 


| h 
(osovo-) ? 2908s, mj, 2116x107 g, a= I 72236 um, 


可 得 出 Ne =1254. 这 与 实验 中 最 大 的 凝聚 体 原子 数 限 制 在 650 
和 1300 之 间 Ol 是 相符 的 . 对 每 一 个 Co > Cou, 都 对 应 具有 不 
同 本 征 值 的 两 个 本 征 态 ， 如 图 5-7 中 Co = -1.033 相 应 波 函数 曲 
线 1 与 6 分 别 对 应 =1115 和 -1.75 ， 这 就 是 我 们 所 谓 的 双 稳 
S. 对 于 能 量 本 征 值 6 > 0.364 的 态 , 凝聚 体 原子 数 随 本 征 能 量 B 
的 增加 而 减少 ， 而 对 于 8 < 0.364 的 态 , 凝聚 体 原子 数 随 5 的 减 
少 而 减少 ， 但 原子 的 空间 分 布 越 来 越 密集 ， 如 图 5-7 所 示 ， 压 
缩 的 波 函 数 曲线 4,5,6 分 别 对 应 5 = 0, -0.75, -175. 这 种 压缩 
的 波 函 数 可 能 与 文献 [26] 提 到 的 稠密 态 有 关 . 本 征 态 解 表明 : 对 
TN « Ne 具有 吸引 相互 作用 的 原子 体系 的 BEC 可 以 形成 . 

另外 ， 我 们 可 以 求 得 处 于 临界 原子 数 状态 时 平均 每 个 原子 
的 相互 作用 能 为 

2xh?à "EM 

m 


Um T Sd? rH = xe Jaxx? ow- 
cm s 


Ne2nh’a 1 1.626x 
M ÉL fax Ox) t= DEUM [axe | D(x) |= 
m 4na^ 
0.804 ic) 
(5-31) 


由 此 可 见 ， 负 相互 作用 能 的 数值 小 于 简 谐 势 阱 的 能 级 间隔 
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hw ， 以 致 于 相互 作用 表现 为 对 理想 气体 的 微 扰 ， 仍 为 气体 
状态 ; 这 也 可 以 从 物理 上 理解 为 什么 吸引 相互 作用 体系 有 一 
个 临界 原子 数 ， 或 说 方程 (5-29) 仅 当 N < N. 时 才能 找到 数 
值 解 . 

对 于 能 量 本 征 值 8 > 0.364 的 态 , 凝聚 体 原子 数 随 本 征 能 量 
B 的 增加 而 减少 .而 对 于 < 0.364 的 态 , 凝聚 体 原子 数 随 5 的 减 
少 而 减少 ， 但 原子 的 空间 分 布 越 来 越 密集 ， 如 图 5-7 所 示 ， 压 
缩 的 波 函 数 曲线 4,5,6 分 别 对 应 6 = 0, -0.75, -1.75. 从 原子 的 
空间 分 布 与 本 征 能 量 的 关系 看 ， 原 子 的 空间 分 布 随 基态 能 量 的 
降低 单调 地 被 压缩 ， 如 图 5-8 P q 随 8 变化 的 曲线 所 示 ，4 是 原 
子 分 布 的 宽度 ， 由 O(g) 700) 定义 . FI ®, O 表示 宽度 为 
9 的 波 函 数 . 
5.7 吸引 性 BEC 的 稳定 性 分 析 

为 研究 上 面 所 求 的 两 类 定 态 的 性 质 ， 即 凝聚 体 的 稳定 性 ， 
由 于 凝聚 体 的 不 稳定 是 密度 涨 落 的 结果 ， 涨 落 将 使 系统 进入 更 
低 的 能 量 状态 ， 因 为 原子 间 有 有 效 的 吸引 作用 . 为 此 ， 我 们 用 
基态 波 函 数 的 宽度 来 表示 凝聚 密度 ， 由 (5-25) 和 (5-28) 式 可 求 得 
系统 基态 能 量 与 凝聚 密度 的 关系 


5 

ENGE doS ,., 2Noa 

E(D,)= d pe ne *x $,(x)*2 Nog 
2 5 


4 
D, (x - 
2ma a a ) 


2 


h^ No 
2ma? 


E(q) = Ny V (q) 
(5-32) 
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所 求 得 的 ECg) 随 7 的 变化 关系 如 图 5-10 所 示 ， 计 算 中 使 用 的 是 
不 同 凝 聚 数 No 本 征 解 ， 即 图 中 的 密度 变化 是 由 于 凝聚 数 No 的 


变化 产生 的 . 

8.0 
75 - m 
70 

S5 

R 60 
5.5 06 09 12 15 

q 


图 5-10 E(9) 随 9 的 变化 ， 内 插图 为 局 部 放大 . 


为 研究 给 定 凝 聚 数 N。， 因 原子 的 空间 分 布 的 涨 落 对 体系 稳 
定性 的 影响 ， 需 要 对 给 定 凝 聚 原子 数 N 的 体系 计算 泛 函 
E(BD,(x)) ， 波 函数 的 宽度 gq 作为 变 分 参数 . 给 定 No 可 根据 (5-32) 
式 求 得 E(q) 随 g 的 变化 ; 对 于 No/N。 = 0.992, 0.993,…,1.0，E(g) 
随 4 的 变化 如 图 5-11 中 的 实 线 所 示 . 虚线 为 各 E 曲线 的 极 
值 点 的 轨 线 ， 它 恰好 相应 十 方程 (5-29) Cy 22N.a/a 的 严 
格 解 B,(x). 由 图 可 看 出 , 8>0.365 的 态 在 E(q) 的 极 小 点 , 称 
它们 为 亚 稳 的 BEC 态 ， 并 用 Bo w CD 表示; ifi] «0.365. 的 定 


态 在 E(g) 的 极 大 点 ， 称 它们 为 不 稳定 的 稠密 态 24， 表 示 为 
Di wo (x). 
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图 5-11 给 定 Xo 时 (9) 随 9 的 变化 关系 实 线 : 
E(q) 作为 9 的 函数 ， 相 上 应 于 No /N。 =0.992, 0.993,…,1.0 
(由 上 到 下 的 顺 订 )， 虚 线 相应 于 这 些 曲 线 的 极点 轨 线 


5. 8 宏观 量子 隧道 速率 的 计算 
根据 宏观 量子 隧道 (MQT ) 速率 的 WKB 近似 公式 吧 和 方 
程 (5-32)， 我 们 得 到 : 
Ty = Aexp[-2N, In fager vro - ve y] ^ j = 
Aol 1225N, [dal E(q)~ Eqo)] ^ | 
式 中 有 效 质量 m =3m/2, qo 是 亚 稳 极 小 的 位 置 ， q, 满足 
E(q,)  E(qo) . 利用 方程 (5-32) 的 E(q) 可 计算 从 亚 稳 态 到 (x) 到 


do 
塌 缩 态 D, (x) POE GE To. 图 5-12 中 曲线 1 XR D/A B 
No/N。 变化 的 规律 . A ECE. EA ROAR ADO G B S B s o 
To ec exp[-0.57(N, - No)] FUR yA ik oR PR CPP SE FERE 


(5-33) 
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HR Dy x expl-(1- No/N,)] MER -图 中 . 很 明显 ， 严 格 解 
曲线 1 比 根据 文献 近似 解 计算 的 结果 曲线 2、3 BÉ No /Ne 的 减 
小 更 快 地 趋 于 零 . 这 意味 着 在 临界 点 附近 亚 稳 凝聚 可 以 很 容易 
地 通过 宏观 量子 隧道 而 进入 塌 缩 态 ( 如 图 5-11 中 标 出 的 9 点 相 
应 的 态 ); 然而， 远离 临界 点 时 ， 比 Ne 少 几 十 个 原子 ， 这 种 隧 


|: exact solution 
2: exp[-(1-N JN. ^] 


3: exp[-0.57(N -N)] 


3 


{554 0.995 0.996 0.997 0.998 0999 1.000 
NJN, 
图 5-12 隧道 速率 1,1 4 BÉ Nuy/N。 的 变化 
曲线 12.3 分 别 相应 于 产 格 解 ，exp[-(1- Au/Ne 和] 和 exp[-057(N。- Noj] 
虐 线 为 亚 稳 凝聚 念 oo ws (x) AREAS Oy y, (x) 之 间 的 交 一 积分/ 


5. 9 跃迁 速率 的 计算 

我 们 根据 Kagan P9 关于 因 涨 落 从 亚 稳 态 到 稠密 态 跃 迁 的 激 
BAUR EGE EMIS MTS, 我们 可 以 应 用 我 们 求 得 
的 严格 的 本 征 解 Dowo (x), Dano (x) 代替 文献 [26] 的 高 斯 近似 解 
pos po KBD nox), Qu GO) 都 是 归 一 化 的 本 征 函 数 ， 即 单 
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粒子 波 函 数 ， No 粒子 波 函 数 应 是 单 粒子 波 函 数 之 积 ， 交 迭 积 分 
1 可 短 为 


IQ =([x° axe oo y, o9)" (5-34) 


(No) 的 计算 结果 如 图 5-12 中 的 虚线 所 示 ， 它 在 临界 点 附近 有 
较 大 的 值 ， 而 远离 临界 点 它 很 快 地 变 为 零 . 这 提示 出 在 临界 点 
附近 的 亚 稳 凝 聚 态 可 以 很 容易 通过 “个 集体 交换 作用 而 跃迁 到 
dus. xn UHR. 事实 上 在 临界 点 附近 的 两 种 态 它 们 几 
乎 有 完全 相同 的 化 学 势 和 空间 分 布 ， 但 它们 的 稳定 特性 是 完 
相反 的 . 然而 ， 在 远离 临界 点 时 ， 它 们 的 化 学 势 和 空间 分 布 相 
差 其 远 ， 这 种 跃迁 是 很 困难 的 . 


5.10 小 结 


我 们 求 得 了 一 维 箱 势 中 具有 吸引 相互 作用 的 中 性 原子 体系 
的 能 级 与 波 函 数 的 解析 解 ， 分 析 了 吸引 相互 作用 对 原子 体系 的 
能 级 和 波 函 数 的 影响 ， 并 与 排斥 相互 作用 的 中 性 原子 体系 的 结 
果 进 行 了 比较 ， 相 对 于 无 相互 作用 体系 的 能 级 ， 排 斥 体系 的 能 
级 升 高 了 ， 原 子 分 布 趋 于 分 散 ; 吸引 体系 的 能 级 降低 了 ， 原 子 
分 布 趋 于 集中 . 由 三 维 球 方 势 和 简 谐 势 阱 中 具有 吸引 相互 作用 
原子 体系 NLSE 的 本 征 解 的 数值 解 ， 得 到 类 似 的 特性 ， 即 原子 
的 空间 分 布 随 相 互 作用 的 增强 而 被 连续 压缩 ， 得 到 凝聚 原子 数 
随 能 量 本 征 值 变化 的 双 稳 态 曲线 ， 并 由 此 得 到 与 实验 报道 相符 
的 吸引 玻 色 - 爱 因 斯 坦 凝 聚 所 能 包含 的 最 大 原子 数 N。， 本 征 态 解 
表明 : 对 于 N < Ne 有 具有 吸引 相互 作用 的 原子 体系 的 BEC 可 以 
形成 . 对 每 -个 N < WN. 存 在 一 个 业 稳 凝聚 态 和 相应 同 -- 原 了 数 
的 不 稳 的 稠密 念 . 在 临界 点 附近 炎 稳 凝聚 可 以 很 容易 通过 宏观 
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当 稳 定 的 . 而 这 种 崩塌 的 动力 学 还 没有 详细 的 理论 描写 ， 因 为 
密度 会 变 得 人 得 足以 要 考虑 多 体 相互 作用 . 然而 ， 当 保持 低 的 
密度 时 ， 凝 聚 体 的 演化 可 以 由 一 个 准 粒 子 在 一 个 有 效 势 场 中 的 
运动 来 描写 . 对 于 7 > 0 的 捕 陷 气体 ， 崩 塌 之 后 仅 初始 不 在 凝聚 
体 中 的 那些 原子 被 保留 ， 所 以 这 时 的 气体 处 于 非 平衡 分 布 . 弹 
性 碰撞 重新 热 化 气体 ， 凝 聚 原子 数 继续 增加 . 如 果 气 体 是 被 有 
效 冷却 的 ， 那 么 凝聚 将 被 继续 填充 和 崩塌 ， 直 到 达到 7=0 或 
没有 原子 被 保留 . 这 就 是 文献 [100] 提 出 的 在 茹 发 冷却 过 程 中 凝 
聚 原子 数 的 生长 和 崩塌 的 演化 模型 . 但 这 个 模型 还 没有 被 证 实 . 
凝聚 原子 数 的 涨 落 已 经 被 观察 到 ， 但 还 不 清楚 它 是 内 在 的 还 是 
由 于 捕 陷 装 入 和 蒸发 冷却 的 变化 引起 的 . 如 果 这 些 变化 可 以 控 
制 ， No 可 以 准确 测量 ， 涨 落 的 测量 可 为 凝聚 生长 和 崩塌 的 动 
力学 提供 更 详细 的 信息 . 它 也 将 为 多 体 量子 理论 提供 更 重要 的 


. 93 。 


2000 年 上 海 大 学 博士 学 位 论文 


第 六 章 ”BEC 形成 的 动力 学 


6.1 引言 


当 温 度 降 得 很 低 时 ， 由 玻 色 原子 气体 凝 成 下 色 凝聚 体 的 过 
程 ， 可 理解 为 一 个 个 孩 色 原子 由 气态 跃迁 到 凝聚 态 . 设 凝聚 体 
中 包含 个 原子 ， 其 基态 为 |n > ， 当 一 个 气态 原子 增加 到 凝聚 
体 |n> 并 使 凝聚 体 推进 到 状态 |n+1> 的 几率 会 由 状态 jn> 中 所 
含 的 原子 数 n 而 得 到 增强 . 这 与 光学 激光 器 中 的 增益 机 制 类 似 
存在 于 某 光 场 模式 中 的 光子 数 将 会 激发 更 多 的 光子 向 同一 模 
式 辐射 ， 亦 即 受 激 辐射 对 于 玻 色 原子 情形 ,这 种 受 激 过 程 便 是 
导致 物质 波 的 放大 ， 即 原子 激光 (Atom laser) 获 得 增益 . 这 就 是 
1997-1998 年 MIT 年 所 进行 的 原子 激光 实验 B%， 与 此 相应 
Gardiner 发 表 了 一 系列 的 理论 研究 B14， 他 们 得 出 描述 凝聚 
体 密度 矩阵 p 的 主 方程 ， 并 进一步 简化 得 出 原子 激光 的 变 率 方 
程 . 解 变 率 方程 便 得 出 凝聚 体 原子 数 ” 随时 间 c 的 演化 曲线 
但 未 能 求 得 的 统计 起 伏 . 本 文 在 阐述 凝聚 体 主 方程 的 物理 意义 
的 基础 上 ， 用 解 微 激光 主 方程 的 母 函数 法 texiog 求 得 主 方程 的 解 
析 解 9， 进一步 计算 出 凝聚 体 原子 数 <n>， 增 率 na2 随时 间 
的 演化 ,以 及 凝聚 体 原子 数 的 统计 起 伏 < (Am > 等 .主要 结果 是 
增 率 在 发 展 过 程 中 有 一 个 平台 ， 而 统计 起 伏 为 超 Poisson 分 布 
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6. 2 BEC 的 动力 学 模型 


62.1 凝聚 带 与 非 凝 聚 带 

当 势 阱 中 有 强 的 凝聚 时 ， 阱 的 最 低能 级 将 演化 为 宏观 占 
据 数 ， 并 修改 阱 中 的 其 他 能 级 . 如 果 阱 是 紧 束缚 的 ， 较 低 的 
能 级 的 分 离 应 是 合理 的 . 但 紧 束缚 的 概念 不 是 绝对 的 ， 对 低 
能 级 为 紧 束缚 势 ， 而 对 较 高 能 级 却 是 十 分 宽松 的 . 其 本 质 区 
别 是 态 密度 ， 简 谐 势 阱 中 的 能 态 密度 随 能 量 的 二 次 方 而 增加 ， 
因此 对 于 某 个 适当 的 能 量 值 ， 它 已 有 非常 稠密 的 能 级 . 因此 
可 以 认为 凝聚 效应 对 系统 能 级 的 修正 仅 对 较 低能 级 才 是 明显 
的 . 将 粒子 状态 分 成 凝聚 带 Rc 和 非 凝聚 带 Rnc， 场 算 符 写 为 
V (r) = 660) + v wc (r). 两 个 带 的 划分 是 借助 某 个 能 量 Ex， 选 
择 这 个 能 量 使 非 凝 聚 带 中 所 有 粒子 的 能 量 都 足够 大 ， 以 致 于 
最 低能 级 上 的 凝聚 体 的 存在 对 它们 的 能 量 没 有 明显 的 影响 . 
因为 凝聚 原子 是 分 布 在 很 窄 的 空间 区 域 上 . 图 6-1 为 简 谐 势 阱 


E 


position 


图 6-1 凝聚 带 和 非 凝 聚 带 示 意图 
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中 两 带 的 示意 图 . 这 样 我 们 用 分 离 能 级 来 描写 E< En 的 凝聚 
带 Rc ， 显 然 它 不 仅 包括 凝聚 在 其 中 形成 的 基态 ， 也 包括 受 凝 
聚 体 存在 较 大 影响 的 那些 能 级 .而 非 凝 聚 带 RNc 中 的 原子 作 
为 热 原子 ， 代 表 作 为 热 库 的 大 多 数 原子 ， 它 为 凝聚 生长 提供 
了 粒子 源 ， 由 通常 的 平衡 公式 来 描写 . 
622 凝聚 体 粒子 数 守恒 的 Bogoliubov 近似 

凝聚 带 Rc 中 的 行为 用 全 量子 力学 处 理 ， 势 阱 的 能 级 因 
凝聚 体 的 存在 而 被 修改 . 在 任意 时 刻 Rc 中 有 一 个 给 定 的 原 
子 数 ， 在 这 种 情况 下 的 能 级 由 粒子 数 守 恒 的 Bogoliubov Jy 
法 描写 "I 因此 ，Rc 中 的 态 完全 由 Re 中 的 原子 数 和 Rc 内 
准 粒子 的 量子 态 来 描写 . 可 将 凝聚 带 的 场 算 符 写 为 


b bn 
$t) = Bl Eyt) y ie em) reds 


= VN 
能 量 为 sy 的 准 粒 子 由 消灭 算 符 b,, 描写 ,而 8+ 是 产生 算 符 ， 对 
于 任意 的 N， 它 使 系统 的 Re 由 w 原子 基态 进入 N +1 基态 . 
En C) 代表 凝聚 体 波 函数 ， 它 满足 定 态 的 Gross-Pitaevskii 方程 . 
Sin (7) 和 8gw(r) 是 能 量 为 sy 的 准 粒子 的 产生 和 潭 灭 的 振幅 在 粒 
子 数 守 恒 的 Bogoliubov 方法 中 ， 原 子 数 是 守恒 的 ， 而 准 粒子 是 
声 子 态 的 混合 ， 这 些 声 子 与 从 激发 的 量子 态 到 凝聚 体能 级 转移 一 
个 原子 的 过 程 有 关 . 这 样 ， 算 符 b,, bi 不 改变 总 粒子 数 ， 而 算 
IF 8 使 它 作用 的 每 一 项 的 总 粒子 数 减少 1. 要 描写 的 过 程 如 下 : 

© RNc 中 一 些 碰撞 会 转移 一 个 原子 到 Rc， 使 N -> N+1， 当 然 
也 有 相反 的 过 程 ， 即 非 凝 聚 带 的 原子 与 凝聚 带 内 的 原子 碰撞 使 一 
个 原子 从 凝聚 带 到 非 凝聚 带 ， 使 N -> N -1. @ 考 虑 最 初 没 有 凝 
聚 的 情况 时 ，Rec ALR nc 的 边界 被 固定 在 -个 凝聚 体 最 终 形成 后 
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的 适当 值 . OMER RA HEIT, Rue 中 的 化 学 势 变 为 非 负 
的 ， 这 是 允许 的 ， 只 要 化 学 势 不 超 过 R wc 的 最 低能 量 . OTI 
HEHA u, Bogoliubov 谱 和 和 波 函 数 对 所 有 的 N( 大 的 或 小 的 ) 
都 是 有 效 的 ， 央 为 对 于 小 的 N 和 ww ， 这 里 的 结果 与 微 扰 理论 并 没 
有 大 的 不 同 . 

通过 极限 过 程 定义 宏观 波 函 数 < N, mt | 6" 0) N -1,m, t> 
~ VNE DO 的 意义 在 于 “由 于 凝聚 体 包含 的 粒子 数 是 
宏观 数量 的 ， 改 变数 1 基本 上 不 影响 系统 的 态 ， 因 此 可 以 说 凝聚 
体 中 增加 ( 移 去 ) 一 个 粒子 的 结果 是 转换 N 粒子 系统 的 态 到 
N e VASE AIR AS, PES nmt» 和 | N + 1,m,t > 表示 两 
个 “相同 的 ” 态 仅 有 系统 中 的 粒子 数 相差 1 的 不 同 ” 但 如 果 考 虑 
凝聚 的 生长 ， 这 种 差别 会 有 不 可 忽略 的 影响 . 
6.2.3 主 方程 的 导出 

应 用 量子 动力 学 理论 "外 可 以 描写 凝聚 体 状 态 的 密度 算 符 的 主 
方程 . 一 个 原子 散射 进入 、 离 开 Rec 的 主要 过 程 有 如 下 六 种 情况 ， 
即 没有 准 粒 子 数 的 改变 ， 产 生 一 个 准 粒 子 和 吸收 一 个 准 粒 子 条 件 
下 的 N 一 N+1. 图 6-2 是 可 能 跃迁 的 示意 图 . 

E 


图 6-2 借助 准 粒 子 的 跃迁 示意 图 
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这 六 种 跃迁 几率 可 借助 函数 RÀ 表示 为 


W'(N)- R* (Ëy, uy / h) (6-2) 
W'(N)- R (yi uy! Rh) (6-3) 
Wa (N) = R' (fn (EX + ul A) (6-4) 
Wy (N) = R (fn: (EN-1 + ya)! B) (6-5) 
Wa (N) = R'(g,, (EN + uy) h) (6-6) 
Wa’ (N) = R (2,.(-68. + uy )/ h) (6-7) 
BRR 被 定义 为 
+ u? 
R*(y,@) = Gata [d rja rx ro) Fj + FW, (r, k) 
(6-8) 
Re) = aml? r[arX To + FU + AAW, r, k) 
(6-9) 
式 中 
dl = dikid kd? k,d?k , 
A(T,w) = ó[ ^eyys(r)(K, + k, - k; - k), 
以 及 
1 . v s TA 
Matos fa? vy (r+ Dy (re 4 (6-10) 


CRUS T UE PEE y(r) 的 Wigner 函数 .请 = F(k,r) 是 相 体 
BA ERR RL BIE. 


Awa (r) = Og, (r) + Og, (r) - Og, (r) (6-11) 
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2,2 


hoy, (r) 2 
ho, (r) = 
j 2m 


+Vea(r) (6-12) 


我 们 可 以 写 出 占有 数 几率 p(N,n) 的 随机 主 方程 ，n = (n,) 
是 所 有 准 粒子 占有 数 集 ( 矢量 ) . 


P(N, n) = 2NW* (N - Dp(N - 1) -2(N + 1W*(N)pN, n) + 
AN + DW (N + Dp(N +1, n) - 2NW" (N)p(N, n) + 
Xn Wi CN -DPN - n - ey) - 20, + DWE (NCN, n] 
Y bo, *DW, (N+ Dp(N + ln +e,,) -2n,W, (N)p(N, n)|+ 


m 


Ebo, + DW; (N -Dp(N - ln + en) -2n,W, (N)p(N, n]« 
X sw; + Dp(N + n - e) - 2(n,, + DW,  (N)p(N, n)| 


(6-13) 
其 中 en = (++ ,0,041,0,0,--°} 仅 在 相应 的 脚 标 m 处 有 不 为 零 的 值 
BR R O o) 为 向 前 和 向 后 碰撞 速率 ， 这 些 碰 撞 使 粒子 进入 
(+) 或 离开 (-) 具 有 能 量 ho 的 凝聚 体 ， 碰撞 一 定 发 生 在 凝聚 体 
的 Wigner 函数 不 为 零 的 位 置 . (6-13) 式 中 表示 跃迁 到 凝聚 体 基态 的 
几率 的 项 ， 展 示 出 按 近 似 w 的 碰撞 速率 而 受 激增 加 ， 这 个 数 N 可 
3107, 但 仅 乘 以 的 跃迁 几率 并 不 大 . 这 样 ， 作 为 初步 近似 , 我 
们 可 以 丢掉 准 粒子 项 ( 它 约 为 碰撞 项 的 IJ， 这 样 可 略 去 相应 的 大 
多 数 凝聚 过 程 . 注意 到 与 相 空 间 的 分 布 函数 FU r) 比较 ， 凝 聚 波 
函数 在 r = 0 ARBE, RTW; (r, k) ,我 们 可 以 在 大 多 数 计 算 
中 用 零 代替 +r ， 它 的 积分 给 出 k 空间 的 凝聚 体 的 几率 密度 
1 Sv(h) 上 .最 后 可 得 到 简化 的 主 方程 
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RUN) = 2NW^(N = Dp(N - D - 2(N DW * QN)p(N) + 
AN + DW (N 4+ Dp(N +1) - 2NW  (N)p(N) 


624 涨 落 与 增 率 方程 

凝聚 密度 算 符 的 定 态 解 是 具有 不 同 粒子 数 的 态 的 混合 (不 
是 量子 力学 纯 态 )， 对 每 一 个 态 ， 凝 聚 体 的 波 函 数 纳 (Cr) 即 是 
与 时 间 无 关 的 N 粒子 Gross-Pitaevskii 方程 的 解 . 略 去 各 种 涨 
落 ， 可 求 得 平均 粒子 数 的 速率 方程 为 ! :0 


(6-14) 


N -2W* (N)((1 - e A487 D) +1} (6-15) 
4 是 蒸气 的 化 学 势 ， 由 统计 力学 决定 ; MN) H N 粒子 凝聚 
波 函 数 的 化 学 势 ， 可 由 入 粒子 定 态 Gross-Pitaevskii 方程 求 
fj. 当 二 者 相等 时 ， 即 凝聚 带 Rc 和 非 凝聚 带 Rnc 的 温度 相 
等 时 达到 平衡 ， 系 统 达 到 稳 态 . 


6.3 占有 数 p 满足 的 主 方程 
为 了 对 粒子 数 表 象 中 占有 数 p 所 满足 的 主 方程 (6-14) 式 进行 
解析 求解 ， 作 如 下 近似 处 理 
W*(n-DsW'(nsW'(ens»), 
W (n-\)>W(n)=W («n»), 


(6-16) 


WH 4A-2W*(«n») B=2W («n»), Wl(6-14)sC n] 55 Jg 


p(n) = A(np(n - 1) - (n + Dp(n) + Biin + D pin +1) - np(n)]. 
(6-17) 
用 凝聚 体 潭 没 与 产生 算 子 hp: KRR (6-17) 可 写 为 
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p= A[2b* pb ~ pbb* ~ bb* p] + B[2bpb* — pb*b — b' bp]. 

(6-18) 
(6-18) 式 是 我 们 熟知 的 激光 与 热 库 相互 作用 方程 ，4, BH 
激光 场 与 热 库 之 间 交 换 光 子 的 机 制 而 定 . 现在 是 原子 激 
光 ，(6-18) 式 第 一 项 表示 的 物理 过 程 为 ， 气 态 原子 1，2 dH 
碰撞 并 转化 为 气态 原子 3 及 凝聚 态 原子 4， 转 化 率 为 

W* x eee Dies BAF BUC 


l "esp ; fd ky Jak; Jk; [5 k,9&, + kz - ky - ky) 


8h35 = aln) / PEG OF; OD  F(k, OIE, (A) , 

(6-19) 

sth, = IH jd! xe^* (x) 为 动量 空间 凝聚 体 波 函数 . 
同样 第 ei 凝聚 原子 4 与 气态 原子 3 相 磁 转化 为 气 
态 原子 1 与 2 的 转化 率 W^ ox és eds A n ds didus BI 


is ea} Quia je kı fa kz 2 fa° ks je k,6(k, +k, =k; - k4)x 


SUVs = ior) PI + Fy ONL Fco.) Fs OF, (0 
(6-20) 
6.4 主 方程 的 解 
参照 文献 [103.104] 解 微 激 光 主 方程 的 方法 ， 首 先 定义 母 
函数 
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G(z) = ELE (6-21) 
上 是 有 

Èn z" p(n- 1) = 2(z È +66), 

Èo *Dz'p9) - 6 «0G. (6-22) 

Vins n+ = Žao). 


0 


用 $e 作用 于 主 方程 (6-17)， 并 由 (6-21)，(6-22) 便 得 G(z) 


n=0 


所 满足 的 方程 


a g Hz) = z- z ak 
G + (I - z)z + B(z - D 326€) = A(z-1)G(z) (6-23) 


(z- Be" B - Az 
z = 一 一 一 一 一 一 一 ， =B-A 6-24 
” B= Az +(e Ae" (6:24) 
易 证 


ð 
Lo G+ DO - Az) yy = 0 


ii. zol ,0 


in e (B~ 4)’ e" 
Or (B-Az*(z-DAe y 
Ma] =] Oz0 ze 
€ |. so oz =I 
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G2z0  2A(B- Ay e" (1- e?) 
à (B- Az + (2 - DAe?y 


2 


C zo| = 9^ zo| m 24 -a et z 

x i - 0, PY, » ris (1-e 7) (6-25) 
满足 (6-23) 式 的 G(z) 通 解 为 

GG = Gozo) exp[[ = 全 dz, Go(l)=1 (6-26) 


事实 上 将 (6-26) 式 代 人 (6-23) 式 得 
Ga - zs B(z- 3x; - 
A (6-27) 
(AQ - 2)2 + BG - ):5 7G = Az - DG 


于 是 可 由 (6-26) 及 (6-25) 求 得 


< n >= Z LLET 4 ge iz 
2-2 


Go Oz B - Az Oz 
a,_A -c 
< >0e Vi ') 
(6-28) 
«(Ai »s«n! >-~<n>=<nn-N>4+<n>-<n>= 
e 
2 
P] +<n>-<n>- 
La m 
(6-29) 


式 中 
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en d A dz 
mo B- Az 
OGo| —— 9" G, Oy. 1 0G, Oz, 
7) =f Gt? + = 
2° La agp | Gy & Le 


< n(n- 1) > €` 


2 


“+ <on> e" 7s (l-e) 
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| 
. A > 
1- ay} 0 o acc 
(ise ) e) 


A 
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@H| ð 4 4 ao 
az? E a: B-Az B-4Azo O |- 
(6-30) 
-(- e") 
(B- a 
由 (6-28)~(6-30) 最 后 得 
<(An) >=<(An)? »o e+ 
AEB tye) cna + (6-31) 
B- 
A 2 
——(-e7)« (l-e“)? 
B-A (B - Ay 


考虑 到 A, B 均 是 <n> 的 函数 ， 故 (6-28), (6-31) 提 供 了 

计算 <n>，< (An)? > 的 迭代 式 . 即将 时 间 间 隔 (0, 
若干 个 小 间隔 ， 每 一 间隔 At n» BR < (An)? > 的 变化 可 
参照 (6-28),(6-31) 式 写 为 


TENET ETE 07.4 (220 
kA 

«(An >k+1= 

< (Any! >p ey AtB @ tM-e td)<n>+ 

By — A 
e AR SN 
5 es) (le 
k = A (B - AY 


(6-32) 
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Au = An») Bea = BA n >p). 


考虑 到 ck 2B, 一 A «0, (632): RACE atA, I 
fl «n» K «(^m >, 4: SUL e 7*7 Fle 700 进行 放大 . 在 研 
究 原子 统计 分 布 性 质 时 ， 可 通过 下 面 的 定义 将 此 放大 因子 
去 掉 ， 即 定义 约 化 的 < 元 >, 及 < (any > 为 


< Ñ >y=< n >> je Em 
2 2 XI 
«(An) >,=< (An) > le ^" (6-33) 
定义 如 下 因子 Q, KREET fn A Be io A: 
Q, =< (AR) >, / «fi», (6-34) 
6. 5 凝聚 体 原子 数 的 增 蔓 与 统计 起 伏 的 计算 
在 上 节 求 解 的 基础 上 ， 并 参照 文献 [38,102] 得 
"- 8m(GkT)? 2u(ng)! KT y(n) H(n) £ 
A=2W USOS e { it Ki iT )| (6-35) 


BI A=W" /W* = exp(u(n) - p(n )/ kT} 
式 中 "为 凝聚 体 原子 数 ，me 为 气体 与 凝聚 体 达到 平衡 时 凝聚 体 
原子 数 . puar En. e) o>, paint 


x 2mo’ 
o=——, PARTE, m= 3.84x10% kg, @=27x415 s", 


2=4.9nm， 可 算得 c =1.817x10%m. 在 此 参数 条 件 下 ， 根 据 我 
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们 的 数值 计算 (4.17) 式 ， 化 学 势 w 与 凝聚 体 粒 子 数 n 的 关系 
可 近似 表示 为 x(n) = 0.124(n + 51? P. 于 是 通过 解 方程 (6-31) 
和 (6-32) 可 求 得 当 温 度 T=2 pK, n =5x10 情形 ， 凝 聚 体 


粒子 数 <n> 及 增 率 y = > 随时 间 的 演化 ， 所 计算 的 


结果 画 在 图 6-3 和 6-4 P. 图 6-3 中 凝聚 原子 数 随时 间 的 演 
化 曲线 与 增 率 方程 (6-15) 求 得 的 结果 是 一 致 的 . 由 图 6-4 看 
出 ， 在 BEC 形成 的 初始 阶段 增 率 很 大 ， 中 间 有 一 平台 ， 到 
末了 增 率 又 下 降 到 零 . 这 平台 的 增 率 即 可 代表 凝聚 体 的 增 
率 ， 6-5 给 出 起 伏 C 随 时 间 的 变化 ， 由 初始 Poisson 4} 4 
Q = 1 急剧 升 到 超 Poisson PEQ = 133. Yen, Ht, BER 
粒子 数 ， 增 率 及 统计 起 伏 随时 间 的 演化 也 有 类 似 于 图 6-3 ~ 
图 6-5 的 变化 曲线 . 


图 6-3 凝聚 体 粒 子 数 n 随 时 间 1 的 演化 
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d log(n)/dt 


图 6-4 HKE y 随时 间 1 的 演化 


Lai 
0.0 0.1 0.2 0.3 04 0.5 


图 6-5 起 佑 随时 间 1 的 变化 
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6.6 凝聚 增 率 与 平衡 原子 数 的 关系 

值得 - 提 的 是 在 导出 (6-35) 式 时 假定 nq 个 原子 是 自由 的 
相互 作用 ， 并 未 考虑 到 碰撞 的 受阻 ， 亦 即 假定 自由 程 是 无 限 
大 的 . 若 考 虑 到 自由 程 有 限 ， 粒 子 间 的 相互 作用 实际 上 只 能 
在 自由 程 线 度 内 有 效 进行 ， 故 GP 方程 中 对 硬 球 模型 散射 有 
页 献 的 粒子 并 非 nc 的 全 部 ， 而 只 是 靠 得 很 近 的 部 分 粒子 7, ， 
设 自由 程 为 /， 则 


Hrs 
fren dr 
0 

zug SUL adr iul 


eq ~ "eq R 5 
hrar (6-36) 


mfi-a + RI (Dye!) 


3” 
n, = 61 or RH = (6n, 1n;)?. 


故 给 定 mis. nao BE AT SL XU uu. 用 元 ,代替 (6-35) 式 中 


的 no， 便 可 算出 修正 后 的 4 与 二 . 参照 文献 [81]， 平 均 自由 
Bagh) 假定 原子 密度 为 "= y=", gg 
FH d=49nm, PAn, =8 rn ~3x 10°. 根据 文献 
[39] 的 实验 数据 ,我 们 假定 相应 不 同 na 的 平衡 温度 T 有 
如 图 6-6 所 示 的 关系 ， 对 于 不 同 na 所 计算 的 增 率 如 图 6-7 
中 曲线 所 示 ，y 随 ma 的 变化 关系 与 实验 结果 EG9 大 体 上 是 
一 致 的 . 


程 /= 
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0.9 


0.8 


0.7 


T(uK) 


0.6 


0.5 
0 100 200 300 400 500 
Na (x10") 


图 6-6 BEC 形成 的 最 终 温 度 T 与 平衡 时 
凝聚 中 的 原子 数 n 之 间 的 关系 


0 100 200 300 400 500 
N,,( x10") 


E67 凝聚 生长 速率 y 与 水 衡 时 凝聚 中 
的 原子 数 nq 之 间 的 关系 
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6.7 小 结 
本 章 给 出 了 描写 BEC 的 动力 学 模型 和 主 方程 ， 用 母 函 数 
法 求 得 玻 色 凝聚 体 主 方程 的 解 ， 并 计算 了 凝聚 体形 成 的 增 率 
与 起 伏 ， 增 率 曲线 存在 一 平台 ， 布 平台 的 粒子 起 伏 则 是 超 
Poisson 分 布 的 . 初始 增 率 y 随 最 终 平衡 原子 数 ms 的 增加 而 逐 
步 减 小 ， 与 实验 结果 大 体 上 一 致 . 
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第 七 章 BEC 自发 与 受 激 牙 迁 几率 的 增强 


7.1 引言 


激光 冷却 和 势 阱 中 原子 玻 色 - 爱 因 斯 坦 凝 聚 实验 的 快速 
进展 促进 了 这 一 现象 的 研究 . 已 经 有 很 多 文章 讨论 有 关 多 玻 
色 子 处 于 同一 量子 态 的 凝聚 体 的 光学 性 质 [453544071， 例如 为 
了 通过 冷却 的 原子 系统 的 光 散 射 进行 BEC 探测 ， 几 个 作者 
在 最 近 几 年 考虑 了 这 个 问题 44 他们 讨论 了 在 低温 
7x0 和 空间 均匀 系统 的 凝聚 体 的 弱 光 散射 . 这 样 一 种 情况 物 
理 上 相应 于 无 限 大 势 阱 极限 . 这 个 理论 已 被 扩展 到 有 限 大 小 
势 阱 中 原子 的 情形 . Javanainen 考虑 了 具有 跃迁 波长 4 的 量 级 
大 小 为 a 的 光学 薄 凝 聚 体 ol. 在 他 给 出 的 这 个 框架 内 ， 散 射 
主要 发 生 在 向 前 方向 ， 散 射 截面 有 一 个 罗 仑 兹 线 型 ， 其 宽度 
比例 于 集体 自发 镶 射 ( 超 辐射 ) 速率 . 最 近 的 文章 2 指出 上 
面 的 结论 在 有 限 尺度 a 在 (1~20)4 的 范围 内 的 光学 厚 势 阱 中 
应 被 修正 . 在 适度 的 原子 密度 以 致使 A/m) x1~10， 人 
们 并 观察 不 到 背 散 射 ， 传 播 效应 以 某 种 方式 影响 散射 过 程 ， 
参考 文献 [108] 的 方法 变 得 不 适用 ， 谱 变 成 侨 罗 仑 效 型 并 在 
共振 处 展示 出 一 个 窗 峰 . 在 文献 [47] 中 You 等 研究 了 来 自 热 
阱 中 冷 玻 色 原子 系统 的 弱 光 相干 散射 ， 他 们 推导 并 详细 讨论 
了 描写 来 自 不 同 原 子 的 多 次 散射 过 程 的 散射 方程 . 他 们 提出 
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这 个 方程 的 两 种 解析 求解 方法 .一 个 是 根据 on-shell 近似 的 
散射 方法 0o21， 另 一 个 是 根据 Glauber 的 广义 衍射 9 他们 计 
算 了 散射 截面 ， 角 分 布 和 谱 分 布 有 必要 强调 的 是 他 们 之 所 
以 能 求 得 散射 方程 的 解析 解 主要 因为 他 们 应 用 了 以 下 假定 : 
OGRE KERI p(k)p(k') mieu! KE. @ 忽略 了 偏振 
Bay Eeu 的 依赖 关系 [和 后来， 讨论 了 激发 态 原 子 穿 过 凝 
聚 体 时 所 产生 的 自发 辐射 人 51. 他们 将 总 的 辐射 速率 Yo 表 
示 为 自发 和 受 激 而 进入 基态 的 辐射 速率 之 和 ye rN, 
KEN BEC 中 的 原子 数 ， 比 例 系数 /与 基态 的 Frank- 
Condon Al FAK, 它 可 在 假定 基态 和 激发 态 分 别 为 大 小 为 
h 和 /的 高 斯 波 函 数 而 求 得 . 最 近 ，Hu Huang 54:69 应 用 
类 似 的 处 理 计算 了 Frank-Condon 因子 ， 最 后 得 到 了 在 BEC 
中 吸收 共振 的 Lamb 移 位 的 修正 5 比例 于 N/ksxo,xo 是 高 
斯 分 布 的 大 小 . 本 章 中 我 们 考虑 在 低温 下 和 弱 光 相互 作用 的 
可 准确 求解 的 N 原子 辐射 模型 体系 . 我 们 导出 自发 和 受 激 辐 
射 速率 的 一 般 的 严格 公式 tn, 但 不 用 文献 [109] 中 自 能 核 的 
计算 中 使 用 的 on-shell 近似 和 两 个 假定 . 从 这 个 公式 我 们 得 
到 辐射 速率 的 增强 与 两 个 因素 有 关 . 一 是 原子 的 密度 分 布 ， 
男 一 个 则 是 德 布 罗 意 波长 a 实际 上 ， 所 求 得 的 自发 与 受 激 
辐射 增强 因子 正比 于 体积 私人 2r2) ( 4 是 原子 的 跃迁 波长 ) 
中 原子 数 »， 并 以 一 种 复杂 的 函数 关系 依赖 于 德 布 罗 意 波长 
dg BL. BA, SAT MER KA HE gy = Von dy sa 
=k, /Vmk7 = 2.45 VW, INE RAR US EI OCA. 
AT 50, Ey — cH, EPR. 男 一 方面 ，Lamb 移 位 的 修正 ， 
Hé oot, BPE. Hésolt, BE -ATBIBI SHE 


113. * 


2000 年 上 海 大 学 博士 学 位 论文 

BEC 形成 实验 的 诊断 中 是 有 用 的 0521. 很 明显 ， 我 们 所 得 到 
的 结果 非常 不 同 于 由 作者 You, Hope and Hu Huang 所 得 到 的 
结果 . 解决 问题 的 关键 在 于 我 们 严格 的 解析 解 考虑 了 原子 的 
热 运动 ， 也 没有 上 面 提 到 的 近似 和 假定 . 


7. 2 单 原 子 辐射 的 线 宽 与 能 级 移 位 

初始 处 于 激发 态 的 原子 自发 辐射 回 到 基态 ， 由 于 辐射 场 对 
原子 的 反作用 ， 会 使 得 原子 辐射 有 一 定 的 线 宽 且 有 能 级 移 位 . 这 
从 原子 波 范 数 随时 间 的 变化 可 以 看 出 来 . 这 就 是 Wigner 与 
Weisskopf 的 自然 线 宽 理论 0 了 1 

电磁 场 量 子 化 给 出 电场 E 经 量子 化 后 按 空间 模式 wx(7) 的 
RAR" 


E(r, n5 ci Y nto, (b+ e? - bi e"? E Ma(r) = 


: (7-1) 

DNE e"! - E; ge! us (r) 
式 中 4 为 模式 指标 ，z; ,为 光 的 偏振 ，y(r) 为 归 一 化 的 空间 模式 ， 
若 为 平面 波 ， 则 yx (7) = 上 2 et. Legen RS, b4,6i 分 
别 为 第 4 模式 的 光子 的 潭 灭 与 产生 算 符 . 由 上 式 给 出 的 E(r,1) 是 一 
个 算 符 展开 ， 只 有 作用 在 态 函 数 上 才 有 意义 . 例如 作用 在 真空 态 
[0; >, 5,10, >=0,b} 10; »41; >, 故 上 式 为 

E(r,1)|0 >= 
i Ly 2mher, (5; e" be YE, oui (r) > 


(7-2) 
PERL US EIOS BOE VEG Met ay go!) 
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«gl; | Hle 0; >= 


[2 , Ls 
i 2mh inu ual- EO) < elr] g> Ero = (7-3) 
à 
(-e)i 区 e'®% a (Weg Exe 
@) 


ER TTAB -DREH "A. Eun DE TIE 
向 末 态 |g > 跃迁 的 HEME TE. 参照 弱 场 微 护 理论 我 们 求 得 原子 
的 初始 状态 |e,0> PCS gz > 的 模 量 aco() ayy, CO 随时 间 
1 的 变 率 方程 为 


da,o(t) 1 " legl 
Sar Tin o Olg nem agy (O) (7-4) 
da,,,() 1 ， -i veg 
di yp S gla Hle» en asol) (7-5) 
由 此 解 出 


1 
dg; T iy fis gla lH'le,0 > e aeo(t')dt'= 


-e [2nh : 

e|2nh ua. . “Hie “at 
ee iy ELS fe Kogo S Id 
h Yo Ll ` ! 


(7-6) 


2 m : 
Veg Ezo y fie ns aM MOTTO )dr 


day o(t) TE 2nhe? ( 


2 
dt h^ gd 0; D 


(7-7) 
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tf a. (0) 作 Laplace 变换 ， 并 注意 到 4a.0(0)=1， 可 得 


1 


Cd e he (EIU 
à(s)-| s (7-8) 
cde ud PER D 
由 Laplace 反 变换 得 出 初 态 模 量 a.0(7) 
MAL 
hem e"ds (7-9) 


1 
aea (D Fy He e ait 2nhe!(V,, Eig) 
h Do (er 0; tis) 
EEE s 复 平 面 虚 轴 的 路 径 积 分 , 再 补 上 左 半 平 面 无 限 大 半圆 ( 贡 
献 可 略 去 )， 便 是 一 包括 左 半 平 面 的 闭路 积分 ， 显 然 对 积分 作出 
贡献 的 即 被 积 函 数 在 左 半 平面 的 极点 R.(s)<0， 当 s 很 小 时 


1 


Og Ca is 
Quy Oa is 
~ (7-1 


kaa PIRE 
(Weg Oa) +S (Mg 0) +s 


inmd(@y~ O 
k ( cg 0) 
0) 


式 中 了 为 取 主 值 ， 代 入 上 式 得 
e"ds ni dae 
=se“ (7-11) 


l - 
aco 7 5 Ss 
21 stzl.ti^o; 


式 中 
(7-12) 


De = aD 
2nhe (Ve Exo) 
2n e Elo) o 
y 2 : (Og - 0;) 
"m Lo; 


e 
h 
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- oru) 
^ 1 2nhe" (v, Eig) P 
[ES 


(7-13) 


2 3 - 
h gdje Lo; OU, — 0, 


7.3 v 原子 的 自发 跃迁 速率 
将 单 原子 的 Wigner-Weisskopf 自然 线 宽 理论 [031 


1 i * e0|H'|g,l; >< gl, |H'le,0> 
a is hag aey REGE! gla iHi Z 
2 h? io Ou Wi ids 


I 
=i gio 


Ache;(Gu, @aris) ia”? Og tis 


2,— 一 2 
Rne? (Vo Ero)? 


(7-14) 
推广 到 入 原子 情形 , 定 态 波 函数 多 , 取 如 下 形式 : 
Pn = Vn ns May s 
Wm =e EMM Pm Pe (Gi) (7-15) 
REPER GA WARP, G) Sel! sy ER i 个 原子 的 
内 部 及 质心 坐标 的 波 函 数 . 可 以 直接 将 (7-14) 式 中 的 厅 代 以 
Sy，， 这 样 便 得 到 : 


i=l 
l +i 
zyti^eo,- 
37 


i Me la ><g,li|H,'le,0> 
K iaija Og @,+is 


(1-16) 
应 用 (7-14), (7-15) 式 便 得 出 相互 作用 密度 矩阵 元 为 


a 1 K 


2000 年 上 海 大 学 博士 学 位 论文 


ika Ì 


<e0{H' [gla >= hgior (Fe 


. BT 
gi pea pw Pro» (47) 
100; 


r'()-e ie Pg) Rih 


i@gitik, Fy 


<glalH'|e0>=hgrcr (7))e 


29: = 
Rig =ie hob” Eio, (7-18) 


È (=e Bg) rj/h 
将 (7-16) 式 写 为 
T _ .2 «rn (7)> 
gr tito “idee m 01 is T 
p - (7-19) 
$E Or (0C)? i qz 
jA Deg Witis e 
注意 到 方 括号 中 的 第 一 项 表示 自作 用 ， 而 第 二 项 表示 N 原 
子 之 间 的 相互 作用 ; n.r 分 别 是 能 级 的 “上 ”和 “下 ” 算 
TR. «el|r'Gr (i)le>=1，<glr (i)r (5)1g>=0. 这 样 


Ser" Gaye 6)» 便 是 个 原子 中 处 于 激发 态 的 原子 数 ， 而 方 
程 (7-19) 恰 是 这 些 原子 的 总 的 线 宽 与 位 移 。 车 我 们 将 (7-19) 式 两 
WH aer Gr E> 规 一 化 ， 并 将 规 一 化 后 的 了 y+iAw 用 


同样 的 符号 表示 
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1 
Dy tide =i Ygi,———————* 


jo” Qu @, +is 


iD Days Gr (F + 25)» n) 
sr Gr G)> (7-20) 


m 1 ki 5j 
LAFI) = Ygi, 一 一 一 一 一 。 Hà fh 
E 


Dy Witis 


ar; = Ñ 


式 中 第 一 项 丛 是 单个 原子 的 线 宽 与 位 移 了 y +iAo ， 第 二 项 为 
NN 个 原子 带 来 的 修正 . 函数 上 (AP ) 可 通过 下 式 计算 !o 


ba 8) iW(G,0),£- kz dr, 
K(5,9)-— iin? e = +( 3cos? 6 x ol 
WOE) => sin? os +(1 3cos? 6" et + £088 


= 


(7-21) 
式 中 |k 2/2, AREER, 0 EBRE i M Ar, 


之 间 的 夹 角 . 现在 的 主要 问题 是 如 何 计算 关联 函数 
«r' yr GN). 由 (7-17),(7-18) 式 得 出 


<r’ (Fr (F + AF) >= cro», 


« r* (rn )r (FR) >< e Pe Px Ari th E 


因此 , (7-20) 式 写 为 如 下 形式 : 
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it +iAw = 
iT +ið@, +} < e Pe BRA (A >L(AF,) 
AF, 
(7223) 
用 ADe ^p, 分 别 表示 Peo Pe 的 热 涨 落 ， 那 么 有 p. =< p, > +Ap,, 
Py =< p, > +, WRI FI Hk, 方向 人 射 ， 由 于 动量 守恒 ， 
WAT « p. »-«D, »-hk,. 然而 ， 如 果 光 子 沿 两 个 相反 的 方向 人 
射 ， 有 时 吸收 右 光子 有 时 吸收 左 光子 ， 平 均 吸 收 的 净 光 子 动量 
是 零 ，<P.>-< Pe。>=0. 因此 ， 我 们 称 前 者 为 单 边 问题 ， 相 
应 于 通过 原子 的 散射 光 ; 而 后 者 为 双边 问题 ， 相 应 于 原子 各 向 
同性 的 自发 辐射 Cir gt ERE SHEE). 下 面 首先 讨论 
自发 辐射 情况 ， 下 列 方程 成 立 : 
< e PePe s/h gi Pe Mg AFM 
" i Ap, AF, NC " 
h 2n? 
ip, Ar, (App: Ar)? 
h oR? 


«I <l- 


>z (7-24) 
loan -H apg aF,/M? 
E2 c2 £ 


ELTE 
e pen 


ROP Ag = V2xh?/mKT7 ， 在 上 式 的 推导 中 利用 了 
^p, [2m = ^p, /2m = KT 12. 进一步 可 假定 在 体积 如 /2x” 中 
有 n= p4n(A/22)) 个 原子 则 
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£< e UP wiih >L( AF, ) = 
My 


2n( MF)? ag 


fod ar, LOAF, )e 
y [£^ dé [sim doe" [K(E,0)-iWE,A)]= (1-25) 
0 0 


2m(y'+id"), 


£p 7 214g /À 
Ln MT -了 tito, + nid") (7-26) 


方程 (7-26) 给 出 N 原子 自发 辐射 跃迁 几率 为 几 ，N 个 原子 自发 
NORTEL sino 恰 是 单 原子 的 二 ye tiso, WME 
增强 项 2y(y'+i6') ,对 于 空间 分 布 不 均匀 的 大 小 为 a 的 高 斯 分 
fi, nance! !” ， 上 面 的 方程 可 以 推广 为 


Y< eU Dg) ^rj!h >L(AF,) = 
AF; 


mo JE dë Jsin0d0 x (7-27) 


e PIR P wee) iwquey- 
2yno(7'+ id") 
对 于 单 边 问题 , (7.24)、(7-25) 式 可 类 似 地 写 为 
« e Pec erah y 


(1-28) 


cel Pe Peraih 、 e'f! M, 
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yen Anh >LCAF ) = 


fodar, LCS )e OP n A 
rpe ag finado zy [[sino, 49, ag, x 
eU HEMT KE 9) iW EO) = (7-29) 
mo j£? a£ fsino aoje FE x 
0 0 
(K(5.0) -iW(,0)] - 
2jymo(y'*ió*) 
这 里 jo( ) 表示 球 贝 塞 尔 函 数 .在 (7-29) 式 的 推导 中 使 用 了 下 面 的 等 式 


fsinodo— f fsing, do, d$, e " 9? P, (cose) = 
D sa l (7-30) 


i" j, (Ky, Ar; ) Pp (COS0)So,m 

如 图 7-1 "PAR ARMY AE 0,6,.0,0, 1 d i EEA cos") = 
cos(8,) cos(8) + sin(@, )sin@ cos(d, — 6). 比较 (7-27) 和 (7-29) 式 可 
以 看 出 ， 两 个 方程 都 是 通过 高 斯 函数 e EE eeu, 和 
a (为 简单 ,下面 仍 用 1 号 表示 Tsa + 1/1a”)， 其 差别 仅 在 于 
(7-29) 式 中 多 了 球 贝 塞 尔 函 数 JO. WWE, REET y's Bi 
如 的 变化 关系 曲线 ， 图 7-2 表示 (7-27) 式 各 向 同性 的 自发 辐射 情 
形 ，(7-29) 式 的 单 边 散射 情形 画 在 图 7-3 中 . 
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Hi(0 .1 ) 


图 7-1 矢量 kj .Ar 和 4 的 示意 图 


图 7-2 根据 方程 (7-14) 画 出 的 因子 y'. OB Ey 的 变化 关系 
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图 7-3 根据 方程 (7-16) 面 出 的 因子 p. Bü Ey 的 变化 关系 . 
7.4 NN 原子 的 受 激 跃迁 速率 
如 果 考 虑 原子 的 受 激 辐射 ， 场 算 符 作 用 于 光子 态 | n; > 得 到 
(b; e? — b; e!!! )|n; >= 
[natlei In; &1»— (7-31) 
dna 71e! |n; -1> 


初始 和 中 间 态 分 别 是 lema > [ging dos | gin, -1>. 类 似 演 
算 导 致 由 于 自发 辐射 和 受 激 辐射 (xm ) 引起 的 单 原子 的 线 宽 
和 能 级 移动 


1 
=yt+iha, = 
3” s 


: 2 n, +1 n; 
id 8ic( ad " E 
xe Qu,—0; HİS (Qu +O; tis 


(7-32) 
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由 (7-24) 式 并 参考 (7-20) 式 的 推导 ， 可 以 写 出 相应 于 N 原子 的 公式 


D SPA = 
2 


ntl n; 
+ 


~~ 
i Vein : 一 ) = 
de Weg— Qj +18 Ou +O, +is 


i2 Laos Gir (i + AF))> 
i = L(Ar,) 
Ver Gr c) i 


n+l x n; e Nu 


LAF) = Dei, ( 


Ae Og — OZ lS Woe ta, tis 
g E 


2 n; +l ikz dF, 
in. Qeg- 0, +i : 
(7-33) 
在 热平衡 状态 下 ， 按 辐射 系数 41,, By, By, 的 爱 因 斯 坦 关系 ， 
可 写 下 受 激 辐射 的 跃迁 几率 与 吸收 几率 分 别 为 : 


3 
ho; 
à 
By Nz = Ayn, Ny = 
ne? 


1 
Sr +i do, + yn(y*ió')]n; Ny 


i (1-34) 
By na N, = Ann N= 
ne 
ai, +i Aw, * n(g*ió')]n; N, 
No Ny 分 别 为 基态 、 激 发 态 的 原子 数 、 n; 是 模式 4 的 平均 


ot TC 
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7.5 数值 结果 及 与 文献 结果 的 比较 


7.5.1 数值 结果 

参照 各 向 同性 的 自发 辐射 情况 的 图 7-2， 增 强 因子 7 在 
£y 7245 有 极 大 值 y'=1.453， HUS ELIT ACH RECS GE, 
并 且 能 级 移 位 因子 5' 趋 于 -1.0. 将 典型 的 钠 原 子 凝 聚 体 的 数 
HE"), m=23x1.67x10g, 4-58«m, p - 10 cm?, T22 t KAKA 


En = xà 14 = 1n 1A ,n= pP Qn), ATO 
mKT 


4g 20434 «245, n=1.035, Ak, 2y 为 单位 的 增强 因子 ， 
n(y'+id") 20.0355 + i0.1136. 然而 ， 若 将 温度 降 至 0.06274K, 
IEH} čp 22.45, 增强 因子 达 于 极 大 值 n7 i0") = 1.505 +10.711. 
图 7-3 给 出 了 (7-29) 式 单 边 光 散 射 情形 的 类 似 结果 ， 当 后 较 大 
时 ， 辐 射 速率 的 增强 y' 随 és 几乎 线性 地 增加 ， 位 移 修正 5' 减 小 
到 较 小 的 值 . 


7.5.2 与 You 的 解 (站 的 关系 

参照 文献 [48]， 我 们 忽略 原子 的 热 运 动 ，Ap。 = ^p, = 0 ， 
并 采用 onshell 近似 ， 对 单 边 问题 fk, 1=|k， |. m 
(ar, siso) 表示 的 方程 (7-20) 的 第 二 项 ， 简 化 为 


1 iDa ee cd Fyre pepe 
(insi) = = 

bi S eR" Gor (F)> 

Y goa 


j———2 7*7 — [pr Xr, d AF, x 
C, g0, 十 is 
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ffdcoso, dé, err i 
Gr +i Ao.) fl dF, SF) d AF, x 
[[dcoso do, Xi" Qn n» 
Jn (ki Ar, )P, (cos 0,) x 
X CD" (2m + Djs (ka Ar; ) Pn (cosy) = 
1 2 ~ 
(Gr tide, yAn J p( AF, X A7, ) d AF, x 
YJ2n + Dj (ky, Or; )P, (cos) 
(7-35) 
其 中 cosy = cos cos, *sin8sin6, cos(ó - ó,). 应 用 类 似 于 You 
论文 9 中 的 记 法 
1 1 ; VEN 
gin = Ani y +i No) oC Nr Yr, Y d AF, ja Cki v) 
(7-36) 
这 样 ，N 原子 散射 辐射 的 增强 (Ar +iAw Ja 可 被 表示 为 


27 十 1 
2 


1 
—AT, +i do )。 = 
G otio )w =>, 


Y, P, (cos 0) (7-37) 


这 个 方程 等 价 于 You 论文 [9 中 的 (510), (54), (56) 和 (57) 式 ， 从 
方程 (7-36) 和 (7-37) 的 推导 ， 我 们 得 出 You 的 关于 光 散 射 的 论 
文 属于 单 边 问题 ， 除了 on-shell 近似 和 两 个 前 面 提 到 的 假定 
以 外 " ,得 到 上 述 散射 截面 增强 的 最 重要 的 原因 是 忽略 了 原 
子 的 热 运动 , 而 我 们 的 属于 双边 问题 的 自发 辐射 解 (7-27) 包 
括 原子 的 热 运动 ， 也 不 需要 on-shell 近似 和 两 个 假定 . 
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7.6 小 结 


本 章 给 出 一 准确 求解 N 原子 辐射 几率 的 模型 ， 所 求 得 的 自 
发 与 受 激 跃迁 几率 的 增强 因子 正比 于 用 2x?)(4 为 原子 的 跃迁 
波长 ) 体 积 内 的 原子 数 n ， 并 依赖 于 德 布 罗 意 波长 ta. BIDS AS 5j 
原子 的 跃迁 波长 4 之 比 &p =V2n Àg/ A = Ak, /VmkT =2.45 
时 , 辐射 速率 的 增强 因子 5 达到 它 的 极 大 值 , AT 0, Ey >o 
时 ， 增 强 因子 5' 趋 于 零 . 另 一 方面 ， 当 上 -0 时 ，Lamb 移 位 的 
BES BTE; 当 & >w% 时 ，#5' 趋 于 一 个 有 限 值 . 
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第 八 章 ”结论 和 展望 


本 文 简要 回顾 了 BEC 研究 的 历史 ,介绍 了 实现 BEC 
的 相关 技术 和 碱 金 属 原子 和 和 氨 原子 BEC 的 主要 实验 结果 和 
JLE BEC 的 理论 研究 工作 ， 较 详细 地 介绍 和 总 结 了 理想 
气体 的 玻 色 - 爱 因 斯 坦 凝 聚 的 性 质 和 规律 . 这 些 都 是 我 们 开 
展 非 理 想 气 体 BEC 性 质 理论 研究 的 基础 . 

我 们 计算 了 一 维 箱 势 中 有 排斥 和 吸引 相互 作用 中 性 原 
子 体系 的 波 函 数 和 能 量 的 本 征 解析 解 ， 给 出 了 解 有 相互 作 
用 情形 中 性 原子 玻 色 - 爱 因 斯 坦 凝聚 的 一 般 方法 ， 估 算 了 
这 种 非 理 想 气体 BEC 的 临界 温度 随 原子 间 相 互 作用 大 小 和 
原子 密度 的 变化 情况 ， 得 出 随 着 凝聚 原子 数 的 增加 或 原子 
加 相互 作用 强度 的 增 大 ,临界 温度 T, RS 

根据 波 函 数 归 一 化 和 收敛 条 件 ， 给 出 了 非 线性 定 态 薛 
定 请 方程 数值 来 解 的 一 种 方法 ， 并 求解 了 简 谐 势 阱 中 具有 


和 激发 态 解 .由 三 维 球 方 势 和 简 谐 势 阱 中 具有 吸引 相互 作 
用 原子 体系 NLSE 的 数值 本 征 解 ， 得 到 类 似 的 特性 ， 即 原 
子 的 空间 分 布 随 相 互 作用 的 增强 而 被 连续 压缩 ， 以 及 凝聚 
原子 数 随 能 量 本 征 值 变化 的 双 稳 态 曲线 ， 由 此 得 到 与 实验 
报道 相符 的 吸引 BEC 所 能 包含 的 最 大 原子 数 Nu. 对 每 一 个 
N < WN. 存 在 一 个 亚 稳 凝聚 态 和 相应 同一 原子 数 的 不 稳 的 重 
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WAS. 在 临界 点 附近 亚 稳 凝 聚 可 以 很 容易 通过 宏观 量子 隧 
道 而 进入 塌 缩 态 ; 然而 ， 远 离 临 界 点 时 ， 凝聚 体 应 是 相当 
稳定 的 . 

根据 最 近 的 BEC 的 动力 学 模型 ， 用 母 函 数 方法 求 得 玻 
色 凝 聚 体 主 方程 的 解 ， 并 计算 了 凝聚 体形 成 的 增 率 与 起 伏 ， 
增 率 曲线 存在 一 平台 ， 而 平台 的 粒子 数 起 伏 则 是 超 泊 松 分 
布 的 , 初始 增 率 y 随 最 终 平衡 时 凝聚 原子 数 ns 的 增加 而 逐 
步 减 小 ， 与 实验 结果 大 体 上 一 致 . 

最 后 ， 给 出 了 BEC 原子 辐射 几率 准确 求解 模型 ， 求 得 
玻 色 - 爱 因 斯 坦 凝 聚 体 的 自发 与 受 激 跃迁 几率 的 增强 因子 正 
HLF A QT). (4 为 原子 的 跃迁 波长 ) 体 积 内 的 原子 数 n， 并 
依赖 于 德 布 罗 意 波长 4. 

关于 EBC 形成 的 动力 学 、 原 子 光学 、 非 线性 原子 光学 
等 方面 的 研究 方兴未艾 . 在 原子 光学 方面 ,凝聚 体 作为 相干 
物质 波 的 干涉 、 衍 射 效应 ,光波 对 物质 波 的 散射 等 都 需要 
做 进一步 的 实验 和 理论 研究 . 原子 激光 的 相干 特性 的 研究 
指出 ,宏观 数量 的 原子 处 于 同一 量子 态 表明 凝聚 体 是 高 度 
相干 的 ， 但 它们 处 于 捕 陷 状态 ， 凝 聚 体 的 相位 和 振幅 都 是 
很 简单 的 ， 不 像 光学 腔 中 复杂 的 相干 性 质 . 另 一 个 就 是 耦 
合 输出 后 波 的 相干 性 质 ， 由 于 光子 之 间 的 相互 作用 可 以 忽 
略 ， 光 波 近似 在 自由 空间 中 传播 . 原子 激光 情况 就 十 分 不 同 
了 ， 因 为 原子 间 有 不 可 忽略 的 原子 间 相 互 作用 ， 耦 合 输出 的 
凝聚 体 的 相干 性 质 就 非常 不 同 于 捕 陷 凝聚 体 的 特性 ， 特 别 
是 它们 的 高 阶 相干 函数 . 文献 [57] 根 据 双 光子 Raman. 散射 
技术 (一 级 Bragg 衍射 ) 产 生 了 高 度 准 直 的 原子 激光 . 在 此 之 

- 后， 物质 波 放大 装置 的 研究 就 是 一 个 重要 课题 . 物质 波 干涉 
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仪 是 获取 高 相干 效率 的 装置 ， 也 是 获知 凝聚 体 振 幅 和 相位 
信息 的 一 种 理想 工具 . 还 有 一 种 物质 波 的 相位 相干 放大 装 
置 ， 其 过 程 依 赖 于 物质 波 的 超 辐 射 效应 . 这 是 自发 辐射 过 
程 导致 自 放 大 的 相干 过 程 . 原子 激光 的 实现 将 进一步 拓展 
物理 学 的 一 个 新 的 研究 领域 一 一 非 线 性 原子 光学 . 由 于 物质 
波 满足 的 GP 方程 中 的 非 线 性 项 有 类 似 于 Kerr 型 相互 作用 
项 ， 物质 波 也 应 有 相应 于 非 线 性 光学 中 的 三 阶 非 线性 效应 
被 观察 到 ， 如 物质 波 的 多 波 混 频 等 . 在 一 定 条 件 下 由 GP 方 
程 可 以 推出 孤立 物质 波 的 存在 ，GP 方程 也 支配 着 孤立 波 的 
传播 ， 排 斥 的 原子 间 相 互 作用 仅 维 持 暗 孤子 的 传播 . 高 度 
准 直 、 准 连续 的 可 调 原子 激光 为 用 强 物质 波源 进行 实验 提 
供 了 可 能 性 . 另外 关于 EBC 的 形成 的 动力 学 特性 还 不 能 很 
好 地 被 理解 ， 特 别 是 凝聚 的 初始 生长 阶段 ; 连续 高 亮度 的 
原子 激光 实现 的 理论 和 实验 技术 ;，EBC 的 激发 以 及 依据 平 
均 场 的 各 种 非 线性 效应 的 研究 . 将 来 更 多 的 宏观 量子 态 特 
性 的 实验 必 将 推动 玻 色 系统 多 体 理 论 ， 特 别 是 非 均匀 玻 色 
系统 理论 的 发 展 ; 同时 理论 和 实验 研究 都 将 会 为 这 种 新 物 
态 开 拓 新 的 应 用 领域 . 
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